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HAUPTAUFSÄTZE 


Turbulenz bei Wärmeschichtung ). 
Von H. Schliehting VDI ın Friedriehshafen a. B. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung, Göttingen.) 


$ 1. Einleitung. An kühlen Sommerabenden kann man bei schwachem Wind zuweilen 
die Beobachtunz machen, daß über feuchten Wiesen Nebelschwaden mit einer sehr scharfen 
Begrenzung dahinziehen, ein Zeichen dafür, daß die Turbulenz des Windes gänzlich aufgehört 
hat. so daß die Luftschiehten ohne turbulente Vermischung laminar übereinander hingleiten, 
Dies ist dadurch zu erklären, daß sich infolge der abendlichen Abkühlung ein ausgeprägtes 
Temperaturgefälle gebildet hat, welches eine Vermischung der wärmeren und infolgedessen 
spezifisch leichteren oberen Schichten mit den bodennahen, kälteren, schwereren Luftschichten 
verhindert. 

Hiermit in engem Zusammenhang stehen Beobachtungen von L. F. Ricehardson[l]?) 
über die Abhängigkeit der Böigkeit des Windes, die ja als ein Maß für die Turbulenzstärke 
angesehen werden kann, vom vertikalen Temperaturgefälle. Er stellte fest, daß die Böigkeit 
sering ist, wenn die Luft unten kälter ist als oben, und daß sie zunimmt, wenn sich das 
Temperaturgefälle dem labilen adiabatischen Gleichgewicht nähert. Überraschenderweise 
wird aber die Böigkeit nicht mehr wesentlich stärker, wenn das instabile überadiabatische 
Temperaturgefälle vorliegt. 

Ahnliche Beobachtungen hat schon 1916 G. 1. Taylor[2] gemacht. Durch gleichzeitige Auf- 
zeichnungen von Temperatur und Windgeschwindigkeit in etwa 40 m Höhe über dem Boden 
während mehrerer Tage und Nächte stellte er fest, daß in denjenigen Nächten, wo infolge 
starker Ausstrahlung das nächtliche Temperaturminimum sehr ausgeprägt war, die Turbulenz 
sänzlich erloschen war. In Nächten mit weniger starker Abkühlung dagegen, etwa infolge 
von Bewölkung, waren die Windgeschwindigkeitsschwankungen fast ebenso groß wie am Tage. 

Das Dahinströmen von Süßwasser über Salzwasser ohne wesentliche Vermischung, das 
man im Meer zuweilen, z. B. im Kattegat, beobachtet, gehört in dieselbe Erscheinungsgruppe. 
Auch die Auffällige Stabilität der Bjerknesschen Polarfrontflächen, bei der «die kalten Luft- 
massen einen Keil unter den warmen bilden, ist auf die Schiehtung zurückzuführen [3]. 
(4.1. Taylor hat übrigens 1927 einen einfachen Demonstrationsversuch angegeben, dureh den 
man mit Hilfe einer Salzlösung den stabilisierenden Einfluß der Diehteschiehtung auf die 
TFurbulenz zeigen kann [#]. 


I, Vorgetragen auf dem IV. Intern. Mechanik-Kongreß, Cambridge (England), 3. 7. bis 9. 7. 34. 
2) Die Zahlen in | | verweisen auf das Literaturverzeichnis am Schluß der Arbeit. 
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Umgekehrt hat man infolge der Konvektionsbewegungen eine vermehrte Turbulenz und 
damit eine stärkere Durehmisechung, wenn infolge starker Einstrahlung die unteren Luft- 
schiehten wesentlich wärmer sind als die oberen. Beobachtungen darüber finden sich in dem 
Buch von W. Schmidt: „Der Massenaustausch in freier Luft[5].“ 

In Göttingen werden seit 1927 unter der Leitung von Professor Prandtl von H. Reichardt 
experimentelle Untersuchungen einer geschichteten Strömung ausgeführt. Ein Luftstrom wird 
in einem wagereehten Kanal zwischen einer von Wasser gekühlten und einer mit Dampf 
eceheizten Platte durehgzeblasen, wobei sowohl die obere wie die untere Platte geheizt werden 
kann (stabile bzw. instabile Schichtung). Professor Prandtl hat 1927, als diese Ver- 
suche sich noch ganz im Anfangsstadium befanden, in Form einer Energiebetrachtung eine 
sehr einfache Theorie dieser Erscheinung gegeben, die das Wesentliche wenigstens qualitativ 
klar hervortreten läßt [6]. Da die hierbei benutzten Überlegungen grundlegend für das Folgende 
sind, geben wir diese Theorie kurz wieder. 

Hs wird eine in warereehter Riehtung strömende Flüssiekeit angenommen, die eine 
Diehteschiehtung in der senkrechten Richtung (‚-Riehtung) besitzt, so daß die Dichte von 
unten nach oben stetie abnimmt. Bei den turbulenten Mischbewegrungen wird dann Arbeit 
dadurch geleistet, daß Schwereres gehoben und Leichteres gegen den Auftrieb gesenkt wird. 
Der Weg, den hierbei ein Teilchen in der senkrechten Richtung zurücklegt, bevor es sich 
mit der neuen Umgebung wieder vermischt, ist der Prandtlsche Mischungsweg ![7], und die 
Auftriebsdifferenz pro Volumeneinheit für ein um den Weg I in der Senkrechten verschobenes 

do 


Teilehen ist gl Die Arbeitsleistung pro Volumeneinheit der verschobenen Flüssigkeits- 


Ay 
masse während der Zurücklegung des Weges 1 ist demnach: 
„ Yyıt! 

d 0 4 ” d 0 
sw Way ray 

y=ıYı 
Um die am Austausch beteiligten Mengen zu erfassen, denken wir eine wagerechte Fläche F' 
auf einem Bruchteil 5, dieser Fläche herrsche Aufwärtsbewegung mit der Geschwindigkeit v, 
und auf einem Bruchteil 9, Abwärtsbewegung mit der Geschwindigkeit r,, so daß die gesamte 
pro Zeiteinheit hindurchtretende Menge Fi(P,r,—+Psr,) Ist. Damit hat man dann für die 
Hebungsleistunge infolge der Diehteschiehtung 
| ,.d OÖ 


L, F(ß,v, +ß,0)5 gl’ 


1 1 / ) dy' (dl). 


Diese Leistung muß bestritten werden aus dem Energievorrat der turbulenten Mischbewegungen. 
Dieser wird erhalten dureh die Arbeitsleistung der Grundströmung an unserem Element, die 
durch das Produkt der scheinbaren turbulenten Schubspannung mit der Verschiebungs- 
geschwindigkeit gegeben wird. Für einen Körper von der Grundfläche F und der Höhe I, in 
du 
dy 
Dabei ist 7 die turbulente Sehubspannung. für die man nach Reynolds r= ou’ v hat, wo u, v’ 
die turbulenten Schwankungsgeschwindigkeiten sind (u = ul) = mittlere Strömungsgesch windig- 
keit, Nach Prandtl ist 15. und nach dem Vorherigen *’ —=p,r,+P,"r,. Damit 
hat man bis auf einen noch zu ermittelnden Zahlenfaktor x, in den der Korrelations- 
koeftizient von #' »’ eingeht, für die turbulente Schubspannung 


dem gerade die oben angegebene Menge ausgetauscht wird, beträgt diese Leistung Fri 


RT 1 \ 
—=xo0 day: vı, TPet)), 
und damit für die Arbeitsleistung 7; der turbulenten Scheinreibung 
> 
Ir or |, UF BUN HE nn ae (2). 


Ob nun infolge der stabilen Schiehtunge die Turbulenz erlischt, oder ob sie erhalten bleibt, 
hängt davon ab, welcher von beiden Energiebeträgen der größere ist. Ist die Arbeitsleistung 
der turbulenten Scheinreibung die größere (1; > L.), dann bleibt die Differenz zur Erhaltung 
der Turbulenz übrig. Ist aber die Arbeitsleistung gegen die Gewichtsdifferenzen die größere 
(1, > br), dann muß die Turbulenz erlösehen. "Turbulenz ist also energetisch möglich, wenn 


LETER 7 do 

Fxol’ B,v, +B,.0v.)> — Flid,v, + B.v)> #P-—, 

® lan) Er 202 Sale aD dy 
oder nach Kürzung du? ydo 
zo| ? ) ‘ 
"\dy 2 dy 
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Für die dimensionslose Größe 

ydo 

o dı 

nn. Seeds re © 
me 

dy 


die für alle Strömungen mit Diehteschiehtung grundlegend ist, möge die Bezeichnung 
‚Riehardsonseche Zahl“ eingeführt werden, da L. F. Richardson sieh als erster schon 1920 
mit geschichteten Strömungen beschäftigt hat [9]. Dabei bedeutet also 9=0 die homogene 
Flüssigkeit. Nach der Prandtlschen Energieabschätzung hat man dann also 


Turbulenz energetisch möglıch für )<2z. 
Turbulenz energetisch unmöglich für >2x. 


Über den Zahlenfaktor z läßt sich etwas Genaues nicht aussagen, doch ist zu vernniten, daß 
er nahe an I liegt. Prandtl wählt z=1 und erhält damit als Stabilitätserenze 9 2. 
G.1. Taylor hat später gezeigt| 10], daß bei Verfeinerung der Prandtlschen Überleeuneen 
noeh ein Faktor 2 herausfällt, so daß man dann nach Prandtl 9 = 1 als Stabilitätserenze 
erhalten würde in Übereinstimmung mit ähnlichen Überlegungen von L. F. Rieha rdson|®). 
Später haben G.1. Taylor[il] und S. Goldstein[12] das Problem einer Strömung mit Dichte- 
schiehtung theoretisch weiter verfolgt, und zwar haben sie nach der Methode der kleinen 
Schwingungen die Stabilität einer ebenen Laminarströmung mit stabiler Dichteschiechtung 
untersueht. Zähigkeit und Kompressibilität wurden dabei vernachlässigt, und ferner beschränkten 
sie sich aus Gründen der mathematischen Einfachheit auf solche Geschwindigkeitsprofile, die 
sich aus Geradenstücken zusammensetzen. Ferner wird in der Störungsdifferentialeleiehung 
nur der Einfluß der Diehteschichtung auf die potentielle Energie (Schwerewirkung) in Betracht 
vezogen, während der Trägheitseinfluß der Diehteschiehtung vernachlässigt wird. In welchen 
Fällen dies erlaubt ist, und welche Vereinfachung es für die Reehnung bringt, werden wir 
später noch sehen. Die Resultate der Stabilitätsuntersuchung von Taylor und Goldstein 
lassen sich in den beiden Dimensionslosen © und //d darstellen. wo / die Wellenlänge der 
Störung und Ö eine charakteristische Länge des Geschwindigkeitsprofiles (Grenzschiehtdicke) 
bedeutet. Für die meisten von Tavlor und Goldstein untersuchten Fälle ergab sich 
keine eigentliche Stabilitätsgrenze für 9, sondern es existiert für jeden Wert von © immer 
noch ein gewisser, wenn auch sehr enger Bereich von instabilen Störuneswellenlängen /. 

Tavlor konnte ein definitives Resultat nur erhalten für den Fall einer nach oben oder 
unten unendlich auseedehnten Flüssiekeit mit linearer Geschwindiekeitsverteilunz und stetieer 
Diehteverteilung, wo sich 9='/, als Stabilitätsgrenze ergab. Die gleiche Stabilitätserenze 
erhielt Goldstein für den Fall einer nach oben und unten unendlich ausgedehnten Flüssig- 
keit mit stetirer Diehteverteilung, wobei die Geschwindiekeit unten und oben konstant ist, 
während sie in einer Zwischenschieht sich linear mit der Höhe ändert, 

Ein Vergleich dieses theoretischen Resultates mit Versuchsergebnissen des Göttinger 
Warmkaltluftkanals ereab durchaus keine Übereinstimmung in bezug auf die Stabilitätserenze 
9 '/,, wobei allerdings zu berücksichtigen ıst, daß die im Kanal gemessenen Geschwindig- 
keitsprofile nieht linear sind, wie in der Taylor-Goldsteinschen Theorie vorausgesetzt. 
Es erscheint deshalb angebracht, die Modifikation der Tavlor-Goldsteinschen Reehnungen 
dureh Berücksichtigung der Reibung zu untersuchen und auch für die Geschwindiekeits- 
verteilung der Laminarströmung ein solehes Profil zu wählen, das den experimentell zu 
verwirklichenden Verhältnissen besser angepaßt ist. Dies schien noch um so mehr Aussicht 
auf Erfole zu haben, als auch bei der homozenen Flüssigkeit die Stabilitätsuntersuchung erst 
dann zu einem befriedigenden Resultat geführt batte, als man die Reibung in richtiger Weise 
berücksichtigte. Dabei zeigte insbesondere W. Tollmien[13], daß die Annäherung eines Laminar- 
profile dureh Geradenstücke nicht ausreicht, sondern daß man Profile mit endlicher, von Null 
verschiedener Krümmung zugrunde zu legen hat. 

In den ersten Stabilitätsuntersuechungen der Laminarströmung einer homogenen Flüssiekeit 
von Lord Rayleigh[i4] wurde sowohl die Reibung wie die Profilkrümmung vernachlässigt, 
sie ergaben nicht die gesuchte Instabilität. Weitere Untersuchungen von A. Sommerfeld[15], 
k.v. Mises[16] und L. Hopf[17] über die Couette- Strömung (lineare Geschwindiekeits- 
verteilung), die teils die Reibung berücksichtigten, hatten ebenfalls nieht den gewünschten 
Krfole: es ergab sich Stabilität für alle Störungswellenlängen und alle Revynoldsschen 
Zahlen. Spätere Untersuchungen von Prandtl[1S] und ©. Tietjens[19], welche die größten 
Reibungselieder berücksiehtigten und Profile zugrunde legten, die sich aus Geradenstücken 
zusammensetzen, ergaben zum ersten Male eine Instabilität, aber noch keine Stabilitätsgrenze. 
kın voller Erfolg, nämlich die wirkliche theoretische Berechnung der Stabilitätsgrenze (kritische 

21° 
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Revnoldssche Zahl) war aber diesen Untersuchungen erst beschert, als Tollmien auch 
noch die Krümmung des Geschwindigkeitsprofils berücksichtigte. 

Die Tavlor-Goldsteinschen Stabilitätsuntersuchungen der Strömung einer geschichteten 
Flüssigkeit bilden das Analogon zu den Rayleighschen Untersuchungen der homogenen 
Flüssiekeit, da sie sowohl Reibung wie Profilkrümmung vernachlässigen. Wir wollen sie im 
folgenden dadureh erweitern, daß wir Reibung und Profilkrümmung in Betracht ziehen (vgl. 
das nachstehende Schema). 


Schema der bisher ausgeführten Stabilitätsuntersuchungen zur Turbulenzentstehung. 





Inhomogene Flüssiekeit 


| Homorene Flüssiekeit s \ i 
en £ (mit Diehteschichtung) 





| ohne Reibung | mit Reibung ohne Reibung mit Reibung 
l,ineares Profil | Rayleigh Sommerfeld, G. 1. Taylor 
v. Mises, Hopf. Goldstein 


P’randtl, Tietjens u.a. 








(ekrümmtes Profil Tollmien Schliehting 




















Mit den hier betrachteten Strömungen mit Diehteschichtung stehen in engem Zusammen- 
hang «die gekrümmten Strömungen einer homozenen Flüssiekeit. Z. B. liegt bei der ebenen 
Strömung zwischen zwei konzentrischen rotierenden Zylindern, von denen der innere ruht 
und der äußere umläuft, eine stabile Schiechtung dureh die Zentrifugalkräfte vor. Prandtl[6] 
hat 1927 durch eine Energiebetrachtung ganz ähnlich derjenigen über die Strömung mit 
Diehteschiehtung für eine dureh Zentrifugzalkräfte geschichtete Strömung eine Stabilitätsgrenze 
angegeben, die in guter Übereinstimmung mit Messungen von F. Wendt[20] steht. Den Ein- 
luß einer solchen stabilen Zentrifugalschiehtung auf die kritische Reynoldssche Zahl haben 
wir in einer früheren Arbeit untersucht, wobei das Resultat ebenfalls in befriedigender 
Übereinstimmung mit dem Experiment steht [21]. 

Kehren wir jetzt zu der Strömung mit Diehteschiehtung zurück, so läßt sich hier die 
Stabilitätsuntersuchung ganz allgemein folgendermaßen formulieren: 

Bei vorgegebener Laminarströmung U U (y) soll für jede überlagerte Störung (Wellen- 
länge 4) die Größe der Anfachung (Dämpfung) für jede Reynoldssche Zahl und für jede 
Schiehtung (Riehardson- Zahl) errechnet werden. Sei 4, die Amplitude der Störung zur 
Zeit #0, so hat man als Amplitude A zur Zeit #:A= A,-e’it, wobei also f; ein Maß für 
die Anfachung ist. Dieses Problem hängt von acht wesentlichen Größen ab, nämlich 

Pi 

), = Störungswellenlänge, 

d  charakteristische Länge der Laminarströmung, 
Maximalgeschwindigekeit der Laminarströmung, 
[#) Dichte, 
Io  Diehteunterschied in der Laminarströmung, 


loe Det-In)krement der Anfachung., 


» - kinematische Zähiekeit, 

g = Erdbeschleunigung. 
Aus diesen acht Größen lassen sich nach dem Prinzip der Dimensionpsbefreiung fünf unab- 
hängige Dimensionslose bilden, welche die charakteristischen Variablen des Stabilitätsproblems 
sind, nämlich 


I); 
— dimensionslose Anfachungsgröße, 
Um 
‚  dimensionslose Störungswellenlänge, 
pe ) ) r 
R=hevnoldssche Zahl, 
2 
y Io 
Vo y 
ir 9=0= hiıchardsonsche Zahl. 
| m 
3) 


lo , > U? lo 
"= relative Dichteänderung oder ER fo "roudesche Zahl). 
oO h) 4 O 
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Die allzemeinste Lösung des Stabilitätsproblems bedeutet also die Auffindung der Funktion 


ni N) ) RO | oO 


=. BB. a: 
) [#} 


Gl‘ 
m 
In einer früheren Untersuchung [22] haben wir diese sehr allgemeine Funktion @ von fünf 
Variablen dahin spezialisiert, daß wir für eine homogene Flüssigkeit (90, Ao/o=0) die 
Anfachung als Funktion der Revnoldsschen Zahl und der Störungswellenlänge berechnet 


haben, also: 

I); 0 1 

(Gr, ( - , )) . R\ 0 5 

Um’ 
Für das Folgende legen wir unseren Untersuchungen die in fast allen übrigen bisherigen 
Stabilitätsuntersuchungen verwendete Beschränkung auf, daß wir nur nach denjenigen Störungs- 
wellenlängen 4 =/, fragen, die weder gedämpft noch angefacht werden (indifferente Störungen). 
Wir nehmen also 5%; =0 an, und bestimmen statt (3) also die speziellere Funktion 


a3, R, oO, r) a 


die nur noch von vier Variablen abhängt. Gegenüber der homogenen Flüssirkeit ergeben sich 

dureh Einführung der Dichteschiehtung sofort zw eı neue Dimensionslose, nämlich © und 1olo. 

Dabei gibt © den Einfluß der infolge der Diehteschiehtung auftretenden pontentiellen Energie 

der Schwere, während 10/o den Einfluß der Diehteschiehtung auf die Trägheit darstellt. Das 
1 do FR 


Verhältnis dieser beiden Dimensionslosen, nämlich i F* ist das Quadrat einer 
Jo yo 


Froudeschen Zahl. Statt © und Ao/o können wir also auch © und F als unabhängige 
Variable wählen, so daß wir an Stelle von Gl. (4) schreiben können 

[, RO, Fl=0 ......2.2.. .(da) 

210 

Das gleichzeitige Hinzutreten von zwei neuen unabhängigen Variablen infolge der Schiehtung 
bedeutet ein außerordentliches Anschwellen der Stabilitätsreehnung, die schon im homogenen 
Fall recht umfangreich ist. Diese Funktion @, werden wir zwar nicht analytisch, wohl aber 
numerisch explizit errechnen. Um schließlich eine Beziehung zu gewinnen, die man bequem 
mit dem Experiment vergleichen kann, eliminieren wir am Schluß noch die Störungswellen- 
länge /,, indem wir speziell nach derjenigen Reynoldsschen Zahl R= R, (kritische Rey- 
noldssche Zahl) fragen, wo gerade noch eine einzige ungedämpfte Störung existiert, während 
alle übrigen gedämpft sind. Wir bestimmen also noch 

G,(R,, 9, F)=V0 
oder RR, RO, Be, Sin U 50 re a tn a 
und daraus eine kritische Riechardson-Zahl ©, unseres Profils. d. i. den größten Wert von 
), bei dem überhaupt noch eine ungedämpfte Störungswellenlänge existiert. Dieser hängt also 
noeh ab von der Froudeschen Zahl F: 9, 9, (F). Nach den Energiebetrachtungen von 
Prandtl und Richardson ist 9, = 1 unabhängig vom Geschwindigkeitsprofil: nach Taylor 
und Goldstein, die in ihren Untersuchungen die Trägheitswirkung der Schiechtung zegen- 
über der Schwerewirkung vernachlässigten, also F’=-0 annahmen, ist ©, '/, für gewisse 
lineare Geschwindigkeitsprofile und stetige Diehteverteilungen. 

Tollmien errechnete für die Plattenströmung bei homogener Flüssiekeit für die kriti- 
sche Reynoldssche Zahl den Wert AR, = (U „0*[v). = P7D, also R,;(9, F)= 9575 für 9-0 
und alle F. Es ist zu erwarten, daß sich für stabile Schiehtung (9 > 0) R,(9) >575 und für 
instabile Schiehtunge (9 <V) R;(9) <td ergibt. 

Wir werden unsere Stabilitätsuntersuchung mit Diehteschiehtung zunächst ganz allgemein 
für ein Profil mit von Null verschiedener Krümmung durchführen und, um Anschluß an Toll- 
mien undan experimentelle Ergebnisse zu gewinnen, als Beispiel dann die Laminarströmung 
an der ebenen Platte durchreehnen. Dabei machen wir noch die wichtige Voraussetzung, 
daß die Schiehtung nur in der Grenzschieht vorhanden ist, während außerhalb der Grenz- 
schicht, wo die Geschwindiekeit der Laminarströmung konstant ist, auch die Diehte als 
konstant angenommen werden soll. Für die Verteilung der Dichte nehmen wir später aus 
(‚ründen der mathematischen Einfachheit ein Exponentialgesetz an: 


O, () 0, me ’ 9, 
Ido 


» eine Konstante Ist. 
0 d „ 


4 


0 daß 
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Ss 2. Die allgemeine Störungsgleichung. Wir wenden uns jetzt den Einzelheiten der 
keehnung zu. Die ungestörte ebene Laminarströmung habe die Richtung der horizontalen 


„-Achse, und sei als eine Funktion der Höhe „ allein gegeben: U = T’(y) (Achse senkrecht 
nach oben). Die Diehte der ungestörten Strömung in der Höhe y sei bekannt: o,= 0, (NY). 


Wır beschränken uns auf den zweidimensionalen Fall, so daß die Beweruneseleiehuneen und 
dıe Kontinuitätseleiehunge lauten: 


[d " OH 7) Op Hal u 
S - v u 7773 
-\o0t Or Q = 0x 0x ( ff 
[d% Ar’ or) Op Pl OF r\ (7a) 
VO u > Ü r 3 a ji >| oA ' , . . .\td), 
r of (\ Y Q ‚' Q „ Q 4 ) ‚ ” / 
(\ [A ( / 
\ = v 
0) 0% 





Die von uns vorausgesetzte Inkompressibilität drückt sich dadureh aus, daß jedes Teilchen 
während der Bewerung seine Dichte beibehält, also 

Do 3) An) 0o 1 

: ‚tu > u 0 a Kur Se ee .(ıb). 

Dt Of Var u7, 
Die ebene Störungsbewegung wird als eine in der «-Riehtung fortschreitende Wellenbewegung 
anzenommen”). Wır setzen also 

7 Ik) + u v—v, 

i j ERSTE BE . (94) 
wobeı H u,(yye:\e N ’ / . , "(pe AN ft) | 


die Komponenten der Störungsbewegung sind. a=27[/4 Ist reell und bedeutet die räumliche 
Kreisfrequenz; 9 ist im allgemeinen komplex; = Pr t#Pß;. Pr Ist die zeitliche Kreis- 
[requienz der Störungsbewegung: pP; gibt die Anfachung oder Dämpfung, je nachdem ob 
positiv oder negatıv. 


Ferner setzen wir für den Druck p und die Dichte o: 


pP Pot) .- p' 9 Oo 0, (y) L 0’ 
| | ’ .(sb). 

It y ple'‘ Y ZiE 0 —=o,(ye'" \ ’f | 

Führt man die Gl. (Sa, b) in (7a,b) ein, so erhält man, wenn man im Sinne der Methode der 

kleinen Schwinzungen nur Glieder erster Ordnung beibehält, 


o . ) rs . 9 N zZ; 
0, tal Pu, tr] sap, Tul—-a’u, ru, ) 
FR pr | ’ tn N 
0, 18 1a PJvt, TIO,l pı ra a ( 
; (I), 
v, rtram, 0 
® u » | z ’ () 
ER PJLLTVYıB 





wo der Strich die Differentiation nach y bedeutet. Dies sind vier Gleiehungen für die vier 
Unbekannten ,, ”,, p,. o,, aus denen sich nach Elimination von u,, p,, 0, eine Gleichung 


für ”, ergibt, welche lautet: 


I 
’ 
0 


(«a U—BP(v, —a®v,)—alalU )Uw’ — IS gar, + U— pr, — ala U—-B)U'v,, 


(10). 





ivlal! 


u 


Dabei wird » 00, als konstant angenommen’). 
Wir führen jetzt noch eine Stromfunktion der Störungsbewegung 


va, oly)ar\“* 1) 
' ‚».0y Oy 
ein. so daß weren u „und " 
( „ 35 
", q (m): v, rayıy). 


') Es bedentet dies keine Beschränkung für die Allgemeinheit der Stabilitätsuntersuchung, denn H. B. Squire 


elrte (Pro Rov. Soe. A. Vol. 142, 1035), daß, wenn eine ebene Strömung geren dreidimensionale Störungen bei einer 


bestimmten Revynoldsschen Zahl instabil ist, sie dann für zweidimensionale Störungen schon bei einer niedrigeren 
tevnoldssch Zahl instabil ist. Die zweidimensionalen Störungen sind also „gefährlicher“ als die dreidimensionalen. 
’ (4) ur FE 

Danach ist also auch die Zähiekeit mit z veränderlich, und zwar so, daß eonst. Die Glieder mit 
099) 0) 


li-ssigen wir jedoch als klein von höherer Ordnung. 








",f.angew. 
und Mech. 


iten der 
zontalen 
'nkrecht 

9, (Y). 
ven und 





.(Ta). 


"eilehen 


iD). 


wegung 


. (Sa) 


ımliche 
Kreis- 
em ob 


a (Sh). 


de der 


(9, 


je vier 
chung 


(10), 


'Syqyuire 
ij einer 
"ieeren 
onalen. 


Nu 
nit - 
OO 


ea dc © a9 


Band 15, Heft 6 ' = oe ’ 
en Sehliehtinz. Turbulenz bei Wärmeschiehtung 319 





’ 


Setzen wir ferner noch e= pla, wo der Realteil von e die Phasengeschwindigkeit der Stö- 
unesbeweeung bedeutet, und führen wir noch dimensionslose Variablen ein, indem wir alle 
(;eschwindiekeiten auf die Maximalgeschwindigkeit 7/,, der Laminarströmung und alle Längen 
auf die Grenzsehiehtdieke 5 der Laminarströmung beziehen, so erhalten wir für die Strom- 
(funktion g der Störungsbewegung aus Gl. (10) die grundlegende Differentialgleichung 





do, 49 ) do ' 
‚ 2 of ’‚ 2 1 rY n ‚ es - | ne ® r 2,| 2 N fr; 
(I c) Wi Ray (I ce) | ‘ 0, dy Un? ‘ O0, dy | erg (I e)l /# 
(11). 
i I’ ) —„‚ ) “ 2; ' 1 0 do, 7 > ’ | 
Au RR‘ ei; u 2 ala Be Sch Ze O, dy u es | 
nn en u Y a 2 en C £ 
Dabei bedeutet der Strich die Differentiation nach 5 . 7 steht für ER: für ER. für «od 
‘ 77 In 


und R Z. B bedeutet die Reynoldssche Zalıl. 

Die Randbedingungen sind: Verschwinden der beiden Störungskomponenten für y 0 
und =», wenn wir eine Laminarströmung annehmen, die bei y„=0 durch eine feste Wand 
begrenzt ist, während sie nach oben unbegrenzt ıst: also gg’ 0 für y=V und „= x. 
Mit diesen Randbedingungen erweist sich die Stabilitätsuntersuchung als ein Eigen wertproblem 
() do, 4 2) () do, 


derart, daß sich zu jedem vorgegebenen Wertequadrupel «a, RR, ein im all- 


ur dy IR, . OD, dy 
gemeinen komplexer Wert von e ergibt, dessen Imaginärteil über Anfachung oder Dämpfung 
der betreffenden Störung in der betreffenden Strömung entscheidet, die durch eine Rev- 
noldssche, Riehardsonsche und Froudesche Zahl eharakterisiert ist. Wir beschränken 
uns Jedoch, wie schon in der Einleitung gesagt, auf solche Störungen, die an der Grenze 
zwischen Stabilität und Labilität liegen, für die also e rein reell ist. Für diesen Fall ist, 
falls in der Laminarströmung kein Vorzeichenwechsel von 1” vorkommt, e<T,, und a=_R. 
ls eibt in der Laminarströmung also eine Schicht, in welcher die Phasengeschwindirkeit der 
Störungsbewegung gleich der Geschwindigkeit der Hauptströmung ist. Diese Schieht wird 
als kritische Schicht (7 = y7;) bezeichnet und spielt im folgenden, ähnlich wie in den früheren 
Arbeiten. eine besondere Rolle. 


$ 3. Die reibungslose Störungsgleichung. Da an der Stabilitätsgrenze aA sehr groß ıst, 
können wir Aufschluß über die Lösungen der allgemeinen Störungsgleiechung (11) erhalten, 
indem wir die reibungslose Störungsidifferentialgleichung 


Ö do, 4 N) N) do, ‘| 
ur dy Be : 0 dy . 


0 


(U co” (p” ap) —(l e)U”g ce)’ (U co)U’ op) 0) (12) 
studieren. Wir denken die Geschwindigkeit der Laminarströmung in der Umgebung der 
kritischen Schicht, wo U ==e, in eine Reihe entwickelt: 


I c=U, y-y)tz Ur y—-y), 


die wir mit dem quadratischen Gliede abbrechen. Damit erhält man aus (12): 





LU: ee, 
+ a en, eat oe 
I 15 Au Un > 
Ky-+l, 7 Y,)7 3 1,’ (+ u‘ | I-+ 1’ (Y-HYn)\9 v (49). 
| | U; 4 
W-W)t5 17, WW Pd 
I Br I 
Dabei ist gesetzt: 
do, yo | q do, //dUN\? 
IN Fr yo = 
O, dy [ m ! Ta O, dy | dy I; 
(Ita). 
/ do, 
ody 
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Damit A und L in der Grenzschieht konstant sind, nehmen wir für die Dichte das einfache 
exponentielle Verteilungsgesetz 


O, Dee 5 Er ee rer 6 
an. Dann ıst 
yy | 
IN (au 
ee 





‘ 


Die reibungslose Differentialeleichung enthält infolge der Diehteschichtung die beiden Dimen- 
sionslosen AK und L, wo RK die Schwerewirkung und L die Trägheitswirkung der Schiehtung 
auf die Störungsbewegung darstellt. G. 1. Taylor und Goldstein vernachlässigen in ihren 
Untersuchungen fast durehwer die Trägrheitswirkung gegenüber der Schwerewirkung, was 
nach der Differentialgleichung (13) nur erlaubt zu sein scheint, falls 7, N, oder nach (14b), 
dU\® ö BE, 
falls | gu 3 s l, wenn also die Froudesche Zahl der Strömung klein gegen 1 ıst. Aus 
EY/R I 09 
unseren weiteren Rechnungen wird jedoch hervorgehen, daß bei gleichem Zahlenwert von K 
und Z der Einfluß der Schwerewirkung (A) auf die Stabilität weit größer ist als derjenige der 
Trärheitswirkung (7). Bei den Messungen, mit denen wir unsere Reehnung vergleichen wollen, 
z. '1,d.h. L>K. Da sieh nieht ohne weiteres abschätzen läßt, wieviel in 
diesem Falle die Trärheitswirkunge ausmacht, führen wir im folgenden die vollständige Rechnung 


Ist soeAar 


mit Schwere- und Trärheitswirkune dureh. 
Wir kommen jetzt zur Integration der Differentialgleichung (13). Die Lösung dieser 
Differentialeleiehune denken wir nach Potenzen von K und ZL entwickelt, 


g=yP+KyW+HLPEWH.. 2. 2.2.22. (16), 


wobei wir die Entwieklune mit den linearen Gliedern in AK und ZL abbrechen. Wir wollen 
also, auszehend von der homozrenen Flüssiekeit, die Stabilität einer schwach geschichteten 
Strömung untersuchen (K = L =0 ist die homogene Strömung). Dabei ist 4°” die schon aus 
früheren Untersuehungen bekannte Lösung der Störungsdifferentialgleichung der homogenen 
Flüssigkeit und g'N) und 9» sind Zusatzlösungen infolge der Schiehtung. Diese Entwicklung 
und insbesondere das Abbrechen mit den linearen Gliedern ist nur zulässig für Werte von 
N und Z, die klein gegen 1 sind (vel. Fußnote S. 328). 

Da wir die Laminarströmune, deren Stabilität wir untersuchen wollen, durch eine kon- 
stante, eine lineare und eine quadratische Funktion annähern werden, integrieren wir im 
folgenden Gl. (15) für konstante, lineare und quadratische Geschwindigkeitsverteilung. Dabei 
eehen die Reehnuneen für die Zusatzlösungen g(R) und ge) Immer ganz parallel. 

Für denjenigen Bereich, wo die Geschwindigkeit der Laminarströmung konstant ist, 
wollen wir, wie schon in der Einleitung gesagt, auch die Dichte der Flüssigkeit als konstant 
annehmen. Man hat dann dort als Störungsdifferentialgleichung 

g" ag 0 


mit der für große y abklingenden Lösung 


7 me @ ay ° . . . . . . . . . . . . . (IT). 


Für lineare Gesehwindiekeitsverteilune erhält man mit der neuen Variablen 


Ya ö 


Un (18) 
um mit any; a, | 
dl? ö ) dq\ 
aus (13 ° _ q)+ Kyg-+ Ly; ıy,e Ru I | | 
») Y, Frr k, 1) yrHuyı 11 7 UF dy) (19) 
Für «lie homogene Flüssiekeit hat man aus 
d’ a‘ 
—. a 0" UV 
Ay? 
die Lösungen 
(0) Zina,y, (0) fd ») 
p, . Ds bo} a,y, -» ’ er UV). 


A 
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Für die Zusatzlösung y‘# erhält man 
oleicehung 


aus (16) und (19) die lineare inhomogene Differential- 


d: (Kı 2 
u / ei PaLL) 1; 
Y, | dy, a, y(K)) ] EEE Er 
und für die Zusatzlösung g: 
dgl | da 
2) 2 or... eh) 
y,| dy? 4,49 I 7 UF dy,) ) 


wobei der homogene Teil von beiden Gleichungen derselbe ıst wie die Differentialgleichung 


für 9”. Man braucht also nur noch 


ein partikuläres Integral von Gl. (2la) und (21b) zu 


bestimmen, um die allgemeinen Integrale y%) und 99 zu erhalten. Diese partikulären 
Inteerale lassen sich aber oline weiteres durch Quadraturen bestimmen. Man erhält eın 
Fundamentalsvstem für die Lösungen y'K) und g9 von Gl. (2la, b) in der Form: 


Yı 





7 (0) 7 (v0) 7 ‚0 y ER * (0) 7 Fa 
7 (K) — —. y Wh \ x a d „, - . \ 2 dy, 
Y:) 4 Ye Mm 
Jo Yo 
(22). 
Un Yi 
» 0) (u) ” (0) / (u) 
/v dq v (0) ltr “rPy .) 
ın.(l) me (” (v Jdı r op, ( | dıy, iv 1,2) 
/v Y/v Yı \ 2 | Y, d y, Yı Tr Ta \ fı Y, d Y, Mi 
Yo Yo 


Die Wahl der unteren Integrationserenze ist gleichgültig, wir wählen als solche den 
Punkt 7 9,. wo wir das lineare Geschwindigkeitsprofil an das parabolische anschließen. 
Für parabolische Geschwindigkeitsverteilung, die wir in der Form 


N 


schreiben (a = gegebene Konstante), er 


UE 


und mit (ko 


u sea rar 
hält man mit der neuen Variablen 

Y%Yk 

(Ad LER 


(a — Yr)a 


aus Gl. (13) als Störungsdifferentialgleichung: 


6 2 (1. r 
Yo" | l 2 + on: Ay | + UR | | 


n Y+Ky 





(24). 
} UR 4 UER ’ day | 
+ IL (a Y7;) 1 > | Y,.)4 uUR 1 2 dy,) 0 
Die Lösung y“ dieser Differentialgleichung für die homogene Flüssigkeit (X © L=0) ist schon 


von Tollmien und in unseren früheren Arbeiten berechnet worden.  Eın Fundamental- 


system ist gegeben durch 





(0) 
| d, Ur fr d, Ya L- d; Ya Be 
a —Yr 
(0) 
(0) 2 1 ) > 
e,te.u.te.y’-+. or Yy,, 
7 > 0 ı Js + » Ms A Y7; N: 
| Be A 
wobeı d, | M d, ee 2 s d, Fr. 6 d, 15 ) 
| - (29). 
d: 7 BULAU 3)d„_, 2a’ dns — a,’ d.-;| n„=2,8 
a,” | u. 
eo | , e, —=0; € si > l e; S gi 15 
| ) vr Bi. 5) S)d 192 Nd | 
a Fe ie On 3) en, +2? e,n_. ar en 3 (an —3)d„_,+2(2n ’ 
n=EB,.. 


Für die Zusatzlösungen y'R) und g'L 
zleichungen 


erhält man aus Gl. (16) und (24) die Differential- 


2 Ur „ j 7 (R) 2 K | Ns K {07} 
y|ı > | dy, a, p\ u +, (1 > |g' ” J 
und Ns (dl? q (L) ö ’ | | (L) | 2 | Y,, dq m. 
(1 2) ET au TYP Tg day) 
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Der homogene Teil dieser inhomogenen Differentialgleichungen für y'$) und y) stimmt wieder 
überein mit der Differentialgleiehung für 4°”, so daß man ein Fundamentalsystem für die 
alleemeinen Lösungen y'® und 97) wieder durch Quadraturen erhält, nämlich: 





I2 VE 
(ı) D (0) (u) cu) . (0) (0) 
f (KA) ! ws L 1 L < 4 4 \ 4 - hi / v ER .) 
2 d; aM . Ya\" Ye MY ’ Ya\" es (v=19 
„ 1 ) | Ns 1 ) 
4 Yo 
/ (0) (0) (0) (1 7) vo a 
Ey‘ /v ‘/ı Ya } „, ‘v dm, (26) 
| A: | = 
Ze 
[7 (u) 
(0) (0) (1 N3)4 ‚ +. ‘ 
Ya \ Ya | | , ‘vr | dy, (y 1.2) 
1 —5 
(ra) | 


jo 


Damit sind die reibungslosen Lösungen der geschichteten Flüssigkeit für lineare und para- 
bolische Geschwindigkeitsverteilung gewonnen. Von diesen Integralen ist nur 9,” im sinzu- 
lären Punkt 9° 77, der Differentialgleichung regulär, während 9," 
y,+ dort Singularitäten aufweisen, indem sie oder ihre Ableitungen Terme mit log (y- 91.) 
bzw.(y7— y7,) ' enthalten. Diese Singularıtäten rühren offenbar nur von der Vernachlässigung 
der Reibung her und fallen fort, sobald man diese in richtiger Weise berücksichtigt. Die 
Integration längs der reellen y-Achse ist in der Nähe des singulären Punktes y=y), zu er- 
setzen dureh die Integration im Komplexen längs eines Halbkreises um den Punkt y= y7. 
Dabei ist es zunächst unentschieden, in welchem Sinne man den Punkt y= y) zu umkreisen 
hat, d. h. welehen Zweig des log man beim Übergang von positivem zu negativem y— y7; zu 
wählen hat. Tollmien hat für 9," diese Übergangssubstitution angegeben, die man erhält, 
wenn man in einer kleinen Umgebung der kritischen Stelle y== 7 die vollständige Störungs- 
differentialgleiehung (11) diskutiert. 


: y,'W) und 7 „(R), pl, 


Ss 4. Verhalten der Lösungen in der Umgebung der kritischen Schicht. Hierzu nimmt 
man eine kleine Umgebung um die Stelle „= y7, (Übergangsgebiet) an, in der man mit 
cenügender Genauigkeit 7 e durch U, (y 97) und U’ und U” durch U,’ bzw. U,” er- 
setzen kann. Ferner setzt man 





s— a=—laRD, I ’ıgsen: : “3 5 8 8 


Dabei kann weren des kleinen Wertes von e selbst im UÜbergangsgebiet (d. h. kleine Werte 
von 9 29,) 7 große Werte annehmen. Für 9 (7) ergibt sich aus Gl. (11) und (27) die Differential- 
oeleiechune 

I. I; 
2iad)—nyle— te an ira) PA VToTr re p—-ngpoy4=0 (28), 
\ Un I: 

I ö do, 
R 0, dy 
sind. Tollmien[13] hat an Hand der entsprechenden Differentialgleichung für die homogene 
Flüssigkeit (N 10) gezeigt, daß die Übergangssubstitution für g,” folgendermaßen lautet: 





wobeı ın Gl. (11) die Terme mit als klein von zweiter Ordnung vernachlässigt worden | 








Ur 
„ In 0 . y ne e, | e, „ -P, “ n— n : y a log ( Y7;) 
. * - ° - I; 
2 | UT .\ 

y—yr<0:gQ, te y+a +... + 77, m" dlogiy-yrl-i7, 

« . “ « -n IR ü a 
ls erweist sich also von den unendlich vielen Werten von log (y - ) +2 ki (k = positive 
oder negative ganze Zahl) derjenige mit A = I als physikalisch reell. Durch eine Rechnung, | 
die ganz analog der für y,"” von Tollmien durchgeführten verläuft, erhält man für unsere 
/usatzlösungen infolge der Schiehtung @,'®%, g,'P) und g,) dasselbe Resultat, daß nämlich 
die Terme mit log (y - 7,) für negatives y - 9y, zu ersetzen sind dureh log 'y— y.  im. | 

Ferner ergibt eine Näherungsreehnung, die in der Differentialgleiehung (25) nur das 

erößte Reibungsglied berücksichtigt dy”’"-+ng" 0), für die homogene Flüssigkeit noch 
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zwei weitere Lösungen 9,” und g,”, die man zur Darstellung des allgemeinen Integrals der 
Störunesdifferentialeleiehung vierter Ordnung (11) nötig hat. Man erhält 


1} r 
(u ) 71. C2) fe .g ‘ 
13.4 \d\ 7 » Hı,, | „een )ld ı, rer . (30) 
FT. HF 


wo JIV»@ die Hankelsche Funktion erster bzw. zweiter Art bedeutet. In ähnlicher Weise, 
wie wir in $3 aus der reibungslosen Differentialgleiehung Zusatzlösungen infolge der Schiehtung 
vewonnen haben, hätte man auch für g, und y, solehe Zusatzlösungen zu berechnen, nämlich 


m) -Ky®) r“ iq (L) 


ara 34 3,4 1 3,4 
® a ' i 0 = RK) 
Wir werden jedoch im nächsten Paragraphen zeigen, daß sich die Berechnung von 4, N und 
(L) 4 . er . 
7. , erübrigt, da nämlich 
(K) (K) (L) (L) 
734 _712 734 1712 £ 
y m und —ıa er: Te 1° 
34 Tı>» 34 Iı» 


Der Einfluß der Diehteschiehtung auf die reibungslosen Lösungen g,,, Ist also sehr viel 
oerößer als auf die Reibungslösungen 9,,,. So daß man für unsere Näherungsrechnung den 
letzteren vernachlässieen kann. Dadurch wird die Rechnung erheblich vereinfacht. Dieses 
Resultat steht in Analoeie zu den Verhältnissen bei einer Strömung mit Schiehtung dureh 
Zentrifugalkräfte. In unserer früheren Arbeit in den Göttinger Nachriehten [20] haben wir 
ebenfalls gezeigt, daß bei stabiler Schiehtung durch die Zentrifugalkräfte der Einfluß dieser 
Schiehtung auf die reibungslosen Lösungen q,,, sehr viel größer ist als auf die Reibungs- 
lösungen 9,,4- 


Wir geben jetzt diesen Beweis: 


$ 5. Der verschiedene Einfluß der Dichteschichtung auf die beiden Lösungspaare. Dazu 
[ragen wir nach einer Darstellung der vier Integrale 4, 9. %, 9%, für die Wandnähe (sehr 
kleine 9) unter der speziellen Annahme, daß die Phasengeschwindigkeit e der Störungs- 
beweeune sehr viel größer ist als die Geschwindigkeit in Wandnähe. Es soll also die kritische 
Stelle außerhalb der Schicht liegen, wo wir die Lösungen betrachten. Wir können dann für 
die Wandnähe U —-e ersetzen durch ce und U” durch U”,,. 
I do, 
o,.dy 
meine Störungseleichung nach Gl. (11): 


Weeen nach Gl. (15) lautet dann für eine wandnahe Schicht die allge- 


’ { | 
J a? ec) 7 (). 


e 


| a 
r . MN ‚ | > r / 
7 "Hlia Re Zar)gq s -a+ta R|t Mi 3 
| 
Die Lösungen dieser Differentialeleiehung mit konstanten Koeffizienten sind von der Form: 
ur ‚=, 234. 
Dabei sind die k, die vier Wurzeln der Gleichung vierten Grades 


ran : nl ns? 
"+(iiaRe-2a)®+|+iah | Un 7 
| | 


“ a? e) v), 


Die beiden Lösungspaare sind gegeben dureh 


P=- - iaRc— 2a?) + | (aRc—2a” — —iaR|\Un + Fr ae), 
woraus folgt 
k' ‚a? ln 
oder, wenn 1." =V, er 
k" a a ! K| ZT } da A a 
SE 
oder Ben N = I 








by 
n 
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Hbenso findet man, wenn man a? eeren a Re vernachlässiet. 


; ’ ‚ Du u 
l- iaRc+K| 
3.4 p 
| K ia /U.\?’ 
I; 2 + ] (dA Rr | ) | | 
oder 23 Rı 1Cc 
Aus 
K /U,r’ 
(0) , RK all a) 
Yı Yı HA Yı e s 
findet man 
y,'®) alte 
"> Pradaaia ) (dt | 
/ 2 1 
und ebenso aus 
K/ı Urn 
(0) ‚ RK Vie ref!  (5;)\ = J, 
7/3 3 ya Ih VE =e A VORR ande 
() ‚(KR [7 4 | 1% | va 
7 a ) Ad f I? r . 
Damit ergibt sich 
(0 
/ a Y Zi >.) 
/ " (0) > . . . . . . ä ’ Mi : R , 2 (32). 
Yı ie, VER li C 


Für die Strömung längs der ebenen Platte ist nach unseren früheren Untersuchungen [22 
mit d als Bezugslänge an der Stabilitätsgrenze näherungsweise: 


7 0,S1: R 100: e—=V42. 


Damit wird 


an.(K wm (% N | 
ph) m,” 3550 29 


Dasselbe gilt für y,'®) und y,'®%), und eine ganz Ähnliche Ableitung läßt sich auch für y,) 
und g,7) geben, womit also Gl. (31) bewiesen Ist. 

Damit ist also gezeigt, daß der Einfluß der Diehteschiehtung auf die mit y langsam 
die sich wie e1@ verhalten, sehr viel größer ist als auf 
"VaRey 


veränderlichen Lösungen 9, und g,, 


die rasch veränderlicehen Lösuneen, die sich wie e verhalten. 


$ 6. Formulierung des Eigenwertproblems. Wir kommen jetzt zur Aufstellung der Rand- 
bedingungen und damit zur Formulierung des Eigenwertproblems. Wir legen ein Geschwindig- 
keitsprofil zugrunde, das vom Wert Null an der Wand ansteigt bis zum maximalen Wert 1, 
und diesen dann konstant beibehält. Für das Gebiet konstanter Geschwindigkeit, wo wir 
auch konstante Diehte annehmen, lautet die Lösung der Störungsdifferentialgleichung nach 


Gl. (IND): y=e *%. Bezeichnet y„=a die Anschlußstelle an das Gebiet mit konstantem 7, 
so hat man in diesem Punkt die Randbedinzeune 
Da +ao DE Pe ‚(83 


ra l id 


Das alleemeine Interral lautet 


‘ ( 9 1 | v .+0,1 ( 7 N 
Die Randbedingungen (y 9’ =0 für y=VO und y=-) lassen sich auf Grund der besonderen 
Kigenschaften der vier Integrale 9,, .. y, erheblich vereinfachen. Zunächst braucht das 


Integral y,, das für große y sehr stark über alle Grenzen wächst (wie e VaRcy) in dem all- 


eemeinen Integral nicht berücksichtigt zu werden, also C,=0. Ferner braucht das Integral 
y,, das sich wie e \«ReY verhält, nur an der Wand (y==0) berücksichtigt zu werden, während 
es im Anschlußpunkt an das Gebiet mit konstanter Geschwindigkeit (ya) mit sehr guter 
Annäherung als Null angenommen werden kann. 
Für „=a lautet demnach die Randbedingung gemäß (35) einfach: 
CYPiattlzY aa allg 1a + (0,9 a) 0 
oder, wenn man 
D,, IN er ET . (33a) 


setzt, 
BELEGE. dr en ar EEE: 
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Für die Randbedingung an der Wand (y=0) ergibt sich 


C Yıw + C,q sw tr Ya w v 
(3tb, e). 


’ 


| 
C,o9,u+ (C,9 . u rU,g sw 0 | 


Damit hat man als Randbedingungen die drei homogenen Gl. (34a,b,e). Falls eine von 
Null verschiedene Lösung y vorhanden ist, muß also die Determinante 


D, 2 PD, 5 () 


/ıw sw Isaw - B , ] - : ‘ 2 (39) 


’ 2 
I ıaw Jaw I saw 


verschwinden, womit sieh in bekannter Weise die Stabilitätsuntersuchung als ein Eigenwert- 
problein der Störungsdifferentialgleichung erweist. Die komplexe Gl. (35) hängt außer von 
den Konstanten des zuerunde gelegten Profils ab von den Parametern «a,e, R, RK, IL, die sämt- 
lieh, insbesondere e, rein reell sind, da wir nur nach den Durehgangspunkten von Stabilität 
zu Labilität fragen. Die Funktionen 9, und P,. Sind im allgemeinen komplex [(vel. Gl. (29) 
und (30). Denkt man sich aus den beiden reellen Gleichungen, mit denen die komplexe 
(l. (35) Äquivalent ist, den Parameter e eliminiert, so hat man noch eine Gleichung, in der 
a. R. K. Iı enthalten sind. Diese Gleichung gibt bei festzehaltenem A und 1 die zu dieser 
Schiehtung (K, L) gehörige Indifferenzkurve in der a R-Ebene, welche die gedämpften Störungen 
trennt von den ungedämpften. Die Lösung des Stabilitätsproblems kommt jetzt also hinaus 
auf die Diskussion der Gl. (35). 


Durch Ausreehnung der Determinante (35) erhält man 


( 
au Jow D, a Yıw PP. a 
(6). 


[4 


’ ’ 
‘D «DD 
saw I aw / Aa Jam / A 


Wir entwickeln jetzt diese Gleiehung nach Potenzen der Schiehtungsgrößen AK und 7, ent- 


sprechend dem Umstand, daß wir, von der homogenen Flüssigkeit ausgehend, die Stabilitäts- 
untersuchung für schwache Schichtung durchführen wollen. Dabei brauchen wir nach den 
im vorigen Paragraphen angestellten Überlegungen nur die rechte Seite dieser Gleiehune zu 
entwiekeln,. welehe ausschließlich die reibungslosen Lösungen enthält, also nieht von der 
keynoldsschen Zahl R abhängt. Andererseits hängt die Lösung g, von der zugrunde 
eelerten Laminarströmung nicht ab. Wir können sie wie in der Tollmiensehen und in 
unseren früheren Arbeiten ohne weiteres von Tietjens übernehmen. Die linke Seite von 
(+1. (36) hängt nur ab von 








rs un = 
Nur Y,\a li N ) er i . } : i ; i u ö ; (4 ). 
Ann | 47 3 W n. riv * ® . rg ryv 
Für Dina) ; erhält man nach Tietjens die untenstehende Tabelle (Tabelle 1). 
EG zu 
Tabelle 1: D(y,.) als Funktion von 7,.. 
In IB \ Y [Zi u /) \ 7 u) 
v 0,702 0.425 ü 3.0 1.100 VD1D 4 
0.5 0.785 V.4ll 3.5 1.180 — 1.130 i 
1.0 0,020 0.389 ; 1.0 0,4160 -—- 1.250 
1.D 1.043 0,297 Gi 1.5 0.0405 0.8080 ; 
>20 1.206 0.147 Gi 5.0 0.0057 -+- 0.3645 7 
95 1.357 + 0.108 i 5,5 0.1913 + 0.2393 i 


Setzt man ferner noch 


| vi 210 »D Ad (7 m PD, 7 


; KlacÄ2) : : : 3. » 188), 
JE T aw Dia Yımw PD, a 
so erhält man aus Gl. (36) 
Dir) 
— E(a,c,K,L) . . u 2 2 + ID. 
Iw 


Diese komplexe Gleichung bildet den Ausgangspunkt für die ganzen weiteren Rechnungen, 
deren Ziel die explizite Darstellung der durch diese Gleichung gegebenen Beziehung zwischen 
den fünf Größen «a, e, R, K und L ist. 
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$S 7. Lösung des Eigenwertproblems. Für die numerische Behandlung dieser Gleichung 
ist es wesentlich, daß von den fünf Parametern auf der linken Seite die beiden Schichtungs- 
erößen K, L und auf der rechten Seite die Revnoldssche Zahl AR nicht auftreten. 

Die Auflösung dieser Gleichung, also die Berechnung der zusammengehörigen Eigen- 
werte a, ce, R, K, L geschieht zunächst analvtisch, soweit es sich um die Entwicklung nach 
Potenzen von AK und Z handelt. Die weitere Behandlung der damit gewonnenen Gleichung 
muß in ähnlicher Weise, wie bei den früheren Stabilitätsuntersuchungen [21, 22], auf 
numerisch-graphischem Wege geschehen. Die ganze Rechnung Ist sowohl im analytischen wie 
im numerisch-graphischen Teil außerordentlich umfangreich, so daß es nieht im entferntesten 
mörliech ist. sie hier einizermaßen lückenlos wiederzureben. Wir müssen uns deshalb auf 
eine kurze Andeutunz des Gedankeneangees beschränken: 

Die Entwieklung von Gl. (39) nach Potenzen von AK und Z läßt sich formal schreiben: 


Din ,,.) FRI RK OK N OH 39.) 
9a c)+-K +1, = (39a), 
ET: oA Of, 
wobeı E° den Wert von Gl. (35) bedeutet, wenn man dort für 9,, g, die Lösungen für die 
homorene Flüssigkeit einsetzt, also 
, (0) op ’ (1) op 
aan | rIw la ur 2a 
BR” (40), 
Y;, ((1) <b\ ’ (il) (u) 
2 Dur la / Lır Aa 
OE OH 
.. . e .. us .. ® ryv . +.) ® v ° ” 
während man \ N un«dl \ erhält, wenn man für sämtliche Terme in (35) ihre Entwicklung 
( \ ( 4 
nach A und Z einführt, also 
(0) . K), (1) 
JR 7 vu I\ Dan I; 4 vw a | 
(v tee. SE FE 
$D p np L np") N | 
ya 7 va ' vAa hie 
(AK) (1.) , N) (2) j en 
Yy , 9, und damit P,, und P,, sind nach $ 3 bekannt. 


Die komplexe Gl. (39) ist Äquivalent mit zwei reellen Gleichungen. Für die linke Seite 
von Gl. (39) ist die Zerlegung in Real- und Imaginärteil nach Tabelle I ohne weiteres ge- 
eeben. Die rechte Seite enthält einen imaginären Bestandteil infolge der Übergangssubstitution 
des Integrals 9, (Gl. (29)). Führt man die angedeutete Entwicklung nach Potenzen der 
Schichtungsgrößen A und L explizite durch, so erhält man nach langer Rechnung Ausdrücke 


\ FE u ? 
= a ' u. GM u, ; Se 
für den Real- und Imaeinärteil von E", Wr und \j als Funktion nur der beiden Variablen 
( \ ( 4 
. f (0) Q FÜ Ö F Han .. ‚ . . 
ande. Die Zahlenwerte von E, | Kol für das noch näher anzugebende Geschwindig- 
() A (! 14 


keitsprofil, mit dem die Rechnung durchgeführt wurde, sind in Tabelle IT angegeben. Die 
weitere Auflösung von Gl. 89) geschieht nun auf graphischem Wege. Man trägt bei fest- 
eehaltenem e, AK, L die linke Seite von (39) mit 7), als Parameter und die rechte Seite mit «a 
als Parameter in einem Polardiaegeramm auf (Abb. 2). Im allgemeinen erhält man zweı 
2 l i 
Schnittpunkte der - Kurve mit der E (a, e, K, L)-Kurve, zu denen die Parameter «a,., «, und 
Hau 
rs» rs gefunden werden. Die zugehörigen Revnoldsschen Zahlen R erhält man dann 
aus der Gleichung 


ww yrlak U, )®. 


) 


Der große Vorteil der Entwieklung nach Potenzen der Schichtungsgrößen AK und L 
macht sieh für die numerische Rechnung dadureh geltend, daß man die Abhängigkeit der ge- 
suchten Durehgangspunkte von Stabilität zu Labilität von dem Schichtungszustand in sehr 
. r Er en m 0E 0E . 
einfacher Weise beherrscht. Hat man bei festem e die Größen E, as x einmal als 

oh’ oh 
"unktion von « berechnet, so hat man die Kurven F=E(a,e, K, L) sofort für die homogene 


Strömung und für alle Scehiehtungszustände. Aus den in Tabelle II angegebenen Zahlenwerten 
N ı N ' 
von. ; une N geht hervor, daß der Einfluß der Schwerewirkung der Schiehtung sehr viel 
( \ ( d 
OH 


I. 

n" 0 
nungen haben wır auch kleine negative A- und L-Werte mit in Betracht gezogen, also 
schwach ıinstabile Schiehtungen, z. B. die Strömung eines gleichmäßig temperierten Luft- 
stromes oberhalb einer geheizten Platte, die ein wenig wärmer ist als die Luft. 


N 


. .+ s ryr a . . ( ‘ .)» 4 
größer ist als der Eintluß der Trägheitswirkung |, vel. S. 332). Bei unseren Rech- 
( 
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S 8. Anwendung auf die Plattenströmung. Wir haben für unsere Stabilitätsuntersuchung 
mit Diehteschichtung als Beispiel für die numerische Durchrechnung die Laminarströmung 
längs einer ebenen Platte gewählt. Die entsprechende Untersuchung für den Fall der homo- 
veenen Flüssigkeit hat früher Tollmien[l3] durchgeführt. Die Berechnung der laminaren 
Geschwindigkeitsverteilung längs der ebenen Platte führte zuerst Blasius[23] nach Ansätzen 
von L. Prandtl[24] aus. Für unsere Stabilitätsuntersuchung können wir darauf verzichten, 
die Abänderung der Laminarströmung durch die Dichteschichtung in Betracht zu ziehen. 
Durch eine einfache Abschätzung’) läßt sich nämlich zeigen, daß diese klein ist, wenn eine 
schwache Diehteschiehtung nur in der Grenzschicht vorhanden ist, wie wir es vorausgesetzt 
haben. 

Für unsere Zwecke brauchen wir eine angenäherte Darstellung des Blasiusschen 
Profils dureh eine Gerade und eine Parabel (Abb. 1). Wır lassen das Profil an der Wand 
(4° 0) mit einem geraden Stück beginnen, das wir tangentiell an ein parabolisches Stück 
anschließen, und dieses wiederum tangentiell an die konstante Geschwindigkeit U = U 
Wir erhalten somit die schon früher für U angegebene Näherungsdarstellung 


ım* 


























„7 ae I „ 
(0 3 -0.1758: r 1.655 y 
Un 
Zn SR. > ae 
<< 7 = | (42). 
‘ m ) 
7 56 
3 1.015: j | 
( IH 
| ra 
| 
y/d 
! | eu 
| 
| 
70 7 
Papnf j Larf 
VESEHW 77 DMEN, 
/empera 
\ 
Fo 
\ 
Ph 
ei A 
ar 2 Ce 
u 50 TC mW nn ; 
Ös x 3 7 7 74 ıE) 7 £ np 
0 05 U/ im 20 R/E) 
7 Dlasıus Profil j My Krs0 
f : Nähe y Ysfunklion Ari i (7 . EEE i i i i i 
) ‚ Die tersuchte Laminar D ’ : br 
Abb. 1. Die BIENEN, Lamin Abb. 2. Polardiagramım Kurven und E («,e, K,L)-Kurven für 
strömunge mit Diehteschichtung. Nun 


c 0,10 und A N, : . - a E 0. 
100 760 107 30 


Man hat also a 1,015 dureh Vergleich mit Gl. (23). Für die Übergangssubstitution 
Gl. (29)] darf man Jedoch nieht im linearen Bereich U” = const und LU” = 0 und im para- 


5), Die dureh die Sehiehtung hervorgerufene Kraft in der x-Riehtung beträgt pro Volumeneinheit: 


Id do do 
(vd) a: )Q d) 
= au üx a Ax 
' i dU 
und die Reibungskraft pro Flächeneinheit 
dy 
/ dl «U 
oder pro Volumenerfnheit e A ö 
J d 17] 02 


Ks ist demnach erlaubt, den Einfluß der Diehtesehiehting auf das Laminarprofil zu vernachlässigen, wenn 


do "Un 


u) 
(] 
ee 7 d 
Nun ist #. di Um? e UppnX dd . j : ) ’ do Be 3 
Eau ) dx UnXx Um‘ dx Um 
a } l I [7 } 
Ks muß also sein 7% * oder —.<]1, 
Um ur Uml9)? 


und das ist gerade die Voraussetzung, die wir schon oben machten. 
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Die Durcehgangspunkte a d* und 


für verschiedene Werte 


der Schiehtungesgrößen A und L. 





















































. 1 | 
N: m I, () N 00 
C ( a" j (aö®), (' m |, (a 6*), ( (' n_ | (a 6*), (' a" \ (a 6*), 
0.20 1.15 - 10° 0.038 0,10 1.05 - 10° 0.053 4.05 - 10° 0.080 
0.30 1.84. 10° 0,OTS 8.75: 10° 0.187 0,20 7.16. 103 0.005 2.57-10% 0,164 
0.35 1.01. 10° 0.101 2.99. 10° 0.232 0,30 1.57 - 10° 0,148 3.88: 103 0.232 
0.375. 7.25-10? 0.121 1.8 - 10° 0.251 0.35 462.10? 0,185 1.854 - 10° 01.256 
0.40 5,86 - 10? 0.139 1.30. 10° 0.266 0.375. 7.85- 10? 0.207 1.26 - 10° 0,267 
0.45 1.04. 10° 0,193 6.08 103 0.282 0,40 1.67.10? 0.251 5.60.10? 0.265 
| 
1 Ä su 
IN 
60 0.10 1,05 - 10° VOGA 05.109 V.OOS 
1.2) 4,27 -10° 0,050 0.20 7.30: 103 0,108 > 11-10 (1.17? 
0.30 1.68 10° 0.003 6.75 - 10° 0.202 0.30 1.66-10° 0.165 2.48: 103 0.234 
1.32 420).10° 0,121 >61 :-10° 0.240 0.35 1.06.1093 0.204 1.69: 103 V.>D8 
0.375 |, 7,15-10? 0,140 1.70. 10° 0.257 0.375 . 897-10? 0.230 1.11- 10° 0.208 
0.40 5.70 -10° 0,159 1.20.103 0.271 | 
0.49 t.>1- 10? 0.233 5.23.10? 0.264 N (0 
0,10 1.19. 10° 0.079 2.24. 10° 0,104 
K I 0.20 808-109 0.125 1.71.10 0.170 
100 0,30 1.00. 10° 0.100 3.14. 10° V.>38 
0,20, 8,22. 10° 0,060 0.35 ‚21-10° 0,236 1.47: 103 0.255 
0.30 1.60» 103 0.106 2.71. 10° 0,211 | 
035 8.96- 10°? 0.136 2,36 10° 0,245 K=. 
0,319 ER AN 0,500 1,89.10° 0,002 0,10 1.4? 10° 0,007 1.68: 10° 0.105 
0,20 9,40. 10° 0,100 1,12.10° 0,208 0,20 4304-10? 0.145 1.44- 10 0.183 
0.42 2.18 10° 0.198 8,28 - 10° 0,20 0.30 > 10-103 0.211 > 73-103 0,236 
0.325 1.66 - 10° 0.232 1.00. 10° 0,244 
K==% 
(1.20 7.20: 10° 0.077 3.76. 10% 0.149 RK = 
0.2) 01 - 10° 0.101 1.20: 10° O,ISS = — 
0.30 1.53: 10? 0,129 1.64 109 0.223 0,20 1,04.10° 0,106 1,22. 10° 0,204 
0.323 1.15-103 0.113 320.103 0.238 0.225 6.51-10° 0,181 7,44. 10° 0.195 
0.3 S.03.10? 0.150 27-10? 0.251 0,25 6,05. 10° 0,186 7,09. 10° 0,199 
0.375 7.36.10? 0,181 1.42: 10° 0.261 Re 
0.40. 633-10? 0.205 1.02-103 1.274 24 
0.42 3.76 - 10? 0.239 7.13 - 10? 01.275 0.225 7.18- 10° 0,190 
K — - I; r u, 
C (Om e (a 6*), ( = a (a Ö®), C * en ei (a6), (‘ m (a 0), 
v f ) 2 vw : v 
0,30 1.94 - 10° 072 0,20 S39.10° 0.065 
0.35 994-102 0.102 3.16- 10° 0,218 0,30 1.07 - 10° 0,113 2.21- 10° 0,207 
0,40 6,49. 10? 0,150 1.38 - 10° 0,248 0.375 7.94.10? 0.166 1.51- 10° 0,20 
0.425 5.20-10? 0.159 9,41-10° 0.254 0.40 6.71-10? 0.187 1.07:10° 0.257 
0.45 4.58 - 10? 0,184 6.38 10? 0.252 
ÄA= 
6) 
1.20 1.17- 10% 0.040 
0.25 1.11- 10° 0.062 
0,30 1.S1- 103 0,0S6 7.46 - 10% (0.189 
0.35 1.03 - 10° 0,114 2.76. 10° 0.225 
0.40 6.22-10° 0,150 1.28- 103 01.253 
0.425 5.29.10? 0.173 8.6310? (1.262 
0.45 4,92. 10° 0,218 5.44 - 10° 0,249 
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Fortsetzung zu Tabelle 3: 





a : 1 
NK 0) I; 10 RK 0 





l 






( | m | | ao” = | | (q Ö* ia r | — ) (a Ö*F), | ) (a Öö*F), 



































>20 6.16- 10° 0.006 1.08. 10? 0.187 (1.20 8.283 103 0.179 1.14 - 10% 0.199 
0.30 1.45 - 103 0,139 4.009.108 0.242 0,30 1.20 - 10% 0.210 2.72.10? 0.257 
0.375 6.47: 102 0.195 1.33 103 0.286 0.325 1.35 - 10° 0.228 1,04 - 10° 0,265 
0.40 5.49: 10? 0.214 975-102 0.296 0.35 1.109.103 0.259 1.50. 10° 0,270 
(1.4?5 1.65: 102 0.251 6.65 - 102 0.300 t 
K Z 
K 1 +) 
\ 
he) 0.20 458.103 V.IS8 1.04-10* 0.195 
> 712-103 0.141 1.49. 10% V.ISS 0.2D 3.82 10° 0,203 2.05 - 103 0,225 
0,30 1.97 10° 0.183 3.10.10? 0,256 0,30 1.91- 10° 0.224 2.58 103 0.252 
0.375 7.8610? 0.251 1.07:-103 (>88 0.325 1.55 10° 0.246 1.76- 103 0.262 
0,39 1.66 - 10? 0.272 8.45: 102 (1.284 
j 1 
K —— 
)+) 
0.25 1.10. 10% 0.212 1.70. 103 0.224 
0,30 206-103 0.235 237-103 0,250 
1 ‘ 1 
N ) L. — RK 
.) 40 
U D\ IPRE Un 6*\ 5 l OF \ 
C pe. | (a 0®), | ao0®), C Bi ao0®), | iur (a 0*), 
v 4 v 2 5 v N v 2 ? 
20 5.26: 103 0.136 1.30. 10° 0.218 0.30 1.77-10? (0.247 > 26.103 0,276 
0.30 1.32.103 (1.157 3.58: 103 (0.269 0.325 1.27: 10° 0.256 1.78 103 0.294 
V.3TD 3.92.10? 0.212 ‚21-10? 0.310 0.35 426-100? V.268 1.43 - 103 0.308 
0.40 1.34 - 10? 0,231 401-107 0.320 0.37 1.92.10? 0.290 1.04- 10° 0.319 
0.45 287.10? 0.272 5.05.» 10? (0.340 ' 
h IN = 
. di 
IN 
u 0.30 1.88. 103 V.2D8 217-103 VD80 
0,25 | 3,88. 10° 0,216 | 4,55. 10° 0,243 0.325  1,35-10° 0.266 1,71-109 0.296 
0,50 1.52. 10° 0,212 2 0.35 1.01: 103 >80 
0.35 S.65- 10? 0.242 
1.375 .6.76- 10? 0.2509 1.014. 70° 0.317 K= 
.),- 
0.40 5.39.10 0285 817-102 0.321 Zn 
.) JN» 3 97 
0.425 5.06- 10° 0.311 5.82.10? 0,326 0,50 2,0210 Ze 
0.325 1.140: 103 0.274 1.65 - 10° 0,204 
0.35 1.04- 103 0,2090 1.26 - 10° 0.312 


holisehen Bereich 2” — eonst. nehmen. sondern man hat U’ /U’ nach der genauen Blasıusschen 
Reihe zu nehmen. Die ersten Glieder ergeben: 


A, 


3 
1.68 I; 3,65 | R Ä 





Un | 

(43). 
u” ./Iy“ 
Di 18,4 | rn) 


Damit sind alle erforderlichen Angaben gemacht, um die numerische Reehnung durch- 
führen zu können. In Abb. 2 ist für einen bestimmten e-Wert (e= 0,10) das oben besprochene 
Polardiagramm wiedergegeben, wobei die verschiedenen E-Kurven den angeschriebenen Werten 
der Schichtungsgröße AK entsprechen und immer L=0 ist. Man sieht, daß es für diesen 


3 ou PA; 
c-Wert etwa für K >, keine Schnittpunkte der E-Kurven mit der -Kurve, d. h. keine 
_* jw 
Instabilität, mehr gibt. Für andere e-Werte ist dieser größte K-Wert, für den noch gerade 
eine Instabilität möglieh ist, ein anderer. Der überhaupt größte A-Wert, für den noch gerade 
eine Instabilität möglich ist (kritischer Wert von A), ist A == '/,,, wenn man immer L=0 an- 
nimmt. Für andere Werte von L(L= io: + "10, + ",) ergeben sich andere kritische Werte 

von K. 
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Bisher hatten wir alle Längen auf die Grenzschichtdieke Ö bezogen (halbe Breite der dem 
Blasiusschen Profil eingeschriebenen Parabel, Abb. 1). Bei der Darstellung unserer Rechnungs- 
T. 


ergebnisse wählen wir die Verdrängungsdicke * fi i dy als Bezugslänge (* = 0,341 0). 
0 

Die von uns berechneten Durchgangspunkte bei den verschiedenen Schichtungen (charak- 
terisiert durch A und ZL) sind in Tabelle 3 zusammengestellt. Die zugehörigen Kurven sind 
in Abb. 3, 4, 5, 6 gezeichnet. Das von den Kurven umschlossene Gebiet ist jedesmal der 
Bereich der ıinstabilen Störungen. Die Diagramme geben uns den gesuchten Zusammenhang 
zwischen der Störungswellenlänge = 2r/a, der Reynoldsschen Zahl R und der Schichtung 
(K und 1), wie wir ihn in der Einleitung in Gl. (4) formuliert haben. Dabei entspricht jedes 
Diagramm einem festen Wert von L(L= — '/,.,; 0, + '/, + '/,), während die einzelnen Kurven 
jedes Diagrammes zu verschiedenen Werten von K gehören. Jede Kurve ist also charakterisiert 
dureh zwei Parameter AK und Z. Die Kurve K 0, L=-0 ist die von Tollmien berechnete 
Indifferenzkurve für die homogene Flüssigkeit. Für stabile Schiehtungen (K >0, L>V) ıst 
der Bereich der instabilen Störungen kleiner und für instabile Störungen größer als für die 
homogene Flüssigkeit. 

Dabei weisen die Indifferenzkurven noch die folgende Eigentümlichkeit auf: Für die 
homogene Flüssigkeit und für instabile Schiehtungen erstreckt sich jede Indifferenzkurve für 
oroße Reynoldssche Zahlen ins Unendliche, während für stabile Schiehtungen die Indifferenz- 
kurven geschlossen sind. Mit zunehmender stabiler Schiehtung wird der von den Indifferenz- 
kurven eingeschlossene instabile Störungswellenlängenbereich immer kleiner, bis schließlich 


e//=-- ef 














07 
2 3 I 40 4 50 
Abb. 3. Die Indifferenzkurven für die Plattenströmung mit Diehteschiehtung. Die reziproke Störungswellenlänge « 0* 
de 
als Funktion der Reynoldsschen Zahl und der Schiehtungsgeröße X, für ZI ". 
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Abb, 4. Die Indifferenzkurven für die Plattenströmung mit Diehteschiehtung. Die reziproke Störungswellenlänge « d* 
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- . Um 0* a 5 a 
als Funktion der Reynoldsschen Zahl " und der Scehiehtungsgröße A, für 1 TR 
1} 
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\bh. 5. Die Indifferenzkurven für die Plattenströmung Abb.6. Die Indifferenzkurven für die Plattenströmung 
mit Diehteschichtung. Die reziproke Störungswellenlänge mit Diehteschichtung. Die reziproke Störungswellen- 
ı)* T ) 
* als Funktion der Revnoldssehen Zahl undder länge « d* als Funktion der Reynoldssehen Zahl " 
| ne 
Schiehtungeseröße K,für 1 Ki, und der Schichtungsgröße A, für L a: 
) IH 
VI na iu: 33 I "1 
v /H III %5 | | 
I l 
| N H+41 \ 
\ | 
‚ 1} 
ITEM 
| / 
| 
|, / | 
A | 
' |; Tb | 
} 
ı 
IA 
7 
; | 
HH | 
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4 “ f ' } 
eu | 
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ee ZN l 4 
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| U) nl 
— in _ mn Vesmumenne .- a mw i 1 EEE 
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Abb. 7. Die kritische Reynoldssehe Zahl der Platten- 
strömmng als Funktion der Richardson-Zahl und der 


Um? 
Froudeschen Zahl AR, Rıie, =]. 


)* ( 


hei einem noeh von L abhäneieen Wert von K die Indifferenzkurve sich auf einen im End- 
lichen in der « R-Ebene gelegenen Punkt zusammenzieht. 
Hier interessiert nun insbesondere die zu eimem bestimmten Schichtungszustand ge- 


’ ER IP 5 
hörige kritische Revnoldssche Zahl. In Abb. 7 ist die kritische Reynoldssche Zahl( = 


In 
Damit ist also die für den Vergleich mit dem Experiment wichtige Beziehung gefunden, die 
wir in der Einleitunze in Gl. (6) formuliert hatten. Für kleine Werte von © ändert sich der 


als Funktion von © aufeetraren. und zwar für verschiedene Werte der Froudeschen Zalıl 
( 


” " ” ” . . r / . « v ® .. 
Verlauf der R,-Kurven nur wenig mit der Froudeschen Zahl je, . Jedoch die Stabilitäts- 
\ ( 4 
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orenze 9,, also der größte Wert der Schiehtungsgröße, bei dem noch Turbulenz möglich ist, 

| UN... ; 

hängt noch etwas von der Froudeschen Zahl ab. Die Abhänegiekeit 9, 9, | ,) ist ın 
( { 


IR 
Sy 
tung gegenüber ihrer Schwerewirkung zu vernachlässigen ist, erhält man ©, 0,0409. Es 
ereibt sich also bei Berücksichtigung der Reibung für das Erlöschen der Turbulenz ein 
sehr viel kleinerer Wert der Schiehtungsgröße als ohne Reibung nach G. I. Tavlor und 
(‚oldstein, die 9.0250 erhielten. 

Nach Gl. (Ita) ist die Schiehtungseröße AK so definiert, daß man das Geschwindiekeits- 
refälle an der kritischen Stelle zu nehmen hat, das je nach der Lage der kritischen 
Schieht verschieden ist. Bei der Angabe unserer Reehnungsergebnisse haben wir es jedoch 
vorgezogen, die Riehardsonsche Zahl mit dem Geschwindiekeitsgefälle an der Wand zu 
bilden, nämlich 


Abb. S wiedergegeben. Für 0, also für den Fall, daß die Trägrheitswirkung der Schich- 


4 do, AdUN\? 
() — > 
ee, ee 
die wir aus AK dureh die Umrechnung 
AT 2 AU > 
(-) u i 
[N | dy j dy 2 


AT! \/dU } , i , 
| immer nahezu eleieh I ıst, sind die 9-Werte nur wenle von den 


dy/,/\dy)ıa 
h-Werten verschieden. 


erhalten. Da | 


Im nächsten Abschnitt geben wir einen Vergleich unserer theoretischen Resultate mit 
Versuchsergebnissen. 
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Abb. 8. Die Stabilitätsgrenze 9, Abb. 9. Eine Geschwindigkeits- und Tempe 


raturverteilune im Göttinger Warm-KRalt- 
luftkanal (nach H. Reiehardt). 


als Funktion der Frondeschen 


ii . 
Zahl| ,. .)- 
! 4‘ 


$ 9. Vergleich mit Versuchsergebnissen. Geschwindigkeits- und Temperaturverteilungen 
in einer Strömung mit Diechteschiehtung, die für den Vergleich mit unserer Theorie sehr 
ceeienet sind, hat schon vor einigen Jahren H. Reichardt®) im Warm-Kaltluftkanal des 
Kaiser Wilhelm-Instituts für Strömunesforschung zu Göttingen ausgeführt. Der Kanal ist 
I5 m lang, I m breit und 25 em hoch. Die Geschwindigkeitsprofile wurden am Ende des 
Kanals mit einer Hitzdrahtsonde gemessen, die zur Kompensation des Temperatureinflusses 
aus zwei Hitzdrähten mit verschiedenem Durchmesser besteht. Der eine von diesen Hitzdrähten 
konnte gleichzeitig zur Temperaturmessung benutzt werden. Untersucht wurde bisher nur 
ler stabile Fall, wo die obere Platte des Kanals mit Dampf geheizt (100°) und die untere 
mit Leitungswasser gekühlt wurde (10° C), so daß man ein Temperaturgefälle von WO" C auf 
>> em hat. Es mußte mit sehr kleinen Windgeschwindiekeiten von etwa I m/s gearbeitet 
werden, da bei größeren Gesehwindigkeiten ein merklieher Einfluß der Temperaturschiehtung 
auf die Strömung nieht zustande kommt. In Abb. d ist eine gemessene Gesehwindigkeits- 
und Temperaturverteilung wiedergegeben. Man erkennt, daß im mittleren Teil des Kanals, 


6) Bisher noeh rieht ausführlich veröffentlieht, vgl. jedoch die Versuchsberiehte im Lit.-Verz. Nr. 26. 
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Vergleich mit Versuchsergebnissen. In Über 
Theorie wurden im stabilen 
instabilen Bereich nur 


Abb. 15. 


einstimmunge mit der 
Bereich nur laminare und im 
turbulente Zustände beobachtet. 


(Geschwindigkeitsprofile). Die gestrichelten Kurven wurden 
legt. Übereinstimmung zwischen Messung und Rechnung 


instabil“ mit „exper.: turbulent‘*. 
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wo die Geschwindigkeit nahezu konstant ist, auch die Temperatur nahezu konstant ist. Es 
ist also die in unserer Theorie gemachte Voraussetzung, daß eine Diehteschiehtung nur in der 
(renzschieht vorhanden ist, und daß die Dichte außerhalb der Grenzschiecht konstant ist, 
ziemlich gut erfüllt. 

Für den Vergleich mit unserer Theorie haben wir die gemessenen Profile dahin aus- 
oewertet, daß wir für jedes Profil die maximale Geschwindigkeit U, und die Grenzschicht- 
dieke $ (und damit die Verdrängungsdicke $*) der Grenzsehieht ermittelt haben. Mit diesen 
wurde dann die theoretische Geschwindigkeitsverteilung der Laminarströmung ermittelt und 
deren Gescehwindigkeitsgradient an der Wand 

dU BI 19 
| 1.68 
Iy/ ) 
erreehnet. In Abb. 10 bis 14 sind für verschiedene Maximalgeschwindigkeiten 7/,, jeweils für 
die Grenzschiehten an der warmen und kalten Platte die gemessene und die für unsere 
Stabilitätsrechnung zugrunde gelegte Geschwindigkeits- und Temperaturverteilung aufgetragen. 
Der Zusammenhang zwischen Temperatur und Dichte ist wegen des konstanten Druckes 
vegeben durch 
0 2 Ido IdT 
De u oter ody Tdy' 


wo T die absolute Temperatur und der Index w die Werte an der Wand bedeutet. Unsere 
für die Theorie (Gl. (15)) angenommene Dichteverteilung nach dem Gesetz 
Yu 


-— > 
Oy = OO, u ( 


entspricht damit einer Temperaturverteilung 
‚m mn ’4Y - 
/ En e . ‘ ° . . . . . ° . . . . . R ° ° \ I»). 


Die Konstante y wurde aus den Messungen so ermittelt, daß die gemessene Temperaturver- 
teilune mögliehst gut mit der durch Gl. (45) gegebenen Verteilung übereinstimmt, die wir 
unseren Reehnungen zugrunde legten (vgl. Abb. 10, 11, 12, 13, 14). Damit kann man dann 
für jedes Profil die Riehardsonsche Zahl 


qaT 
9 Tdy 
(aD) any 

dy | dy " 


U I, 
die Froudesche Zahl F? _ und die Revnoldssche Zahl R —" ausreehnen. Dabei wurde 
( ‘4 ° y 


für die kinematische Zähigkeit » der Wert gewählt, welcher der mittleren Temperatur in der 
Grenzschicht entspricht. In Tabelle 6 sind diese Zahlenwerte zusammengestellt. 

Die Entscheidung, ob das gemessene Geschwindigkeitsprofil laminar oder turbulent war, 
wurde einerseits aus oszillographischen Aufzeichnungen der Spannungsschwankungen des 
Hitzdrahtes gewonnen, konnte aber andererseits auch unmittelbar aus der Form des Ge- 
schwindigkeitsprofils ersehen werden. Für die lJaminaren Profile stimmt die gemessene Ge- 
schwindiekeit trotz der endlichen Kanalbreite und der Diehteunterschiede recht gut mit der 
Blasiusschen Plattenströmung überein, während die turbulenten Profile nahezu die Form des 
'/-Potenzgesetzes haben. 

/um Vergleich unserer Theorie mit dem Experiment sind in Abb. 15 die Meßpunkte in 
ein R, ©-Diagramm eingetragen, wobei für © ein logarıthmischer Maßstab gewählt ist. Jedem 
Profil der Abb. 10 bis 14 entspricht ein Punkt in diesem Diagramm, wobei die laminaren 
und turbulenten Profile mit verschiedenen Bezeichnungen eingetragen sind. Die ausgezogene 
Kurve ist die von uns berechnete Stabilitätsgrenze, welche die kritische Reynoldssche Zahl 
als Funktion der Schiehtungsgröße © für die Froudesche Zahl F--0 gibt. Die Kurven 
für die anderen Froudeschen Zahlen sind fortgelassen, da sie sich, wie man aus Abb. 7 
erkennt, von der Kurve für F=0 nicht wesentlich unterscheiden. Die auszezozene Kurve 
eibt die Grenze zwischen den theoretisch stabilen (laminaren) und den theoretisch instabilen 
(turbulenten) Zuständen. Außerdem ist der Wert 9,4 0,250 nach Taylor-Goldstein und 
der Wert A, = 575 für die Plattenströmung bei homogener Flüssigkeit eingetragen. Der Ver- 
Un 
4 
zwischen 0,56 und 3,82 liegt (vgl. Tab. 6), ergibt nun, daß im theoretischen stabilen 


sleich unserer Theorie mit den Versuchsergebnissen, bei denen die Froudesche Zahl 
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Gebiet nur laminare Zustände und ım theoretisch instabilen Gebiet mit 
einer einzieen Ausnahme nur turbulente Zustände beobachtet wurden’). 
Dabei ist noch zu beachten, daß es keinen Widerspruch gegen die Theorie bedeutet, wenn 
nahe oberhalb der Stabilitätserenze R,= R, (O9) noch laminare Zustände beobachtet werden. 
Diese Kurve bedeutet nämlich diejenigen Zustände, wo die Anfachung der Störungsbewegung 
verade gleich Null ist. Für ein wirkliches Eintreten der Turbulenz müssen die kleinen 
Störungen jedoch beträchtlich angefacht sein, was erst bei etwas höheren Revnoldssehen 
Zahlen der Fall ist (vgl. hierzu unsere frühere Arbeit |22]). Es besteht also völlige Über- 
einstimmung unserer Theorie mit den Versuchserzebnissen. 

Der von Tavlor und Goldstein errechnete Wert 9,4 == 0,250 liegt weit entfernt von 
denjenigen Werten 9, wo nach den Reichardtschen Messungen die Grenze zwischen 
laminarer und turbulenter Strömung liegt. Dabei ist allerdings zu berücksichtigen, daß Taylor 
und Goldstein ihren Stabilitätsuntersuchungen andere Geschwindigkeitsprofille zugrunde 
legten, als sie hier im Experiment verwirklicht sind. Unsere Theorie dagegen, deren Voraus- 
setzuneen bei den Messuneen recht gut erfüllt sind, ereibt nach Abb. 15 eine recht be- 
(riedigende Übereinstimmung mit dem Experiment. 

$S 10. Zusammenfassung. In dieser Arbeit wird mit Hilfe der Methode der kleinen 
Schwingungen die Stabilität einer geschichteten ebenen Laminarströmung untersucht (= und y 
reehtwinklige Koordinaten in der Ebene; .- Riehtung = Hauptströmungsriehtung, %- Riehtung 


senkrecht nach oben). Die zugrunde gelegete Laminarströmung UT (y) ist eine „Grenz- 
schiehtströmung*, die vom Wert 7 =0 an der Wand zu dem für große Wandabstände kon- 
stanten Wert U= U, ansteigt. Kine Dichteschiehtung wird nur in der Grenzschicht 


(Dieke = d) angenommen, während außerhalb der Grenzschicht, im Gebiet der konstanten Ge- 
schwindiekeit U,,, die Dichte als konstant vorausgesetzt wird. Für die Stromfunktion einer 
Partialschwingung der überlagerten periodischen Störungsbewegung setzt man: 


X ch) 


ft 


i(aX it 
vr, D=giy)e —=g(y)e 


Dabei ist @ immer reell und bedeutet die räumliche Kreisfrequenz der Störung 
/=2a/a= Wellenlänge der Störung). p ist im allgemeinen komplex, und das Vorzeichen des 
Imaginärteils von 5 entscheidet über die Stabilität oder Labilität der Störung. Die Stabilitäts- 
untersuchung ist ein Eigenwertproblem der sich für die Störungsamplitude g ergebenden 


Störungsdifferentialeleichung: 


-apl-V-JUpg+Kpg+LW-gGI-W-VPYy+Ugl] 


Bi ec) |q rm I? 7 ’ nu a" l; (4 2 Y')| 


{ 
dA R 
mit den Randbedingungen 9 9’=0Oander Wand 0, und in großem Wandabstand y=x. 


0 d 0.40 NV d Oo 


R bedeutet die Revnoldssche Zahl: AK ‚ „(Rıehardsonseche Zahl) und L = 
y O dy [ m“ 0 d „ 


L 


sind dimensionslose Schichtungsgerößen, | ra Um] dg Ist eine Froudesche Zahl. Für eın 
\ 


vorgegebenes Wertequadrupel a, R, K, L hat diese Gleichung nur eine Lösung für einen be- 
stimmten Wert von e. Wir beschränken uns auf rein reelle Werte von e, bestimmen also bei 
vorgegebenen Werten von K und L im a R-Diagramm die Kurve der Durchgangspunkte von 
Stabilität zu Labilität (Indifferenzkurve). Es werden stabile (K >00, L>OV) und ıinstabile 
Schiehtungen (N <0, L<0O) untersucht. Als Resultat ergibt sich, daß bei festgehaltener 
"roudescher Zahl mit stärker werdender Schiehtung (wachsendes K) die kritische 
heynoldssche Zahl erößer und der Bereich der instabilen Störungen in der a R-Ebene 
kleiner wird, bis schließlich bei einem kritischen Wert der Schichtungsgröße K, der noch von 
der Froudeschen Zahl abhängt, völlige Stabilität der Strömung für alle Störungswellen- 
längen und Reynoldsschen Zahlen eintritt. 

Die Rechnung wird durchgeführt für die Strömung längs einer ebenen Platte, für die 
schon früher W. Tollmien die Indifferenzkurve für homogene Strömung berechnet hat. Für die 
kritische Schichtungsgröße O9; : g« | ‚) ‚ für welehe die Turbulenz völlig erlöschen 

odyf\dy/e 


> 


muß, ereeben sieh Werte von 0,0409 bis 0,029 für Froudesche Zahlen T 0 bis 5, 
( 
während G. I. Taylor und S. Goldstein durch eine ähnliche Rechnung, bei der die Flüssig- 


”), Bei dem herausfallenden Profil ließ sieh aus den oszillographischen Aufzeichnungen keine eindeutige Ent- 
scheidung zwischen Jaminar und turbulent treffen. 
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keitsreibune und Krümmune des Profils vernachlässiet wurden, ©, !/ erhalten hatten, und 
l. A 


nach rohen Abschätzungen von L. F. Richardson und L. Prandtl sich 9, =1 bzw. 9, —=2 
ergeben hatte. 

Die wesentlichen Voraussetzungen unserer Untersuchung (Plattenprofil, Diehteschichtung 
nur in der Grenzschieht) sind recht gut erfüllt bei Messungen von H. Reichardt im 
Göttinger Warm-Kaltluftkanal, bei denen die obere Platte mit Dampf geheizt und die untere 
mit Leitungswasser gekühlt wird (stabiler Fallı. Die Entscheidung, ob ein gemessenes Ge- 
schwindigkeitsprofil Jaminar oder turbulent war, wurde aus oszillographischen Aufzeichnungen 
der Spannungssehwankungen eines Hitzdrahtes gewonnen. Der Vergleich der Messungen mit 


unserer Theorie ergibt eine recht befriedigende Übereinstimmung. >40 
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Die Bestimmung derKnicklast gedrückter, rechteckiger Platten. 
Von E. Trefftz in Dresden '). 


welche am Rande durch Druckkräfte in der Plattenebene belastet ist, ist mehrfach 
behandelt worden; in dieser Zeitschrift z. B. durch Herrn Taylor (The Buckling Load for a 
reetangular Plate with four elamped Edges, diese Zeitschrift, Bd. 13, 1953, S. 147 bis 152) 
und neuerdings durch Herrn Fax en (Die Knickfestigkeit rechteckiger Platten, diese Zeitschrift, 
Bd. 15. 1935, S. 265 bis 277). Diese Arbeiten werden durch die folgenden Ausführungen 
insofern ergänzt, als es mir gelingt, nachweislich untere Grenzen für die Knieklast anzugeben. 
Da obere Grenzen stets mit dem Ritzschen Verfahren ermittelt werden können, kann die 
Knicklast auf diese Weise in praktisch ausreichend enge Grenzen eingeschlossen werden. 
Ich stelle meine Methode am Beispiel der eingespannten, quadratischen Platte dar, auf deren 
Ränder der allseitig gleiche Druck p wirkt. An diesem einfachen Beispiel wird die Methode 
so deutlich, daß eine Darstellung für den allgemeineren Fall nichtquadratischer Rechtecke 
bzw. verschiedener Drucke an beiden Ränderpaaren sich erübrigt. 


D' Problem, die Stabilitätsgrenze für eine eingespannte, rechteckige Platte zu bestimmen, 


1. Die Bestimmung der Knicklast als isoperimetrisches Variationsproblem. Es sei a die 
Kantenlänge, k die Dieke der quadratischen Platte, £ der Elastizitätsmodul, m die Quer- 
Eh’ m? 
12 (m? — 1) 
Flächeneinheit der eingespannten Ränder, bei welchem das Ausbeulen der Platte eintritt. 


kontraktionszahl des Materials, N = der Plattenfaktor; ferner p der Druck auf die 


Wir legen ein £»7- Koordinatensystem so in die Platten-Mittelebene, daß die &-Achse in 
die untere, die n-Achse in die linke Kante des Quadrates fällt. Dann genügt «die beim Aus- 


beulen der Platte eintretende Verschiebung »v der Differentialgleichung 


‚jw „ dw el (O’w oO’) 


N = u Er h oo po (dh) 
10 &* Ro oO) Porto 7”) 
mit den Randbedingungen 
Om | 
»„=0 und —=0O für E=0 und !=a (2) 
O8 
und 
Ow 
w—=0 und =O ir n=0 a He... re 
ol) 


Es wird nach dem kleinsten Wert von p gefragt, für welchen das dureh die Gl. (1) bıs (9) 
gegebene, homogene Randwertproblem eine von Null verschiedene Lösung hat. 
Zur bequemeren Rechnung setzen wir 


” ax dy 
= R Ks . : . . i r . > . R ; i j i ( 1) 


JT JT 


und erhalten mit den Abkürzungen 


e „.phe 5) 
u tume »: er » 
0x Oo? N n? 
die Differentialgleichung 
Bu > ae . re 5 A 
mit den Randbedingungen 
O0 Ww ni . 
w=(0 und = ME u u. . rt ih 
Öx 
und 
Om - | 
w--0 und =0 für y=-0 und y=n....:.:.:..RD). 
0y 


Dieses Eigenwertproblem ist mit dem isoperimetrischen Variationsproblem identisch, 
unter allen zweimal differenzierbaren Funktionen, die den Randbedingungen (7) und (8) 
genügen, diejenige zu ermitteln, welche bei der Nebenbedingung 


t), Vorgetragen in der Mathematisehen Gesellschaft in Göttingen am 5. 2.35. Zusatz bei der Korrektur: Herr 
Willers (Freiberg) hat mieh darauf aufmerksam gemacht, daß Herr Weinstein im Journal of the London Math. Society 
(Bd. 10, 1935, S. 184) eine Arbeit „On a Minimum Problem in the Theory of Elastieity“ veröffentlicht hat (eingereicht 
am 14. 2. 35), in der er das gleiche Problem behandelt. Herr Weinstein geht ebenfalls von dem Courantschen Prinzip 
aus und gelangt auf etwas anderem Wege als ich zum gleichen Ergebnis. 
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dedy ] i : : (9) 


das Integral 
MM ARE ı = 3 a a a . (10) 








zum Minimum macht. (NB. Die Doppelintegrale sıind hier und ım foleenden über die Fläche 
des OQuadrats zu nehmen.) 


"ormuliert man das Eigenwertproblem auf diese Weise, so erscheint der Eigenwert 4? 
als Laerangescher Faktor für die Nebenbedineune /=1. Soll nämlich bei der Neben- 
beilineunge .J | das Interral M ein Minimum werden, so hat man nach der Vorsehrift von 
l,arrangze den Ausdruck 


i ..E [Wo w\? 0 ımw\? 
n=M-2IJ=\\(Amf®drday-a\\ I) Hl )Naray..... cl, 
1; | Now oy/) | 
in welehem zu dem Minimalinteeral \M/ das mit dem Laeraneeschen Faktor / multı- 





plizierte \ebenbedinenunesinteeral ./ hinzueefüet ist. ohne Rücksicht auf die Nebenbedineune 
zum Minimum zu machen. Setzt man die Variation 9 L eleiech null: 









N, >\\ EI dr dady 27 \\ 


. . 


(0 woNnMm RETEen)) m) 





lv dy 0 i : r (12), 







( 
Or 0 oy 04) 


so erhält man dureh partielle Inteeration unter Berücksiehtieune der auch für Se zeeltenden 
Randbedineunzen (7) und (8) 








nes owlA4Aw+ KR AnNdzdy . -» :» : 2 nu ir. KM. 









Soll AL für alle dw verschwinden. so muß 
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sein, Damit ist die Identität des Lagrangeschen Faktors /#° mit dem Eıgenwert der 
Differentialeleiechung «darzeetan. Wiehtie ist noch «die folgende Bemerkung: Setzen wir ım 
(Gl. (12) für das beliebige 18 speziell dw = em, so wird oL=2eh. Da oL=V0 sein muß, 
so folet L=0,. d.h. der Minimalausdruek 17 verschwindet, wenn man in den Integralen die 
lösung = des Vartationsproblems einsetzt. Damit ergibt sich 

M 


I; M-32J=0 bzw. #4 Bra ee (19), 






















d.h. der Eizenwert 2° wird. da .J | ist. eerade eleich dem Wert des Minimalinteerals 4. 





2. Abänderung des Variationsproblems. Verkleinerung der Eigenwerte. Das Variations- 


problem zur Bestimmung der Knicklast ıst für die eingespannte Platte nicht ın einfacher 






N 438 
r . 5 N 0. . > . Ban . ® ch ’ "ur 
Weise Jösbar. Für die frei gelagerte Platte, für welche die Randbedingung | 0 fortfällt, 
0 7/ 
ist die elementare Lösung bekannt. Die Eigenfunktionen werden wa, Cu. sin ar sin v9, 
die Eigenwerte Zu,» a +r”, wo u und » beliebige ganze Zahlen sind. Der kleinste 








KEigenwert ergibt sich mit ur =1 zu #2. 














Um für die eingespannte Platte untere Grenzen für 4° zu bekommen, benutzen wir den 
(Gedanken von Courant., «die Bedinzuneen zu mildern, die der gesuchten Funktion w auferlegt 
sind. Dabei sinken die Eizenwerte,. Gelingt es also, die Änderung des Vartationsproblen:s 
so einzuriehten, daß nach Milderung der Bedinzeungen das Problem numerisch lösbar wird, 
so hat man untere Grenzen für 2° gewonnen. Wir verfahren folgendermaßen: Wir halten 
daran fest, daß die gesuchte Funktion, die wir jetzt v(., y) nennen, am Rande verschwindet, 
0 


wir verlangen aber nieht mehr, daß an jedem Randpunkte 0) ıst (Einspannung), sondern 


on. 


ersetzen diese Randbedineunze dureh die schwächere Forderung, daß 


s 
: 0% . ' 
ürr z2=0 md z=n Km snydy=V . .. x... . (16), 
x yar 

m 

R 
) | E "si / 0 17 
tür / ( une / 7 sın.r (dx Pe + 

J J ey 
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sein soll. Wir suehen also unter allen am Rande verschwindenden Funktionen, die den 


Bedineungen (16) und (17) genügen, diejenige, welehe bei der Nebenbedingung 


do r\? N9»\? 
| 


| 
\.) au) dedy-| 
das Integral 
M=$y((de’daedy 
zum Minimum macht. 


Ist "(.r, y) die Lösung «des gemilderten Problems, »*(., y) die Lösung des ursprünglichen 
Problems und sind Mir) (de? daedy und Mi) (Te) daedy die gesuchten Minimal- 
werte des Integrals M in beiden Fällen, so ıst sicher 


En 5 


Denn die Funktion #w zenügt den Bedingungen des zemilderten Problems: sie ist also ent- 
weder Lösung auch des gemilderten Problems, dann ist v(@, = mw(a,y) und in (18) steht 
las Gleichheitszeiehen: oder »e ıst nieht Lösung, dann ist Moe) nieht Minimalwert des 
oemilderten Problems, der wahre Minimalwert Mir) heet also unter Me) und in (15) steht 
das < Zeichen. Da die Minimalintegrale mit den KEigenwerten identisch sind (siehe Abschnitt I), 
so folgt also, daß der Eigenwert des gemilderten Problems tiefer liegt als der des ursprünglichen. 


3. Lösung des gemilderten Variationsproblems. 
Wir setzen die Lösung ın der Form 


Das abgeänderte Variationsproblem ist in 
einfacher Weise lösbar. 


"(w,y)=  Nawsinursinvy. . | a, . + (9 


an. womit (gleichmäßige Konvergenz vorausgesetzt) die Bedingung erfüllt ist, daß e (we, yp) an 
den Rändern des Quadrates verschwindet (zw und » durehlaufen alle ganzen Zahlen). Die 
voemilderten Randbedinzungeen (16) und (17) verlangen, daß für 


.- \IX\ 


tAursCOSUZSINDTrU . . .  . 80), 





a das Glied mit siny und für 


(\ y' \ 7} \ a} 
07] ACT A 


bei. 0 und .r 


’omvBINMECOSIYU . .: 2: 2 2 2 0. 0. Kl, 


bei y=0 und == das Glied mit sin.e verschwindet. Das gibt die vier Gleichungen 


y y 
Nudaıı ), Nr dueos ua 0", 
— | | | . 


\'ydtıy vd. 


\'yaı»cosrvr rd, 


die man auch, wenn # die ungeraden, 9 die geraden ganzen Zahlen bedeuten. in der folgenden 
Form schreiben kann: 


5 a 
Mad V u En : a er (22), 
2 u 
Hd 
\ a 
ut (du 0 > D 33) 
nz u 
„Hl 
\' 7° 
Jayız V eu irn 
En — 
( 
Yigg 9 Erz Eee KR 
nz ut 
4 
ee 2” a i 
(Die Größen ‚, usw. sind Lagrangesche Faktoren, die nachher benutzt werden.) 
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Das zum Minimum zu machende Integral M wird 


M ( Ir” dr dy 7 \ \ m. 


das Nebenbedingzungsintegral 


Furlu? + vr?) a 


Jydeay N \ ( 


Y E A 


Die Koeffizienten «a. , sind nun so zu bestimmen, daß mit den Nebenbedingungen (27) und (22) 
bis (25) M ein Minimum wird. Zu diesem Zweck fügen wir zu dem zum Minimum zu 
machenden Ausdruck M die sämtlichen Ausdrücke, deren Werte vorgeschrieben sind, mit 


lagrangeschen Faktoren multipliziert, hinzu: 


tlg 





“ + y 
‚) \ 
» 
nz 
‘d] 


und bestimmen die Koeftizienten «a., aus den Gleichungen 


Damit erhalten wır 


und wenn «>| und 


Pr en MN ao .)%) 
- v)(u — »r° % 1 VE 4). 


Setzen wir die Werte der a, aus (30) bis (32) in die Bedingungsgleichungen (22) bıs (25) 


ein, so erhalten wir zur Bestimmung der Lagrangeschen Faktoren a, pP, a’, 5 die Gleichungen: 


A 


y u? 
\ 7 a ‚ A834). 


?+-1)wW +1 


g 
tlg’ +1 


" 
+’ —| 


Das sind homozene. lineare Gleiehuneen. die nur von Null verschiedene Lösungen haben. 
wenn Ihre Determinante verschwindet. Das Verschwinden der Determinante bestimmt den 
Wert von 7°. Es bestehen die folgenden Lösungsmöglichkeiten: 
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BR | ’ u? 
l. a p beliebie. a" pP’ v, rn — \ - - eu .„ „- ; KB. 
2 (2 — 4) <<. (w Iu+|l ist 
(26), . 
u 2 | ’ u” 
Il. a P beliebig, a’: B: ), DZ — \ r - —- —(V) ,„ . .(39. 
22 —4) Z_. (W+-l)(W-+1 4”) 
u 
: er 2 ER ’ [Ti 
om Il. «=p=0, .« beliebig, pP’ beliebig s =- =D , „20% 5 AM. 
(24) u )( +1 4”) 
‘ 
nd (22) In den Fällen I bis III wird aw»=0, wenn «>1 und vr>1 ist. 
um zu Hierzu kommt als vierte Möglichkeit 
d, mit 
IV. a=$=a=p=0, alle au, verschwinden außer für | 
” , , | Ve: (tl). 
u=h>1, v=k>1, ?=W®+M>8 
Den kleinsten Wert für /? erhalten wir ım Falle I. was hier nicht ausführlich vorgerechnet 
werden soll, aber leicht zu bestätigen ist. In der Gleichung für 4 läßt sich die Summe zerlegen 
(28) 
u u? nr ”—1 \’ | 10 
<< _ ((W’+-V)W—+] er, Fire a A" DE 4 
[71 „H [Z1 
Mit der Formel für die Partialbruchzerlegung des Tangens 
nz 4z “7 Ä 
ig = > _ en Eee 
(29). 2 Tan U? — 2° 
[7 
erhält man hieraus 
(30). ° u” 3.23. #°-]I 7 ayi® — 1 
gr z hz z Er er Fr tg ns ME 
<<. (wW+-l)(W"+1—/) #4 2 “ 4yW”—] 2 
„H 
> 
(Sl), Damit wird Gl. (35) nach Multiplikation mit 24° 
A ıYy/ i,.ay# —I R 
s=— 14.7 = ty £ a TE 
(32) "2 292 2 2 
Diese Gleichung hat als kleinste Wurzel 
i..’=B518 ee 
.)°) i u 2 3 u. “= . . 
(33). die also tiefer liegt als der Eigenwert des ursprünglichen Problems der eingespannten Platte. 


1s (25) . a E92 
BER 4. Verfeinerung der Methode. Es liegt nahe, die hier angewandte Methode in folgender 
u: Weise zu verfeinern: Das ursprüngliche Variationsproblem wurde oben dadurch abge- 
A le ra ÖWw _ 
ändert, daß die Einspannungsbedingung pr (16) und (17) ersetzt wurden, 
nen 


welche verlangen, daß auf den vier Rändern die Glieder mit sin. x bzw. sin y in der Fourier- 


0 dureh die Gl. 


(4), 
Or De 2 
entwicklung von _ _ verschwinden. Verlangt man nun, daß nicht nur diese Glieder, sondern 
On 
(39). Ze \® ' | 
die ersten h Glieder der Fourierentwicklung von \, verschwinden, daß also für o Ibis h 
Ä y [ 
i ( R 
.);* . . pn 
(36), bei =0 md z=n \ et . 5 ee A 
[Our ne ; 
M 
dv 
12 Yumı . . 
(4) bei y=0 und yenr \ Frhr Koh un Er 5 
| p J 
u) 


ist, so sind die Bedingungen für die Lösung des Variationsproblems immer noch gegenübhe:ı 
aben, . . v: pn . . . : .. 

dem Ausgangsproblem gemildert, so daß die Eigenwerte tiefer liegen als bei wirklicher Ein- 
Aber die Milderung ist schwächer als ın dem oben durchgeführten Falle, man 


den 
Spannung. 








rue 


.. 
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hat sich dem Falle der Einspannung genähert, so daß die Eigenwerte verbessert sind. Die 
Durchführung der Reehnune ist analor zu der Reehnune des dritten Abschnitts. Für h=2 
erhält man 





ı 2 PE 
Ir: — I.) 








(49). 
















5. Numerisches Ergebnis. Um den gesuchten Eigenwert in Grenzen einzuschließen, 
braucht man noch eine obere Grenze für den gesuchten Eigenwert, die am einfachsten nach 
hıtz bestimmt wird. Macht man für die Lösung des isoperimetrischen Problems den die 
Randbedingungen befriedigenden Ansatz 










Very) sin.esın ytasın esıny+ Psin@sinay-+ysindesıny+-osnaszsindy} . (0) 
















und bestimmt «die Koeffizienten a. 5, y. d nach Ritz, so erhält man für 2° den sicher zu 
eroßen Näherungswert /,” 9,52. Als Endergebnis erhalten wir also /,? = 2° = 7,” d.h. 









2.30 >52. Ki: ur ‚ al). 








6. Beziehung zu den Ansätzen von Faxen und Taylor. Wenn man nieht nur, wie es im 
dritten Abschnitt geschehen ist, den Eigenwert des gemilderten Variationsproblems betrachtet, 
sondern die Funktion r (Gr, y) selbst berechnet, indem man die Werte der Koeffizienten «,, 
und «,, (alle anderen verschwinden) aus den Gl. (30) bis (32) in den Ansatz (19) einführt, so 
erhält man 













1’ nz E 

Er. \ u sın a. esın ne \ vsın.esın 09 en 

"ie, dd = e: > ) = — Pr \ (+). ) 
ne +1 7") j (+) +1 u 

„Hl 


bzw. dureh Partialbruchzerleeune 





asınyy) \' vsın ae \ "sin ur] 





(1 7) Ts 


PA inne Hi} 













sin ad \ sın "m \ 
—] 7") “1 j 











lelı verdanke Herrn Faxen den Hinweis darauf, daß man die hier auftretenden Reihen 
summieren kann. Im Intervall von O0 bis 7 ıst 














ur 7 t _ .kan "uesinn.r 
Coikie ——) z bo] , \ 2 

















Damit wıred 











(sivi-2(e 3) sole: 
. )D] A > 0] zB 
Lasın y 2) 2 
"ke, ry) -_ ü 
’ ta N as es u 
| Ro yi—4 bo] , 
(99). 
u = 7 = 7 
2 bo ‚\y 2 bo} I > 
5 sın.r " 2 \' 2 
nA? u a m ch 
Bay — 4 bo} 
Das ist zeerade die Form des Ansatzes, die Herr Faxen für eine erste Näherung braucht 





(diese Zeitschrift, Bd. 15. 1955. 8. 270) und die auch als Sonderfall in dem Taylorschen 
Ansatz (diese Zeitschrift. Bd. 13. 1933, S. 151) enthalten ist. Diese Übereinstimmung ist nicht 
zufällig. Beide Autoren approximieren die Lösung des Eigenwertproblems durch Funktionen, 
welehe der Differentialeleiehunge zenüzen, aber nieht den Randbedingungen. Daß auch ın 
unserem Falle die Funktion * (, ‚) der Differentialgleiehung genügen muß, folgt daraus, daß 
sie als exakte Lösung des gemilderten Variationsproblenms natürlich die entsprechende | 
Ditferentialeleichung erfüllt, wie auch Gl. (55) erkennen läßt.  Entsprechendes gilt auch für 
die nach der verfteinerten Methode (Abschnitt 4) ermittelten Lösungen. >46 
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(49). 


ılıeßen. 
n nach 
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(0) 


her zu 


(1). 
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achtet, 
en Ay, 
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Zum Problem des FErstarrens für den durch 
zwei parallele Ebenen begrenzten Körper. 
Von Aurt Lachmann ın Berlin. 


1. Das Problem. Bei einem durch zwei parallele Ebenen begrenzten Körper ist der 
Übergang vom flüssigen zum festen Aggregatzustand bisher nur für den Fall behandelt worden, 
daß die Anfangstenmperatur des Körpers gleich der Schmelztemperatur ist (Abb. 3). Für die 
Beurteilung technologischer Vorgänge wäre es erwünscht, den Temperaturverlauf auch dann 
verfolgen zu können, wenn die Anfangstemperatur des Körpers oberhalb der Schmelztemperatur 
lieet. Dieses Problem soll im folgenden behandelt werden. 


2, Die Greensche Funktion. Wird die Entfernune von der Iınken Endlläche mit X, 
die Temperatur mit 9 und die Zeit mit T bezeichnet, so lautet die Differentialgleichung der 
Wärmeleitung 

ALz) Aue 
31 Ad at: ’ 5 i (hl), 


wenn «ad  econst die Temperaturleitfähigkeit ist. An den beiden Endflächen soll die Temperatur 
der Umgebung herrschen, welche gleich Null gesetzt wird, ©* sei die unveränderliche Tempe- 
ratur des Körpers zur Zeit T=0 und I der Abstand der beerenzenden Flächen: mit 


X ‚_at u. . 
4 I E ume ! 75 a ae I) 
wird 
PO), a7, PR 
(+). 
Vu of 


Die Randbedineuneen sind 9-0 für #0 und 1, die Anfangesbedinzunge I für 0. 
In den dimensionslosen Einheiten (2) ist die Temperaturleitfähigkeit, der Abstand der End- 
flächen und die Anfangstemperatur gleich Eins geworden: die Schmelztemperatur 9, hegt 
zwisehen Null und Eins. 

Wir gehen von der partikulären Lösung 


(x Sye 


l pP tt r) 
HERIES II u 9-5 | Bu | | . 4 
ayn | f T 
von (3) aus. Diese Temperaturverteilung wird in einem unendlich langen Stab hervorgerufen, 
l JS 
wenn zur Zeit r an der Stelle £ die Wärmemenee e (Dir, £&: t.r)dr=e augenblieklieh zu- 
v4 
eeführt wird. Hierbei ist der Querschnitt zwischen den Endllächen gleich Eins gewählt und 
die spezifische Wärme e==eonst auf die Volumeneinheit bezowen. 
Wird der Wert & ın (4) an den beerenzenden Flächen #0 und x I wiederholt ve- 


spiegelt und werden die entstehenden Ausdrücke immer mit Vorzeichenweehsel addiert, so wird 


() . Er2 nt x +£&)°\ 
| / v „ \ ) v ’ 
var \ 0 Id 7) \ e u EG T (>) 
n 4 H SH 


die Temperaturverteilung in einem endlichen Stab mit den Temperaturen Null an den Enden 


» VO und 1, wenn an der Stelle € zur Zeit r die Wärmemenge lım  SeVtr. sı Krn)de € 
! r—>( 0 
auzenblieklich zugeführt wird. 
Trägt man jetzt der Symmetrie des Problems Rechnung und spiegelt in (5) die Größe S 


an der Mittelebene „= , , so wird 


G(lr,&: tr) =Vle, &; tr) +Ve(er, 1 Er 


} L Y f 4 
N | . (nm A 5) . (2n 1I+X s)? \ . (’n \ 5) 
\ Pe (8 Tr) \ v’ Ss } ve’ If } 
) 4 
] 7] ' ! a — — r 
n % n % n fi (dt) 
FT. 
fZ til r) / 7 

EZ 
N 4 
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die Temperaturverteilung für den endlichen Stab, wenn zur Zeit r bei Zund 1 £ je e Wärme- 
einheiten augenblicklich zugeführt werden. 






wird die Thetafunktion 





Zur Vereinfachung 










mt 


eingeführt; die Greensche Funktion wird dann mit (9) 








Gae,E: ED. = Hd -97 lc +9; rei-—r 14. Er rinnen (RR, 













3. Die Temperaturverteilung. Soll zur Zeit =0 die Temperatur Eins entstehen, so ist 
für die Längeneinheit des Stabes die Wärmemenge e zuzuführen; dies ergibt mit r=0 die 
Temperaturverteilung 








l 


KENSIWNREIBE . a A Ver er A 


u 











Zur Zeit == fiz) sei der Körper zwischen z=z und z=1- z flüssig (Abb. 1): die Kurve 
Fi=\[0° 2°. ,„|) werde als Schmelzkurve bezeichnet. Sind £& und r=fi£) die Koordinaten der 





Schmelzkurve zur Zeit 7° # (Abb. 2) und rücken in der Zeit dr die Schmelzpunkte um d£ 
nach innen, so sind an den Stellen £ und 1 £ die Wärmemengen qdS zuzuführen, wenn 
wieder der Quersehnitt 1 zwischen den begrenzenden Ebenen betrachtet und die Schmelz- 


wärme q== eonst auf die Volumeneinheit bezogen wird: mit 


















fürt =fiz) >0 Pie “ ‚ Schmelzkurve IT 
‚3 En ! } L-fiz, fest 


ı E % 
flüssig fiz) 
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iz re zz 

































wird daher die Erhöhung der Temperatur infolge der freiwerdenden Schmelzwärme 












(f) 
Be.deakTü BENBEE - : : » SL Ei en wen DE He BR 





Mit (11) wird die gesuchte Temperatur an der Stelle = zur Zeit f- fi2) 


IP (t) 
add (a, hd) +9. (, = Sala, E;t,0)ydAE+KS Ga, fe). . . (IA. 


- 


0 u 


Das Problem hängt unter den gemachten Voraussetzungen, daß die spezifische Wärme und die 
Schmelzwärme konstant sind und daß an den Endflächen die Temperatur der Umgebung 
herrscht diese Annahmen sind in Wahrheit nur angenähert erfüllt —, allein von den beiden 
Konstanten 9, und k ab: die Schmelzkurve f=fi(z) ist so zu bestimmen, daß Vz, f) = 9, wird 





> r = 
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Wärme- 
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(9) 


(10), 


I. so Ist 


—() die 


(il). 


» Kurve 
ten der 


um de 
I. wenn 


schmelz- 
‚g 
(12) 
(13) 
(It). 
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neebung 
ı beiden 
". wird 
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4. Die Verteilung der Temperatur, wenn die Schmelzkurve die Gleichung ?- , hat. Ist 
A 


die unveränderliche Temperatur des Körpers zur Zeit =0 gleich der Schmelztemperatur 
und daher 
BE ae Dee et EEE 


so läßt sich die Lösung des Problems in geschlossener Form angeben: es ist für 


| 




















" | r 
)<z2: 5 und <<, 2 : dla. E »\ | ‚|; 
D| 5 
| | 04 (16) 
0<2<s- und .Sı=2 5:0 Bi, 
| Fr u 
= : ı : ’a,d= il —z,D 
worin 
«x 
P(.r) Vn\® du u a re u. 
das Gaußsche Fehlerintegral ist. Die Gleichung der Schmelzkurve ıst 
= fi) — (IS), 
2 
wobeı 4 als 
.) Ad 
u zab\,)e' a er: ER | u + A 
= 
zu bestimmen ist! Abb. 3 zeigt die Temperaturver- 
teilung für 
g* 
2 
Nach (14) ıst 
aVt 
ade let, )ae 
. u 
" 
und daher 
eVi 
G 2,8; 48, s_ dE (0a, la, A). 
Fr A h 
ı 
Mit (19) wird 
a\t er 
. i 5 j 4 A [= Rt 
leles: a )d5 ; ab| Du LAAL ARE) BEN 5 nn (20). 


Mit Hilfe dieser Formel, die für beliebige Werte .r, t und «a gilt, kann die Temperaturverteilung 
auch für den Fall 9. <1 zefunden werden, wenn die Schmelzkurve die Gleichung 


Re Er eG + . Ta 
hat. Nach (14) ıst 


‚ Enzyklopädie der mathematischen Wissenschaften, fünfter Band, Physik, erster Teil, S. 204, 








’ ‚ En Ztschr. f.angew. 
348 ,achmann,. Zum Problem des Erstarrens begrenzter Körper Math. und u: 














\/ 


/ 
Vt 


lır. fi) Va (#.1)-+-k (ilır, 4 
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ı 
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J 
m — 
vr 
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und mit (20) 





> 


Jlır, #) ’), (a. + h ‚ i 7 | £ \e (dla, (a, t)), 





folelieh mit (16) und (21) 










kyr leer +1 - "Vrao[2)er)d,[x,) für <a 
92, m Be e ERS (22), 
ei ’V? 32 a\ — 2 | 
I ad! = je* Hi ‚ a b| > je! dl. e} für > „ 2 

























Für einen Punkt der Schmelzkurve ıst mit .r 


Y 


e) hi J nn D| = lei fi h ' u A D\| z le! \%,|\z, =) 1 
J ( 


dA _ _ = - 





5. Ermittlung der Temperaturverteilung, wenn = f(z) gegeben ist. Nach (14) ist 
Mer. (a, -+HA,(@,t). Der erste Ausdruck ist die bekannte Lösung des Problems für 
den Fall. daß keine Anderunz des Aggreratzustandes stattfindet und daher ?=P0 ist: 








Tr 
sın(2?n+1l)ıx 





worin 








-) 


D(v)= /nJ® "du. 5 3 ’ - ; u e : i . Im 





das Gaußsche Fehlerintegral ıst. Die Reihe (24) konvergiert sehr schnell für große, (25) für 
kleine Werte f, so daß die rechnerische Ermittlung von %, («, t) keine Schwierigkeiten bereitet. 






Bei gegebener Schmelzkurve kann 9,(#,f) nur durch Integration von G (x, £:t,f(£) 
eefunden werden. 










a) It z=z, so wird mit f=f(z) und r=f(£) der Integrand nach (6) wie 


(2 )? 
R. je Ua) FO) und daher wie , = n 
2y alfi2)— fiS] Zynf (ke)(2—}8) 


Auswertung des Integrals daher sehr unbequem. 






unendlich, eine unmittelbare 







Zur Ermittlung von #(z, t) betrachten wir neben 2 == f(z) eine zweite Schmelzkurve 












/ r IE R 3: 
au H 2 > j ’ : ; i : > ' ; h ‚ N R - ; (ZU). 
Nach (14) ıst 
2 a | 
r\ | ı).|z || & _. -IdE | 
A Aa r (1 (1” 


und daher 





(27). 














(14) ist 


ems für 


24) 


(2D). 


(17) 


(25) für 


bereitet. 
e; u, FE) 


(6) wie 


ttelbare 


uUrVve 


(26). 
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Wird a gemäß 


22 
f (2) B (28) 


ee 


cewählt, so bleibt der Integerand in (27) auch für £>z endlich, weil die beiden (reenschen 
| 


: 5 - unendlich werden. Daher kann das Integral in (27) dureh 
aynf (e)(2—2) 


Funktionen wie 


> 


zeichnerische Integration ermittelt werden: %, [z, —] und 9, (z, fiz)) werden nach (24) bzw. (25), 
A 
) (2, — | nach (23) berechnet. 
Br 


( 


b) Ist ==z, so kann 9 (., #) nach (14) ermittelt werden, da der Integrand an der oberen 
Grenze gleich Null ist. Um aber die Ordinaten des Integrals zu verkleinern und daher die 
Genauigkeit zu erhöhen, empfiehlt es sich, die Temperatur 9 (, t) wieder für die Schmelz- 


’ 

4 
Eu 
7 


kurve =, zu betrachten, wobei a nach (28) zu bestimmen ist. Ahnlich wie für . 
dA 


wird mit t=f(z) 


ce 


Ka de, 7)+9 9x, - k| IG(a,&; fe), FO) — Gla, 5; , zu)|dE (9); 
4 


a” i u 


m « m. 


u) 
- 


5 


[6 


D 


=) ist nach (22), 9, (w,f) und 9, (.r, 
( ( 


3, rl nach (24) bzw. (25) zu berechnen. 
a D 


6. Näherungslösung für ? - f(z), wenn z hinreichend klein ist. Für die Schmelzkurve 


f 2. V<z: (30) 
af . je 2 . . . o . N : . P . . s i . +) 
ist nach (23) 
f =" h ] "T 0 = | b | . 0 3 ( = ) 
ı) E a) > d, D| 5 \e | > D| 5 e » E a8) (>) 
Nach (25) ıst 
hd 7 2 ze A) 
lim »,|2, z AL P(7)- 
> (0 0 
Wird daher «a, aus 
k JT Kr | r | I T Ao\ - (ku 
"; 4 Ko P\ .) e' -| | 4 ku D| .) |e PD .) SE BEE Ze . (82) 
bestimmt. so ist für (30) 
lım #9 lz. = | ®.. 
> N an / 
Hiermit ist gezeigt, daß f= s »„ für hinreichend kleine Werte z eine Näherungslösung für 


(ko 
die gesuchte Schmelzkurve t=f(z) ist, wenn «a, nach (32) bestimmt wird. 


7. Tafeln und Tabellen. Zur Ermittlung von # (z. f) nach (27) und 9 (.r, P) nach (29) werden 
zweekmäßige eine Reihe von Tafeln und Tabellen entworfen. 


Für die Berechnung von d, (z, f(z)) und #, (2, - nach (25) wurde eine Tafel des Gaußsehen 
a” 


Fehlerintegrals ® (x) mit 6 Dezimalstellen benutzt?) und zur Interpolation die ersten und 


zweiten Differenzen berechnet. 


2), Z.B. EmannelCzuber: Wahrscheinliehkeitsreehnung und ihre Anwendung auf Fehlerausgleichung, Statistik 
und Lebensversicherung, 1. Band, S. 455 bis 457. Unveränderter Nachdruck der vierten Auflage, B. G. Teubner, Leipzig 
und Berlin 1032. 
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Für die Ermittlung derselben Werte nach (24) wurden # 2, 3Baz und 5#: für 2001: 
0,02: ...0,49 ım Wiınkelmaß ın einer Tabelle zusammeneestellt. Ferner wurde mit Hilfe 
einer Rechenmaschine lg e”?' r lege” fürr- 1:2:... 100 mit acht Dezimalstellen berechnet. 
so daß lg e’’r, ]g e!”*r und lg e?>”*r für beliebiges r ohne große Mühe eefunden werden können. 


y- 


Für die Berechnung von P(#,y) genügen vier Dezimalstellen (S. 352). Für die Er- 


mittlung von 9, (a, y) und 9, |, =; nach (25) kann daher eine vierstellige Tafel des Gaußschen 
l 


K“ehlerinteerals®) benutzt werden. Für die Bereehnune nach (23) tut die bereits erwähnte 
Tafel für le e”’’ gute Dienste, ferner wurden sinz.r, sind r.«, sind. und die Lorarithmen 
dieser Werte für z—=0,05: 0,10: ... 050 mit vier Dezimalstellen berechnet. 


Für die Berechnung von Al, ",) und #[z. „) nach (22) und (23) wurde mit p 
«A 


(1 


5 der 


Ausdruck Jr P Pipe” für die in Betracht kommenden Werte 5 mit vier Dezimalstellen er- 
mittelt, wobei die Unterteilung der unabhäneizen Veränderlicehen so gewählt wurde. daß nur 
die beiden ersten Differenzen berücksiehtiet zu werden 
brauchen. 

Die erößte Mühe bereitet die zeiehnerische Er- 
mittlung der in (27) und (29) auftretenden Integrale. 
Da die Gestalt der Grenzkurve mit zenüzender Ge- 
nauigkeit ermittelt werden kann, wenn f(z) für 20,01: 
0,02: 2... 0,50 verbessert wird und «da die Integrale 
hinreichend genau bestimmt werden können, wenn der 
Integrand für £<z eleieh O0: 0,01: 0.02: ... bekannt 
ist, S0 Ist eine Tabelle erforderlich. alls der ), 1‘ mit 
vier Dezimalstellen für #0: 0,01: ... 0,49 und be- 
liebiges r entnommen werden kann. 

Für #»< 0,015 sind, soweit mir bekannt ist, die 
Werte 9, r bisher noch nicht berechnet worden. Für 
diesen Bereich ist daher die Tabelle neu zu entwerfen, 
wobeı r sehr stark unterteilt werden muß, wenn man 
mit den beiden ersten Differenzen auskommen will, 
da #,"® nach Abb. 4 sehr veränderlieh ist. Aber auch 
für r >015 führt die Berechnung der für die Tafel 
notwendigen Werte nach unseren Formeln (7) und (8) 

wobei man mit den drei ersten Gliedern der Reihen 
auskommt schneller zum Ziel als die Interpolation 
mit Hilfe der bereits vorhandenen Tafeln, weil diese 
entweder nach a’=are sınk'\, K:K’'’) oder 5?*®) 
angeordnet sind: es muß daher zunächst r durch diese 
(Größen ausgedrückt werden, und bei der doppelten 
Interpolation nach diesem Wert und nach » wären 
Differenzen höherer Ordnung zu berücksichtigen, was 
sehr unbequem ist. Aus diesem Grunde wurde die vier- 
stelliee Tafel für 9, ri" =0: 0,01: ... 0,49) völlig neu berechnet. Die Unterteilung der un- 
abhäneieen Veränderlichen r wurde so gewählt, daß nur die Berücksichtigung der beiden 














\bb. 4. 


ersten Differenzen erforderlich ıst. 


‚ıZ.B. E. Jahnke und F. Emde: Funktionentafeln mit Formeln und Kurven, S. 485. Zweiter unveränderter 
Naehdruck der ersten Auflage vom Jahre 1909, B. G. Teubner, Leipzig und Berlin 1928. 

ı, J). Hoüel, Reeneil de Formules et de Tables Numeriques, 3. ed., Paris 1901, S. 59 (abgedruckt in Jahnke 
Emde a. a.0. zweite neu bearbeitete Auflage, S. 117) gibt 9se für # = 0,05: 0,105... 0,50 und e=9", 180, ... S1 mit 
bis 4 Dezimalstellen. Die unter Leitung von G. Greenhill berechneten Tafeln von R. L. Hippisley (Smithsonian 

la 


Miseellaneons Colleetions, Volume 74 [Druckfehler 73 auf dem Titelblatt] Number }, Publication 2672: Smithsonian 


Un . r 
Mathematieal Formnlae and Tables of Elliptie Funetions, Washington, 1922) enthalten 430 und , 7 mit 10 Dezimal 


1 .) ng 9% er E ce 2 N 
stellen für ı . us — und HD. BRD... ; BOR TU ION Dem erößten Wert « 90 entspricht der 


Isı" 180 ISO 


kl ınste \\ er! ı 0.025. 
5) Die Reports of the Meeting of the British Assoeiation for the Advancement of Science, London (1912: Dundee, 
Birmingham und 1919: Bournemouth) enthalten dieselben Werte wie A) nit vier bis sechs Dezimalstellen für 


K: K'’ gleich 3V3: 5: 3 V2: 4 usw. Dem größten Wert K:kK'’=3V3 entspricht der kleinste Wert r z 0,015. 
Keiichi Hıvashi (Tafel der Besselschen, Theta-Kugel- und anderer Funktionen, S. 85 bis 10%, Berlin 1930, 
ringer und Fünfstelliee Funktionentafeln, S. 148 bis 154, Berlin 1930, Julius Springer) gibt die Thetanullwerte 930 
vol® 2.2.0500 mit 8 bzw. 5 Dezimalstellen, die Tafeln von Hautaro Nagaoka und Sadazo Sakurai 
rs of the Institut of Physieal and Chemical Research, Vol. TI, Pages 1-67, Table Nr. !, December 1922, 
Tokvo) enthalten die Logarithmen von #%,0 für k?= 0: 0,001: ..- 0,500. Dem größten Wert k? = 0,500 


- . 


kleinste Wert r z U,080, 








y 


.angew, 
ıd Mech. 


0.01: 
Hilfe 
schnet. 
önnen. 
lıe Er- 
schen 


vähnte 
tımen 


7 
„ der 
len er- 


5 nur 
‚erden 


he Er- 
orale, 
r (re- 
0.0]: 
errale 
ın der 
kannt 
» mit 
id be- 


t, die 

Für 
erfen, 
ı man 
will, 
"auch 
Tafel 
nd (8) 
eihen 
latıon 
diese 

k?®) 
diese 
velten 
vären 

was 
 vier- 
r Un- 
eiden 


ıderter 


hnke 
Ss] ımit 
Isonlan 
Isonian 


zıimal 


ht der 


undee, 
en Fir 


n 1930, 
te JIs0 
akurai 
r 1922, 


0,500 





a ir 





Band 15, Heft 6 
Dezember 1935 





Lachmann, Zum Problem des Erstarrens begrenzter Körper 3dl 





Wird für r <=0.035 nach (8) 
V).U = le a ge ir _e DE ER 9 DOREEN En 


vesetzt, so gilt für den Absolutwert AR des Fehlers 


| (2 — v)? (?2+v)? (3 — r)? (3 + v)? 
I? en > le ir B= P 17 I P ww Zr; tr 
ayrır 
N 1,52 32 2,92 3? 
r e tr 1 P tr 4 [Z tr n- P ht? un ) 
y+ Tr 


Der letzte Ausdruck wächst für r <1 beständie mit r und ergibt für r= 0,035 den 
Wert 0,00000016. Der begangene Fehler ist daher kleiner als + 0,16 Einheiten der sechsten 
Dezimalstelle. 


v2 (1-9) (14m? 120,01? 
P ir E r E tr P hı 
Die Logarithmen von und wurden gemäß la 
2yar dyrır Zyrr 2yrr 
N 0,012? 
le „lee #r für n=1,2,... 150 mit Hilfe einer Rechenmaschine ermittelt. 
2yaı 
0,012 


Hilfstabellen : Ig 


- und lege #r für die in Betracht kommenden Werte r. 
syrr i 


Für r > 0,035 erfolgt die Ermittlung von 9,» nach (7), wobei wieder nur die drei ersten 
(Glieder berücksichtigt werden: 


2cosaırv  2eosdrr  2cosdar 
( [} 
VUZ er + — r a (4). 
- P 7 ep! teY P-) te? 
Für den Absolutwert R des Restes gilt 
.) ») ») .) 
ER er 4 = EL. ee D - - ‘ 
R< et 7°r eSl ?r DI 04972 .0,055 eS1l 72. 0,035 ALLE 


Der Fehler bei der Ermittlung von 9, # nach (34) ist daher für r 0,035 kleiner als + 0,09 
Einheiten der sechsten Stelle. 


2ecosar 2cosdrarv 2cosaae, , ilfe der Reel 
») » N . 
. 2cosıı e . BRRR 
maschine gemäß |g le(2cosar)  rl&e"” usw. gefunden werden.  Hilfstabellen: 


e’ r 


e2>cosır, lei2cosdar und le!2cos5ärer. für vr =0,01: 0,02: ...0,49. 


= 


8. Die Lösung des Problems für kleine Werte %. Die Betrachtungen sollen für das 
gewählte Beispiel 
"),—=0S und k 01 
durchgeführt werden. 


a) Die Schmelzkurve. Die Werte 9,08 und k 0,1 werden ın (32) eingesetzt; 
dureh Interpolation wird a,= 1,6195 ermittelt. Die Näherungslösung für die Gleichung der 
Schmelzkurve lautet daher: 


/ — — (),38128 2?, 
17 Pre 


die Temperatur auf dieser Kurve nach (23): 


2 \ at 
az, ) 0,20679 + 0,7983219, (2, 1... en. (BB). 
17 Pi Au) 
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Wird #,lz. “ | nach (7) berechnet und in (35) eingesetzt, so ergibt sich, daß bis z=021 die 
GFn 


Temperatur 9 \z.  ,)° 0,0000 ist, wobei der Fehler kleiner als eine Einheit der fünften Dezimal- 
2 

stelle ist. Für 2> 0.21 soll die Schmelzkurve soweit verbessert werden, daß die Temperatur 

für jeden ihrer Punkte um weniger als eine Einheit der fünften Dezimalstelle von dem wahren 

Wert #9, —=0,S abweicht. 


Da nach (28) @® dureh f’iz2) bestimmt ıst, ıst es zweckmäßig, nicht die Schmelzkurve 
unmittelbar, sondern ihre erste Ableitung f’(z) zu verbessern: diese Kurve wird gezeichnet, 
während die Werte ftz) nur berechnet werden. 

h Naar: % a 2.0.21 
Bisz 0,21 ist f’ (z)eine Gerade dureh den Ursprung, für z=0,21 ist f’ (2) 0,16014 
Go 
(Punkt A in Abb. 5). Für z=0,22 ergibt die Verlängerung der Geraden den Wert f’(z) 
>) ,() () 22? 
0.1676 (Punkt O): nach (23) ıst 30,22: — 0.799965. 
A 


E” 


Auszehend von 4 wird nun die Kurve f’(z) etwas nach unten abgeboren, z. B. so, daß 
für 2=0,22 die Ordinate fiz)=#’=0,167 ist (Punkt 1)%. Mit der mittleren Ordinate 0,1640 
0.2] 


zwischen z=0,21 und 0,22 und #(0,21) —, 0.016814 ıst für 20,22 der Wert #, 0, 016814 
A, 
-0,01-0,1640 = 0015457. Die Temperatur 9 (0,22: £)=#, ist nach (27) zu berechnen, wobei 
die in Absehnitt 7 erwähnten Hilfsmittel benutzt werden und das Integral zeiechnerisch zu 
, 22 2” BZ 0,22 er ” 
ermitteln ist: nach (25) ıst f (2) ‚„ und daher +, > 0,167 = 0,018370. Es wird 
4 & a 


I, OST. 

Die Temperaturen 9, und 9, werden als Ordinaten zu den Abszissen #, und #, auf- 
eetraeen (Abb. 6). Bei dem gerinzen Unterschied der Werte darf aneenommen werden, daß 
sich # linear mit f’ (2) Ändert. Wır verbinden daher die Punkte O0 und 1 in Abb. 6 durch 
eine Gerade: für den Schnittpunkt 2 mit 9=0,5 ist f’(z2)=1t, = 0,1675. Dieser Wert wird 
in Abb. 6 einzeetraren (Punkt 2). und es werden die Punkte A und 2 durch eine Kurve ver- 
bunden. welehe die Kurven 10 und Al angenähert verhältnisgleich teilt: die mittlere Ordinate 
der Kurve 12 zwischen 2=021 und z=0,22 ist 0,1638, so daß 


f10,22) = #, = O,016814 + 0,01 0,1638 — 0,018452 


ist. Für 2022, fi)=t, und f (z)=t#, ist die Temperatur 9, nach (27) zu ermitteln und 
zu untersuchen, ob dieser Wert tatsächlich um weniger als eine Einheit der fünften Dezimal- 
stelle von dem wahren Wert 9,085 abweicht. Da dies der Fall ist, darf Punkt 2 in Abb. 
als wahrer Kurvenpunkt angesehen werden. 

Um Rechenfehler zu vermeiden, Ist es zweckmäßig, 9, (z, f(z)), in Abhängigkeit von t = f(z) 


‚|e »,\z2.,)) in Abhängigkeit von f’(z) „ aufzutragen. Die den 
J Ei 


_ _ (7 


Aa 2° 22 


une ei Ku a»| 


Punkten O0, I und 2 entsprechenden Werte müssen angenähert auf einer Geraden liegen. 


PO) En EEE 
A 530.8) 0,167 f (2) 


Abb. 6. Abb. 7. 


und I zu wenie voneinander unterscheiden, ist Abb. 5 nieht maßstäblich gezeichnet. 
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Die Kurve A2 (Abb. 5), wird nun freihändig bis zum Schnitt mit z= 0,23 verlängert 
und für diesen Wert das Verfahren von neuem zur Anwendung gebracht. Ebenso ıst für 
»— 024: 025: ...050 die Verbesserung der Schmelzkurve durchzuführen. Da sich die 


Neieune der in Abb. 6 gezeichneten Kurve für aufeinanderfolgende Werte z nur wenig ändert, 
so kann man den zweiten Wert f’(z) stets so wählen, daß für ihn 9 etwas kleiner als #, 
wird. wenn 9 für den ersten Wert f’(z) größer als 9, war und umgekehrt. Hierdurch 
wird erreicht, daß 9 für den nach Abb. 6 interpolierten Wert f’(z) genügend genau mit %, 
übereinstimmt. 

Um Abrundunesfehler zu vermeiden, ist es angebracht, sämtliche Rechnungen mit sechs 
Dezimalstellen unter Benutzung einer siebenstelligen Logarithmentafel durchzuführen. Bei der 
Ermittlung des Interrals in (27) hat man den Integrand für genügend viele Zwischenpunkte 
zu bereehnen und die Zeiehnune in möglichst großem Maßstab anzufertigen, um das Ergebnis 
mit genüzender Genauigkeit zu erhalten. 

Es braueht wohl nieht betont zu werden, daß das Verfahren sehr langwierige ist und 
mit peinlicher Sorgfalt durchgeführt werden muß, wenn die endgültige Kurve f" (@) (Abb. 7) 
einen glatten Verlauf zeigen soll. Die ermittelten Werte f(z) und f’(z) wurden in der nach- 
stehenden Tabelle zusammengestellt. 




















(z f (2) ? !’ (2) / f / 

0) 0) ( 

0.01 0.0076 V.VOOO4 0.2] 0.1601 0.0681 0.41 0.1519 0.0376 
0.02 0.0153 0.00015 0.22 V.1678 0.01845 0.4? 0.1472 005526 
0.03 0,0229 0.000234 0.23 0.1750 0.027016 0.43 0,1445 0,0671 
0.04 0.0305 V.OOO6GI 0.24 0.1817 0.021095 0.44 0.1427 005815 
0.05 0,0581 0,00095 0.25 0,1884 0.025380 0.45 0,1437 VUDODS 
0.06 VOADS 0.001537 0,26 0.1940 VH2571 0,46 0.1469 0.061102 
().07 0,0534 000187 0.27 0.1983 0.02767 0.47 0.1513 VOH251 
VOUS VOHTO 00244 0.28 0.2020 0.029657 0.48 0.1564 0.0640 
0.09 ).OHS6 0.003009 0.29 0.2035 VO3170 0,40 0.1610 HIERTER | 
0.10 0.0763 00381 0.30 0.2043 0.03374 0.50 0,1658 0.067127 
0.11 0.0839 0.0046] 0.31 V.>VUP>8 VOBDITS 

12 0.0915 V.OOS4S 0,32 0.2006 V.OBTS0 

0,13 0.0091 000644 0.33 0.1970 VOBITS 

0.14 0.1068 0,00747 0.34 0,1920 004175 

0.15 0,1144 VOOSDS 0.39 0.1872 0.043683 

0.16 0.1220 000976 0.36 0,1807 0.0447 

0,17 0.1296 001102 0.37 0.1751 0.04724 

0.18 0.1373 0.0123 0.38 0,1690 0.04896 

0.19 0.1449 0.013376 0.39 0.1624 0.05062 

0.20 0.1525 0.01525 0,40 0.1565 0.05222 


Zu beachten ist, daß f’(z) für 2 0,5 verschieden von Null ist, so daß die Grenzkurve an 
dieser Stelle eine unstetige erste Ableitung besitzt, während f’ (z) für k=0 und z2=0,5 gleich 
Null und daher die Ableitung stetig ist (Abb. 7). 

b) Temperaturverteilung. Da die in die endgültige Abb. 9 einzutragenden Werte 
der Temperaturen nur auf drei Dezimalstellen genau zu sein brauchen, genügt es, sämtliche 
keehnungen mit vier Dezimalstellen durchzuführen. 


Wir bereehnen zunächst 9 («, t) für «= = 0.05: 0,10: ...0,50 und die 20,05; 0,10: ...0,50 
entsprechenden Werte t=f(z). Für 20,20 ist nach 5.352 =, und daher kann (a, ®) 
Oo 


f 57 


„. 
) (* | nach (22) berechnet werden. Für 2>0,20 ist es zweckmäßig, Pd, # nach (29) zu 
Ko 


ermitteln. 
Die gefundenen Werte (x, f) werden in Abhängigkeit von f aufgetragen, wobeı .r als 
Parameter aufgefaßt wird (Abb. 8). Da Pla, t) bei konstantem f an der Stelle, wo da, M 
"0,8 ist, eine unstetige erste Ableitung besitzt, wurden in Abhängigkeit von f diejenigen 
Werte .z aufgetragen, für welche 9 (x, =», ist. Zu diesem Zwecke braucht man nur die 
Werte z der Tabelle in Abhängigkeit von f(2)=t aufzutragen (Abb. 8). 


\ 
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Für 0,05: 0.10: ...050 und t=0,005: 0.01: 0,02...0.06: 0.06727 werden dieser 
Abbildung die Werte 9 (#,#) und für dieselben f diejenigen Werte x entnommen, für welche 


"(ae,t\-=9. ist. Es können dann die ersten Kurven in Abb. 9 gezeichnet werden. Für z2=05 


ist nach S. 353 
se . 5 ee rn ah 


für diesen Wert seı 
Bia,f)= ala) 


die Temperaturverteilung, wenn der ganze Körper erstarrt ist. Diese Kurve ist in Abb. 9 


dureh Striehelung hervorgehoben. Für f> #* haben wir mit 


EB) A er 


/ 
das folgende Problem: zur Zeit n=0V ist die Temperaturverteilung 9 (#,0)= g(x) gegeben. 
Für y>0 ist JO, )= #9 (l,y)=0d, während die Differentialgleichung in den von uns gewählten 

",) Vu DE 
Kınheiten (2) a lautet. Es ıst 
/ 
e 4 
d(a,n)=2,) e Fr sinnaelgö)sinnneädE. 


[7 | 0 


Mit y(1 ° H)=yg(£ und 
sinn ze für gerades n 


sınnall -&) ws 
sinnzc für ungerades n 


1/ 


” 


"sinn +l)ar\ygd)sin@n+l)a&d: 


ı 


(n+]1 


r. . {1 » . 6 R BEN. ;, n F zu. 
"sinzc+t4J,e"""sindrc-H4J er """sindnz 


pe 
In Sg sinn +i)nidE 
ı) 


(39) 


Es wurden zunächst J,,./, und ./, dureh zeichnerische Integration ermittelt; J, und 
alle weiteren Integrale in (38) sind kleiner als + 0,0001 und brauchen daher nicht berück - 
Nach (35) können dann die Temperaturen 9 (r, 1) für #===0,05: 0,10; 


siehtiet zu werden. 
" .... d.h. t=0,07: 0,08: ... berechnet und ın Abb. 9 ein- 


050 und 7 — (). 07 t: O.OS 
eetraeen werden. 
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Wegen der geringen Schmelzwärme (k#=-0,1) ist das Bild der Temperaturverteilung 
anzenähert dasselbe wie bei der bekannten Lösung für k®=10 (Abb. 10). Nur sind die 
Temperaturen zu derselben Zeit # etwas höher und für £<#* ist die Neigung der Kurve # (.r, #) 
an der Stelle, wo 9 (r,#) = #9, ist, unstetig, was in Abb. 9 besonders für #=0,05 und 0,06 zu 
erkennen ist. 





= 


< 
. 
m 


—n 
ID 
m 


7 





»lt,t) für k=-0 weich Mt) 








Abb..10. Abb. 11. 


9. Die Lösung des Problems für große Werte %. a) Die Schmelzkurve. Während 
die Schmelzkurve für kleine Werte k bis 2-05 verbessert werden muß (Abschnitt 8), wird 
diese Verbesserung für große Werte % dadurch wesentlich erleichtert, daß es neben der 


.) & ; 
Geraden “,, der sich f’ (z) für kleine Werte z nähert (S. 349), eine zweite Gerade — , gibt, 
Ay A, 


von der sieh f’ (z) für große Werte 2 immer weniger unterscheidet (Abb. 11). 
Es gilt nämlich folgender Satz: Läßt sich für einen festen Wert z eine Konstante «a, 
so bestimmen, daß 





, \ dı?2 ı2 2 ö ) Yr? M) 
u Va le en kV fa) St ee) 
e-"?1Kı) r a, b| )e te a1? Hr > ey Tl i2ı) E a, b| y et e ad,* 
7 y | 7 y Ye 
i 12 
. . i = (2ı? s?) z £ 
B= l k \ a 12 (f(z,) (= P u Eu sın 7 & d E ( 0) 
0 
7 072 } 
f 7 b 3 
Lik\l\e rt) ID) —e *°r sindr& ds 
0 


kleiner als & ist, wobei + e die größte zugelassene Abweichung der Temperaturen von den 
wahren Werten ist (S. 352), so darf für 22, die Ableitung f’(z) der Schmelzkurve gleich 
Fa a,? 
Pa . k} T ck, 1 ( s pe Ne io ( (m 

„ gesetzt werden, wobei a, aus \, a P|\5 )e », zu bestimmen ist. Ferner ist A (a. fi2)) 
," r Z 


23,) 


2} 


/ 


für 22, und zz gleich »,, während für «2 die Temperatur 9 (a, f(z)) = U 


vesetzt werden darf. 
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Dieser Satz läßt sich folgendermaßen beweisen: Mit f’(z) —, und 
Pe 


ky: a,z\ 


T 
”(a,f(z)) „a,d| ar \e für sSz 
I(a,f(z)) Y Be 5% (41) 
( ‚ k] T A, Jr 
Jr, / (z)) 5a, D| 7‘ e*+ für x 


Qu 








ıst nach (22), (24) und (29): 


ar? 12 2° 


2 7 7 = 
Aca#,f(z)) le er. & D|)e! e “* \sin.o 


MN 9... 
(4,2 u712>2 
| . 7 } . - 


a7? |sin Ink 


sınrzsinm£d:s 


sın3rzsindn£ds 














"ür z2=z, sınd die Absolutwerte der Ausdrücke in (42) kleiner oder höchstens gleich den 
entsprechenden Ausdrücken in (40), es ist daher A (w,f(z,))|<e und nach (41) angenähert 


D) 


kya 2 "ya a\ 
+ (a, fl2,)) 3 V a,b| > — \e für 2=3,, 4 a, D| 5s\e* für & (45). 


RR 


Qu 


Für #=z, muß 9 (2,.f(z,))=%, Sein, es ist daher 
kya (a) & 
a a,P\Z)e' 2 Eee 


m 


und la, 
P\ = — 
”(a,f(z2,)) 3 —- 9, für zS32,, RER *: 3 1 as MB 
(k, 
»\ z 
Für einen beliebigen Wert z, der größer als z, und kleiner oder gleich 0,5 ist, werde 
). 
f (2) Vi EEE EEE 
a, 
vewählt, d.h. die Gerade O0 1 ın Abb. 11 über A hinaus verlängert. In (42) kann dann als 
obere Integrationsgrenze z, gesetzt werden, da für 2, = £=z die Differenz fe) -f)= , 
A 


ist und daher für dieses Intervall die Integranden sämtlicher Integrale verschwinden. Wird 
nach Abb. 11 








f@)-fe)= —— ee A 


eesetzt, so Ist mit (42) 


x & dir T22,° 
Ia,f@)) Ä 0-3 '(e 12 / (2) h j = TR "u | 


7 4 
’ ‘ 9 n22,2 
| ee aan h: | 7 a, 2 a 
—ur "site PEN - a» je e ”,’ sind 
+») .7 \ Fe 2 / 
. 12 ] (48). 
-+ke-"s/) je-PV/W-IC) —e *ı sınzaesınnazSdE 
u 
4 IE 2 ” 
ısındrzaesindaSds 
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Nach (40) ist !A (@,f(z)) <e- "#7 <e, daher darf auch für z>z, nach (41) mit einer Ab- 
weiehung, die kleiner als + e "?+/e ıst, 
































kya (a,.\ kya a, a 
Am I (a,t) > 9 a,0\5 : )e ı für 2S8#, > a,b 5 )e I! für zz (49) 
oesetzt werden. Mit (44) ıst für z>z, wie für z=z, 
( a, d" ( ... ( d .. - 
(x2.D > | ae \9, für 2S2, ", für Bed (U): 
hiermit ist der Satz S. 355 bewiesen. 
Das folgende Beispiel wurde gewählt: Kupfer mit dem spezifischen Gewicht y=-SW0 kg/m’. 
Temperaturen: Anfangstemperatur 1330°, Schmelzpunkt 1080°, Temperatur der Umgebung S0" 
daher nach (2) 
40) , _ 1080-80 
). ), IS. 
1330 — S0 
Spezifische Wärme: e= 0,12 keal/kg’, Wärmeleitfähigkeit 4 312 keal/m st", daher Temperatur- 
RE /. 312 R 
leitfähigkeit a = u 0,292 m/st. Schmelzwärme 4 42 keal/kg. Abstand der Be- 
| 0,12. 8900 | 
r En a a; N 
erenzungsebenen 7= 0,05 m, daher X = T==0,05.- ın m: 7 000556 Fin st. Schließlich 
( 
ist die Konstante 
h den q {> . 
ähert y [ 0.12 RE . 
(43). Nach (32) ist a, 0,067252 für 9,08 und k® 350 und daher 
I 
f, (2) — 442.182 ee 
In 
(44) Er 
wenn z hinreichend klein ist. Für genügend große Werte z ıst nach (44) «a, == 0067585 
und daher 
.) » 
(45). f, (2) = — = 4831,86 2 (52) 
A, 
werde Die Werte f, (z) und f, (z) unterscheiden sich so wenig voneinander, daß die zeiehnerische 
Integration von f’(z) zu ungenau werden würde. Aus diesem Grunde wurden in Abb. 12 die 
(46) 
dt 
ın als 
A 
Wird 08 
(+7) 
035) [754.735 
S 
020 8 
Ss 
>\ 
0) | 
5 /z) / [z) 
(48). 00) 
| 
005, 
0 un 000 0020 0 z 0m Da Da un 0 27 . zu. 
Abb. 12. 
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Abweichungen f, (2)  f’(z) aufgetragen und verbessert. Durch zeichnerische Integration wird 
f,(z) fiz) gefunden: zieht man das Ergebnis von dem berechneten Wert f,(z) ab, so erhält 
man f(z). Sonst ist das Verfahren wie für kleine Werte %, Abschnitt S, durchzuführen. 


Die Zeichnung läßt erkennen, daß bereits für z=0,06 die Ableitung f’(z) der Schmelz- 
kurve gleich fj'(z) gesetzt werden darf. Für 2=0,06 ist f(z) — f,(z) = 0,00193, so daß für 
0.06 


f(2 f, (2) + 0,0193 = 218,95 2° +0,00193 und F(e)=f\ (2)=437,856z2 . . (53) 


ist. Die wahre Sehmelzkurve f(z) verläuft zwischen den Näherungskurven f,(z) und f,(z). 


b), Temperaturverteilung. Für 20,06 werden die Temperaturen wie für kleine 
Werte k, 8. 353, ermittelt. Für 2=0,06 ist (40) mit 2=-0,6001 erfüllt, die Temperaturen 


rn * - but V.. a “)»)L) 4 
können daher für 0,06. 2-05 nach (50) berechnet werden. Für w = ,„ 0,0338 darf 


-) 


P(w)> .. m gesetzt werden, es ist daher mit (50) 
JT 


Jr, - ». Mir 2 B. y E Eur  . 2... (24) 


und 0,06 °.2°.05, Abb. 13; 54,755 ıst wie auf S. 354 die Ordinate der Schmelzkurve 
für 2.05 und g@r) die Temperaturvertellung für denselben Wert z. 


Mit / “ y und gar) 1,68 nach (4) wird mit 89) 


,.0 
, 


[1.6 8sin (2 .Ny2d} 


und «daher nach (38) 


ir .r 005: 0.10: ...050 und 7 — 94.05 u": 4.0 n” 
und 55 97" wurden die Temperaturen ermittelt und in 
Abb. 13 eingetragen. Das Bild der Temperaturverteilung 
entsprieht annähernd Abb. 3 (#350 und 9,1 an 
Stelle von 0,5), ist aber sehr verschieden von Abb. 9 
(01 und 9, =08S). Daß die Schmelzwärme das 
Abklingen der Temperatur stark verzögern kann, läßt 
auch Abb. 14 erkennen, in welcher die Temperaturen 
Mitte Stab für 9,05 und #=0: 0,1 und 350 in Ab- 
hängiekeit von f aufgetragen wurden. 





Wegen des geringen Unterschiedes von f,(z) und 
f,(2) (8.557) ıst > f,(z2) In dem untersuchten Beispiel 
"„, 05 und k=350 bis z=05 eine brauchbare 
Näherune für die Schmelzkurve. Die nach (23) er- 





mittelten Temperaturen 9 (2,  ,) liegen zwischen 0,8 und 
A, 

v7. Das Bild der Temperaturverteilung entspricht 

Abb. 13. nur sınd die Zeiten etwas zu eroß,. es Ist z. B. 

2 »».2 für z=05 an Stelle des wahren Wertes 


4.190 90 





vr 


f.angew. 
nd Mech. 


ın wird 
‚ erhält 
n. 


-hmelz- 
laß für 

(29) 
I f.(2). 
kleine 
raturen 


S «darf 


(54) 


zkurve 


(DD) 


(36). 
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Beitrag zum Problem der Spannungsverteilung in Gurtplatten. 
Von Erich Reissner in Berlin. 


ie Spannungsverteilung in den Gurtplatten von T- und Kastenträgern ist bisher nur für 

Belastungs- und Stützungsbedingungen behandelt worden, wie sie sich bei Spannungs- 

zuständen ergeben, die entweder in der Längs- oder in der Querriehtung periodisch sind') °). 

Der auch technisch nicht unwichtige auf nicht periodische Spannungszustände führende 

Fall eines sehr langen nur in einem endlichen Bereich belasteten Trägers ist jedoch eben- 

falls der Berechnung zugänglich, indem man den ebenen Spannungszustand in der Gurtplatte 
dureh Fouriersche Integrale darstellt. 

jr’ 


F 
LEICE) Pat ] de 177 


Abh. 1. 











Die hier mit diesem Hilfsmittel durchgeführte Berechnung eines unendlich langen, auf 
zwei Stützen ruhenden, durch beliebige Kräfte, insbesondere dureh eine Einzelkraft in der 
Trägermitte belasteten Trägers (Abb. I) mit unendlich breitem Gurt wird praktisch besonders 
einfach dadureh, daß sich die wiehtigsten Größen durch tabellarisch vorliegende Funktionen 
angeben lassen. 


1. Allgemeine Aussagen über den ebenen Spannungszustand in der Gurtplatte (Abb. 2). 
Die Spannungen eines ebenen Spannungszustandes lassen sich bei Abwesenheit von Volumen- 
kräften bekanntlich wie folgt als Ableitungen der Airyschen Spannungsfunktion schreiben: 


>) 


























N2 7! 12 7’ 
ie a K a F MA fläche 
& N 42 3 )] ur we ' mn a 
07 UM t h 
a | in... n %p 
Po F ' r 7 nn 1 
T 
7, 
Mit Hilfe des Zusammenhanges zwischen Spannungen und 2 fi 
Y .. _ r® 
Formänderungen 
! 
a7 | | | | Ba Er 
pr N, 6, 
“= Er | | 





( " | | | ) | | 
E, 10, Öx ur (-) 
v7”oy KE\Y m? u I 3 
ou, or | | " 
y= 1” -T 
: oe 7] Or (1 Dr en 








erkennt man, daß F der folgenden partiellen Differentialgleichung Asıeı) 
eenüren muß: 
Abb.2. Der Sprung in der oy-Fläche 
0’ F 0” F 0’ F rihrt bekanntlich von der Nicht 
> 9 v0) (9). berücksiehtigung der Querkontrak 
tion im Steg ber. 


IAF 


Or?  "or°0 y“ 7 


Die Lösungen der Gl]. (3) haben außerdem gewissen, von Fall zu Fall verschiedenen, 
Bedingungen an den Gurträndern zu genügen. Nur die Bedingungen an der Anschlußstelle 
Stee Gurt können ein für allemal formuliert werden. Für einen T-Träger folgt erstens aus 
Symmetriegründen, daß die Verschiebung quer zur Stegachse verschwindet. Also 


Be ne anne a re te ee A 


I), T’h.v. Käarmän: A. Föppl Festschrift 1924, S. 114 bis 127. 

?) G. Schnadel: Jahrb. d. Schiffbaut.-Ges. 1926, S. 207 bis 29. WRH 1928, S. 92 bis 102. 
*) W, Metzer: Luftfahrtforschg. 1929, S. 1 bis 20. 

°ı, MH. Reissner: ZAMM 1934, S. 312 bis 318. 


5) E. Reissner: „Stahlbau“ 1934, S. 206 bis 208. Norske Vidensk. Selsk. Forh. 1055, S. 15 bis 16. 
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Kine zweite Bedinzung liefert die Aussage über die Dehnung längs der Trägerachse. Wenn 
vorausgesetzt wird, daß im Steg eine lineare Normalspannungsverteilung herrscht, muß sein 











Q Mur) 5 
Exrir,V) = EW in aan Ar: sure Ver (>). 



































Dabei bedeutet W dr) das obere Widerstandsmoment des Trägers einer mittragenden Breite 
> /kr) und Mr) das angreifende Biegungesmoment. Zwischen W und / besteht (bei Vernach- 
lässieune des Trärheitsmoments des Gurtes um seine eigene Schwerachse) die folgende 
Beziehung: 
hı 
N rn kArBäR . Ah) Ä 6). 


’ 


Als Definitionseleiehune für 4 hat ıman (s. Abb. 2) 





I) 
Gr dy 
: er | v. | | A 
Mir) 
Wa) 
Aus (6) und (7) erhält man 
13 Ic v 
a dh MW \ va dy) So 


Dies in Gl. (5) eingesetzt liefert als zweite Randbedingung 


h) 
Sol th) 


hd \ EA day h°« 


„a. (S). 


u 


Für die tragende Breite ergibt sieh aus (6) und (7) die folgende Formel: 
I 
Wd\o,.dy 


ı 


A 1, 


bj WM (.r) Ra) lı () Öx dıy 


[IE 


(*)) 


2. Der Träger auf zwei Stützen mit unendlich weit auskragenden Enden. Die Spannungs- 
funktion F werde unter Voraussetzung von Symmetrie um = 0, mit einem partikularen 
Integral cos a. fiy, a) von (3), in der folgenden Form angesetzt: 


Fi 


Fe \cosar- fma)da. 


/ufolge (3) muß fir a) der totalen Differentialgleichung 
of Zar fi | fiv) -f) 
venüzen, und man erhält 


FF \cosawt[|Ata) + Bla) ay]e=*,+|Cta)+Dia)ayje'"Hda . . . . (10), 





worin Ita), Dia), Cta), Dia) beliebige, durch die Randbedingungen zu bestimmende Funk- 
tionen sınd. 


Nimmt man weiterhin an, daß der Gurt unendlich breit ist (b= x), so muß wegen des 


Verschwindens der Spannungen für y=x sein 


(a) Dia) Bi Deu em  Marsrwmränneiat A FE 








Wenn 


ıuß sein 


(>). 





n Breite 
'ernach- 
foleende 


(6). 


(IS). 


()) 


INUNgS- 
ularen 


10). 


Funk- 


nn des 


(il). 
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Mit Hilfe der Bedingung (4) kann Dia) dureh Aa) ausgedrückt werden. Die Spannungs- 


funktion nımmt dann die folgende Gestalt an: 





v4 
m | m | 
l \ Italcos a.r \ ıyle (de (12). 
. -) . 
| Zzm im 
ı 
und die Ausdrücke für die Spannungen werden 
I 
3m | mn I 
Ör \ a: Ala) cos a. ) ) (9 „da 
i _ 4 iM 
() 
74 
De im | in | 
[2 \ u la) cos a. (U) dla (1,3) 
7 .) -) . 
} Im Im 
ı 
SI 
ai . [IE - 1 
T \ a? Ala) sin a. | a are la 
| Im 
0 
Die Bedinzune (8) des stetieen Übergeangees von Ster auf Gurt wird 
I 
\ Ä 3m’+-2m -|1, 80 ’ 6 v 
a .tıta) |) a = m cos a. da Mila) , .‚.ISal. 
| er it“ dd hd 


Um aus dieser Gleiehune die Funktion Afa) zu bestimmen. ist es nötie, Mr) wie folet 


dureh ein Fourierinterral auszudrücken: 


FT FT. 
») 2 F 


Mir — osiax Miz\cosizaldz) da. 
[(.r) \e (du IN. I 1) 


0 (u 


Damit ergibt sich, wenn man noch zwecks Abkürzung setzt: 


241 m? I6 9 m’ 
> r 3. Hl... r wer in 
hr dam’ + 2m 1) h dS m? + 2m 1) 

JS 

m, i 

Aa) \ Miz) cos (za) dz 

aa+m;,). 

0 


Die vorstehenden Formeln mören aneewandt werden 


Träerers mit Einzelkraftbelastune. Man hat dann 


für 02 I M(z) Pıl 2), 
für I<sSz : Mi2=%b. 
Daraus ergibt sich 
4 ] 
| eos al 
\ Mız) cos (za) dz P\ (l—2) cos(za) da—=P FR 
N N 


Also wird nach (16) 


Pım,.(il -eosal 


Ala) | 
aa + m,) 


(Ih). 


- 
(1) 


(16). 


auf den Fall des 


(Ita). 


(ba). 


24 
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Für die Bereehnune der tragenden Breite nach (9) braueht man 


4 SF # 


r | 
( IF ecosazrıll cos al) z 
6 OS: » . der er (1). 
\ ‚dy r = 2 \ a Aa) co ' / nn. \ (at m) f 


Ks ist für die numerische Auswertune zweekmäßie, diesen Ausdruck folzendermaßen weiter 
zu zerleren: 


SF 4 P4 


of „, Ei ecosazxrıll cosal) cos all eosal cosarıll COS al) 


| 1.) 
 \ Fr da \ a da \ ee ( a) 


. l 
( ) ) 


Die ersten beiden Inteerale in der Klammer lassen sich dureh elementare Funktionen aus- 
drücken, das dritte mit Hilfe der Funktionen 


A LK 
2. sın u en cos H 
Stdn) \ dl U: Url) \ du . 
u Nu 
u x 


welche „Inteeralsinus* bzw. „Inteeraleosinus* eenannt werden und für welche Tabellen vor- 
handen sınd®), ”) dureh foleende Umformune 


I I 
“costa 2) cos (a, 2 u) 
\ da=\ du. 
\ata, " 
ı (2 
Setzt man noch zweeks Abkürzune 
. : JT an: 
cos x Ü t(.r) — sın r > SI (2)j == Dr). 
so erhält man mit = .r/! und mit 
nl Ib m? / Ri) 
im - »,: oO . a 
hd3m?’--2m | 1 - Jı dl 
wobei die Größe o als „charakteristische Trägerzahl* bezeichnet werden kann. für || | 


OF 3Pl {n n | . 
n,(1-E)+In = Pin,‘ Pin, (1+S5)+ Pin, \1- El; (17a). 
R 7 | 2 2 N J, ", | ) 1 | e2 1 .) | | | r (lu } 


- — 


Damit wird nach Formel (9) für © e< ] 


& | | 
1—£+- ur „» Plin dr 5 [Pin d+H)+ Pin, (li :))]\ 


tt m a. 

a FE 9 y Ä t I (da). 
.) Ill" T ps n ° “ 

‚+ Pin, d)) 5; [Pin 1 +D)+ Pin, (li Jar 


1—£ 2 


"ür «die tragende Breite im gefährlichen Querschnitt (re =0) ergibt sich daraus bei Beachtung 
der Tatsachen, daß für £° 1°) 


Cild)=EC+Iins; SIE) TE; Ba ;; 
») 
2 | U+Inn, — Pin,))] 
iQ m’ In, ’ 19 
U) - - i ae TE 
Sn t+-2m—| C+Inn, Pin, | 


‚Jahnke-Emde: Funktionentafeln, 2. Aufl., S. 78 u. f. 


') Man erkennt übrigens, daß für alle praktisch vorkommenden M (x) die Lösung mit diesen und den elementaren 
Funktionen gegeben werden kann. 
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Ind wenn ın 10/3 


MIR in er u 
| [0,577 + In (4,55 0) — P (4,55 o)] 
.2 0,445 \ 14 
TH In D(455 0) ale 

Für o- 1,5 ist mit einem Fehler kleiner als 1" 

0,14 
In(S.14 0) 
/(0) OHR 0 . (91 
21 In (5,11 0) Ben 


Aus (19) ersieht man mit Hilfe der Reihenentwiekluneen für P (.r)®). daß 


= 7.0) | 
lür 0>V 9j Funimmt wie In 
u y #) 
= (0) | 
für 0x \ abnımmt wıe 
. 3] Ino 


Aber auch im Bereich der praktisch vorkommenden Werte von o ergibt sich eine nicht 


’ ‚0 i 
unbeträchtliche Veränderlichkeit von >] mit o (Abb. 3). 











/ 
\bb. 3. Die reduzierte trarende Breite Im 
3] 


005 | | | | refährlichen Querschnitt Für 
I. den Träger mit unendlich-weitüberkragen 

: den Enden bei Kinzellast, 
| Il. den an den Stützen aufhörenden Träger mit 


L u ” | | \ | | starren Endquersehnitten bei Einzellast *S), 
Mn a 9 ni L6 u &0 2 111. denselben Träger wie II bei eosinusförmiger 
)=—+-:7 Momentenfläche. 
Pa 


Ks möge noch der Spannungsverlauf an der Anschlußstelle Steg-Gurt (y 0) betrachtet 
werden. Nach Gl. (13) und (16) wird 


ST. 
Sm+llcosaxil eosaN 
Gr (r, VÖ) m. J’ : \ da 
zm ala mt,) 
u 
I 
.) N . 
in 3m | | | 
- r ) | leosawıl COS ade. 
ni, m Na «a rm, 
( 
+.) = >) 
Sm | & | Pi 
0x (x, 0) = - z——— iin >; Pin, EI) +—[P (nn, | I+E) + Pin, |1- El] vr - (13a). 
Sımol 1 &* p | | WW, 


Dabei bedeutet W, = h? d/6 das Widerstandsmoment des Trägers ohne Gurtplatte. 


L 
: ’ /(N) ] | N , ] , . 
°) In diesem Falle ist ‚wo \ — ‚ wie zuerst v. Karmän a.a.0. 
21 to (k 1) l 12 ur 1,7 


son +l)?ljir < (?nH I) 
| ( ) 


festrestellt hat. 


Brit 
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Für «ie Schubspannung längs y=0 ergibt sich 
4 v4 
sınas(l- cosal) In, /1 | 
Tr I2:0) m, l’\ da P\| Il eosalsınardea. 
ala m,) nt, A arm, 
ı Mm 
l.ın Teil des Interrals kann 
von der Form 


elementar ausinteeriert werden. 


der Rest führt auf Funktionen 
’ JI er : 
Pix) Silke) -z )eosc + like) sin... 


Nach einigen Umformungen erhält man 


lür © S I rix 


« 
1 
— 


II 


r 
Yin, d9+5 |’, +5) - Pin, (1 3))] Ä 
für 1 Pin, d+- 


zoW, 


. (13b). 
| N 
g Il, 7 


oW, 
(Abb. U. 


IY’(ın, +) + ln, \ 
Die Formeln (13a) und (13b) wurden für einige Werte von n 





zahlenmäßie auseewertet 


—,__ 








03 
Sa 
SIIQ 
IS 
SISYA 
| | 
07 
Ö L / v 5 74 4 
AS1 ) 
Ä [ie 
ia Kr. 
I ) 11:11 rer sılJlırs ’ıS » ses S ir J 1) ng (/,0 M 
\ pa ven ki Fü de: Anschlusse tee-Gurt fuı Abb. 5. Einspannungsgrad als Funktion von o. 
verschiedene Werte von "1. (U) 





Man ersieht daraus das im ersten Augenbliek überraschende, aber bei näherer Über 
legung plausible Resultat, daß die Fortsetzung des Balkens über die Auflager hinaus als eine 
Art elastische Einspannung wirkt (s. a. Abb. 5). Der Spannungszustand außerhalb der Auf 
lager kann werzen Fehlens dort befindlicher Belastungen als Selbstspannungszustand bezeichnet 
werden. Ansehaulich kann man ihn so deuten, daß der Gurt den Steg gewissermaßen 
„heranholt*. 


Is ist zu vermuten, daß die Besonderheiten, die bei dem hier betrachteten Belastungs 
zustand des gelenkiz gestützten Trägers mit unendlich weit auskragenden 


linden auftreten: 
Kinspannungswirkung an den Auflagern 


vererößerte tragende Breite im gefährlichen Querschnitt gegenüber 
dem Fall des Trägers, der an den Stützen 
querschnitte besitzt, 


endet und starre End- 
sich auch bei anderen Belastungen ergeben. 


»Is 
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Vorträge der Hauptversammlung in Stuttgart 
der Gesellschaft für angewandte Mathematik und Mechanik. 


Von den auf der Versammlung in Stuttgart vom 26. bis 28. September gehaltenen 
Vorträgen brinzen wir im folgenden die von den einzelnen Vortragenden selbst zur Verfürung 
eestellten Vortragsauszüge. Die vollständige Veröffentlichung einzelner Vorträre bleibt vor- 
behalten. 


Bemerkungen zur Theorie statisch unbestimmter Systeme. 
Von W. Flügge ın Göttingen. 


Dem üblichen Verfahren zur Berechnung statisch unbestimmter Systeme mit Hilfe eines 
Hauptsystems haftet der Nachteil an, daß es bei hochsymmetrischen Systemen nieht immer 
gelingt, ein Hauptsystem zu finden, das alle Symmetrieeigenschaften des ursprünglichen 
Tragwerks bewalırt. Um hier Abhilfe zu schaffen, hat man eine Erweiterung des bekannten 
Gruppenlastverfahrens herangezogen, die darin besteht, daß mehr Schnitte geführt werden, 
als dem Grade der statischen Unbestimmtheit nach notwendig wäre, und daß aus den in 
diesen Schnitten anzubringenden zusätzlichen Kräften Linearkombinationen gebildet werden, 
von denen ein Teil zur Stützung des beweglichen Hauptsystems herangezogen werden muß 
und nur der Rest als überzählige Größen in die Elastizitätsgleiehungen einzeht. Die Auf- 
stellung soleher Gruppen bedeutet neben einem Verlust an Übersiehtliehkeit häufig auch eine 
weitere Belastung der Rechnung. Es soll deshalb eine neue Formulierung der Theorie statisch 
unbestimmter Systeme skizziert werden, die durch vollkommenen Verzieht auf den Begriff des 
Hauptsystems eine größere Beweglichkeit in der Methode ermöglicht. Zur Vereinfachung «des 
sprachlichen Ausdrucks werden wir uns dabei auf Raumfachwerke beschränken, obgleich 
sich der vorgetragene Gedankengang an sich auf jedes statisch unbestimmte System über: 
tragen läßt. 

Kın Raumfachwerk besitze % Knoten und m System- und Stützenstäbe. Im statisch be- 
stimmten Falle m 3% ist die Zahl der Stabkräfte gleich der der Gleichgewichtsbedinzungen 
an den Knoten. Sind die Lasten null, so ist dies Gleiehungssystem homogen und unter der 
mechanisch selbstverständlichen Voraussetzung nicht verschwindender Determinante sind auch 
die Stabkräfte null. 

In einem statisch unbestimmten Fachwerk haben wir einen Überschuß von » Stabkräften. 
Wir wollen » Unbekannte auswählen und alle Glieder, in denen sie vorkommen, auf «die 
reehte Seite schaffen. Dann können wir über diese » Unbekannten beliebig verfügen und die 
Gleichungen nach den übrigen auflösen. Wenn wir die Auswahl so treffen, daß die Nenner- 
determinante des verbliebenen Gleichungssystems nicht null ist, so ergeben sieh auch bei 
unbelasteten Fachwerk Werte für die Stabkräfte, die endlich und nieht sämtlich null sind. Ein 
solches Lösungssystem der Gleiehgewiehtsbedingungen, das zu einem unbelasteten Tragwerk 
gehört, nennen wir ein Eigenspannungssystem, während wir ein Wertsystem von Stabkräften, 
das für eine Belastung an jedem Knoten Gleichgewicht herstellt, ein Lastspannungs- 
system nennen werden. In einem statisch unbestimmten System sind also Eigenspannungs- 
systeme möglich, und ein elementarer Satz aus der Theorie der linearen Gleichungen besagt: 

Es gibt in jedem »-fach statisch unbestimmten Fachwerk » und nur a linear unabhängige 
Kigenspannungszustände. Alle übrigen sind daraus durch Linearkombination abzuleiten. 

Solche linear unabhängige Eigenspannungssysteme sind z. B. die Zustände A, I ın 
der üblichen Theorie. Man kann geeignete Eigenspannungssysteme Jedoch aueh 
noch auf andere Art gewinnen, und in der einfachen Ausnutzung dieser 
Möglichkeit liegt der Vorteil der hier vorgetragenen Darstellung. 

Kin weiterer einfacher Satz aus der Theorie linearer Gleichungen bezieht sich auf dıe 
Lastspannungssysteme: Jedes zu einer gegebenen Belastung gehörige Lastspannungssystem 
ist darstellbar durch Superposition irgendeines andern mit einem Eigenspannungszustand. 

Bezeiehnen wir also die Kraft im Stabe Nr. « bei n verschiedenen, linear unabhängigen 
Kigenspannungszuständen mit Sur, Su2, ...8ur, 2... 8um und mit Su. und Sao zwei Last- 
spannungszustände, so ist S.„ darstellbar in der Form 


N, — Sı. u" F: Nu v4 f . . . , . . . . . . R R (1 |. 


Wir können nun unsere Aufgabe, nämlich die Berechnung des statisch unbestimmten 
Systems, so formulieren: Unter allen dureh (1) dargestellten Spannungszuständen ist derjenige 
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auszuwählen. der wirklich eintritt. Als Kriterium zu seiner Auffindune benutzen wir den 
Satz von Castıieliano, der in unserer Terminologie heißt: 








Unter allen Lastspannungssystemen ist dasjenige, das wirklich eintritt, dadurch aus- 
gezeichnet, daß die im Träger gespeicherte elastische Energie ein Minimum ıst: 















| Sur Su 
NV. ) h) \ u 0 
E u EN 
(s Stablänee, EF Dehnunessteifiekeit). Führt man die Variation aus und setzt (1) ein, 
so erhält man Gleiehuneen des üblichen Tvps 
X ar l Br N ' N); a u gr A. N) ) (at RT ] Au Ö; n 1 dı 0 0 
mt 
v R Y S 
2) ' \ N } N ı? Or0 \ N } N ı . 
I Z lY 1 ee F, l 


wobei aber jetzt die d,, nieht mehr als Formänderungen in einem Hauptsystem definiert 
sind. sondern dureh die eben angeschriebenen Summen und die zugehörigen Eigenspannungs- 


systeme, ar 


Knickung von Rechteckplatten mit streifenweise konstanter Dicke. 
Von Konrad Ludwig ın Hannover'). 


Zunächst wird eine rechtecekizee Platte behandelt. die aus zwei Streifen mit den breiten 
ce und be und den konstanten Dieken f, und f, besteht, auf allen Seiten gelenkig gelagert 
ist und dureh konstanten Druck ın der Längsrichtung bean- 














7 1 T er 
B | 4 | sprucht wird. Die Differentialgleichung ıst 
| 
| | ’ 9 
k: r-— N Om Om Om / Om 0 
‘ { a 2 = 5 i ni, Pr . 
Ä j | VBT 0x0 rs Q Yy B "32 
Io 
| | FE f" 
B =. 
el l I2(1 4°) 


Die Randbedinzuneen an den zeedrückten Seiten werden durch den Ansatz 


N PER ER nn ER 
w‚—sin — IN, ini, y+ K,60] 4, y+ KR, <ın /),y+ K, ol 4, y], n 1.2, 


.. 
ww 


erfüllt. worın 





die Wurzeln der charakteristischen Gleichung sind. Die Randbedingungen an der Seite y=0 
ergeben A, 0, A,-0,. Entsprechend wird für den diekeren Streifen 


’ rn IT X e Fu FE 
n,— sin IN, <imtvr (by)-+K,Smv,(b— y)] 
2 | 


gesetzt, worin die Anzahl » der Halbwellen dieselbe ist wie beim dünneren Streifen, da die 
Stetigkeitsbedingungen für y> ein.r identisch erfüllt werden müssen. Die stetigen Größen sind: 


die Durehbieeune m". 
N 


ur 0 
der Drehwinkel, _, 
ce 


k 





Prof. Dr.-Ing Schleieher, Technische Hocehsehule Hannover, und Prof. Dr. Kaufmann. Technische Hoch 
schule Miinehen, danke ieh für wertvolle Ratschläge. 
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wir den . | Mm m 
das Biegungsmoment pro Längeneinheit G,-- blu 
| Gr Op) 
‘ch atıs- 
’ ol, ua oO’ m 0’ m] 
die Querkraft pro Längeneinheit N,+ , Be —- Win rl 
UaeT 2 | POLE OY oy| 


Diese vier Stetiekeitsbedingzungen ergeben für die vier Integrationskonstanten A, K,, A, und 
K. die vier linearen homogenen Gleichungen 
(1) ein, 





Sini,c-K, +&ind,c-K, — Einv,(b—c)-K, — Sinv,(b—c)-K,=d0, 
1, Eofi,c-K, +4,Cofi,c-K, +», Cojv, (b—c)-K, +r,&ofr,(b— c)-K,—=0, 
| | 
n SE ; u. m ıT ee a : 
rl ul +4) Sini,cK,+t’l— ul +4.) Sini,e-K, 
A |  ° A a ’ 
definiert 4 n ıı\? Br SS | n.ıaı\° Eee 
“ Hf,? | |—r,’ Sn, (b—c) KR, +, | v1 Sinv,(b—c)-K 0, 
INNUNgS- : A | 5 Zw | 2 
>47 
” nn\ , u RR . x n.nı\® a 
t.°’ (2 u)! A, —A’|ECofiA,c-K, +1,’ (2 «)| | / /,\boJ /,c- A 
| z \ da | Ad ’ 5 
icke. | er 


| nr? Ku | ı\” KR 
rt.” (2 | f |», —r’|Eofj», (b—c)-K,+t,’|(2 | N „vr, 1809,15 —c)-A , 
| < | a | j ”4q 


ı Breiten 











gelagert deren Koeffizientendeterminante 
\e bean- 
Igi,c Ig4,e 
A. /, 
‚| nn” nn 4 BZ rl Pier na 4 Ber 
mn 1) j oO zuA,C wi 1) - 0.1 204A,C 
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sein muß. Mit der dimensionslosen Veränderlichen 
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Thermo-elastische Platten-Gleichungen. 
Von K. Margwerre in Berlin-Adlershof. 


Die geringe Beachtung, welche die Probleme der Wärmespannungen im Rahmen der 
alleemeinen Elastizitätstheorie gefunden haben, mag z. T. daher rühren, daß eine systematische 
Übertraeune der sehon von Duhamel und von F.E. Neumann gegebenen drei-dimensionalen 
Spannungs-Dehnungs-Gleichungen auf die der Rechnung zugänglicheren zwei- oder ein-dimen- 
sionalen Kontinuen bisher fehlte. In einer (dem Verfasser erst vor kurzem bekannt gewordenen) 
Arbeit von Kent), auf die wegen weiterer Literatur verwiesen werden kann, ist die Wärme- 
beanspruchte Zylinderschale diskutiert: allgemeinere Schalenformen erfassen die Ansätze von 
Kichelbere?: in diesen beiden Arbeiten ist aber von den Neumannschen Gleichungen 
nicht Gebrauch gemacht. Die folgende Ableitung der thermo-elastischen Plattengleichungen, 
die den Zusammenhang zwischen den bisher unvermittelt nebeneinanderstehenden Gleichungs: 
systemen herstellt (und die ohne prinzipielle Schwierigkeiten auf die Schale übertragbar ist), 
läßt zueleich die Voraussetzungen, unter denen die Näherungseleichungen gelten, hervortreten. 


Bezeichnet (w,.r,.w,) den dureh die elastischen Kräfte hervorgerufenen Verschiebungs- 


vektor, (,. ",. ,) den dureh die Temperatur d hervorgerufenen, so gelten (mit 9; = div (a;, e;, m ;)): 


2 G 7 
Pr ' 1 . . . : 
i£ 59, +20, und zwei weitere Gleichungen, 
m y (.r 
(du, ‚0a ns 
xy (i | \ de und zwei weitere Gleichungen, 
0 VA 
DET 6; ou, ou, a 
z = @, - =D, und sechs weitere Gleichungen. 
(a .) 04 04 


Diese beiden Gleiechungssvsteme lassen sich mit (n,, 0, m,) = (t,,r,,w 


0 1) 1 9 N Him,, Us, ',,) ZU° 
sımmenlfassen zu: 


l 


2 G Q 0% 0 
0, « - 2 U ‚U, ISW, 
. 1 2 ” Or I 
ch). 
7 (} (* H ( a) USW 
® 1 2. - . 
v/ ( „ ( ae) 





In dem (Neumannschen) Gleichungssystem (1) mögen die elastischen Konstanten E,G, m und 
das Produkt aus dem kubisehen Ausdehnungskoeftizienten e und dem Kompressionsmodul: 


( m+I em 
ee et 


> m —& 3 m — 
als von der Temperatur unabhängig angesehen werden. 

Dann gestatten die 4 ersten Gleichungen dieses Systems zusammen mit den unver- 
änderten Cauchyschen Gleiehgewichtsbedingungen und den übliehen, an die Stelle der beiden 
letzten Hooke-Neumannschen Gleichungen tretenden, Hypothesen der Plattentheorie: 


Om nETE h 
[A [A Fr r f " 4 . . ; ä P ; \ 1 (2) 
’ V.r ® a7, 
die Herleitung der thermo-statischen Plattengleichungen. [u,”, 0 sind die von 2 unab- 
häneieen Verschiebungen der Plattenmittelebene z2=0.| 
. h ich dr, 
löst man die dritte Neumannsche Gleichung nach | auf: 
Ü< 
Om, pP m—2 q | (‘ u .0% m —?2 
r ! T er . 0 
> 2>Gm—| nm - 1\dx 0y > Gm A) 


und setzt diesen Wert in die beiden ersten ein, indem man bei dieser Substitution die Normal- 
spannung o, (die in diesen Gleichungen nur den Charakter eines Korrekturgliedes hat) streicht, 
so lautet das System der drei für die Plattentheorie verbleibenden Neumannschen Gleichungen: 


I) Thermal stresses „.. Trans. Am. Soe. Mech. Eng. 1951, S. 167. 
>, Temperaturverlauf ... Verbrennungsmotoren. Forschungsarbeiten VDI, Heft 265 (Berlin 1023). 
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einzuführen (27 Plattendieke). 


bedingungen erhält man am einfachsten aus den Cauchyschen Gleichungen durch Integration: 









und unter der 
3) und 16): 


man dureh 


„ersten" bzw. 











m | 


der örtlichen Spannungen 
Mittelwertsbildungen 


M, 
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loxzdz, 
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Plattentheorie 


G. 
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(4) 


Die zwischen diesen „Resultanten“ geltenden Gleichzewichts- 


\ "dz, 


h 


. N 


Sind Ober- und Unterseite der Platte spannungsfrei, 
einfacher schreiben: 


x 


Q” 


ARHNETEE 


„zweiten Mittelwert der Temperatur“ ein. so Ist aber 














Om 


vermöge (2 


die beiden ersten Gln. (5) 





rleiehen Voraussetzung ergibt die Elimination der Querkräfte 9, und Q, aus 


wi, 


‚’ 
(9) 


und (4): 
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V jr m | R N) ") 
N. 4 ® Me 1 pP : 
a 2m Ooy dw 
) 5) Y N (5): 
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M, E—W|  ;,+— x, >), 
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m 0’ w Il 0°w ce > 
M, E—kW|- +; >), 
' N 7, mor0oy Ah 
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und aus (7) und (5) ergibt sich das System der thermo-statischen Plattengleiehungen in der 
Form: 





Ö & u. od ") m 10 (“ » du ce 09 
Ooadr 0yY 2m Oy\x 0 > 3 0#' 
L N, 
Q ? u 0 5 m—1 0 (° vo "\ 7, 
oyoz ody 2m dal Om 30Y 
"mw, ,„ 0'w Om 7‘ 2 N) 10 
ig a 5) 1) ce in &- ı" d \ ! . . . . . . . ° . ( % 
0. Oro oO‘ h 


An Stelle der in diese Gleichungen eingehenden Mittelwerte 9, 9 werden im allgemeinen 
nur die Temperaturen T, und T, auf Ober- und Unterseite der Platte unmittelbar gegeben 
sein. Mit Hilfe der Fourierschen Wärmeleitgleichung 4,9 ==0, die für die stationäre Wärme- 
beweeung im Innern eines (dreidimensionalen) Kontinuums gilt, lassen sich aber leieht Nähe- 
rungsdifferentialgleichungen für den Zusammenhang zwischen diesen Größen herleiten. Durch 


h h 
Ausübung der Operation |...dz und J...zdz auf A,#=V0 ergibt sich: 
h h 
h h h 
af a. 7? N) 
( ( ed 
‚\Vdz + _ \ Vılz 0), 
u \ . r 02 
h h h 
h h 
32 ( y2[ 9 Ih 
( ( 0A 
‚\9zdz+ - \ Vzdz-+- 1: DA 0 
Or? \ Omy?. | 02 | h 
h hı 


Gilt zwischen der Plattentemperatur 9 und der Außentemperatur T die Newtonsche UÜber- 
gangsbedingung: 
00 EN 
ho al 9) 
!N o,u 0, U 


(7 Leitwert, a Übergangswert, » äußere Normale), so kann man schreiben: 


h -AAdB-ad,+9)=-alT,+T,); | 


, (ih), 
N -AAd-Utahl9,—9)=—-ah(T,—T,) | 


In einer hinreichend dünnen Platte kann der Temperaturverlauf mit z als näherungsweise 
linear angenommen werden’); dann gilt: 


9, +9.,=29 und 9, —-du=69, 


3) Dureh diese Annahme (und die aus ihr folgenden GlIn. (12)) verziehtet man auf die Befriedigung der (räumlichen) 
"ouriersehen Gleiehung JS, I 0 renau so, wie man in der Plattentheorie den Gln. (2) zuliebe die >. und 6. der 
Iiookeschen Gleiehungen aufgibt. 
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und die Gln. (11) gehen mıt der Abkürzune : . über in: 
/ 











[0 | 
Sind die Außentemperaturen T mit = und y nicht allzu stark veränderlieh I 4 „ 7). so 
(X ! j 


eeoen das zweite vernachlässieen : 








kann man das erste Glied auf der linken Seite von Gl. (12. 2 







(13): 


die Platteneleiehunge (10) lautet dann: 






OT) A (1). 









mit den Gln. (8) für die Biegxemomente. Gl. (12, 1) (die niehts anderes Ist als die auf die Platte 
erweiterte „Wärmeleiteleiehune des dünnen Stabes") und mit ıhr die Gln. (9) lassen eine 







ähnliche Vereinfachune nicht zu. 
Wie an anderer Stelle eezeiet werden soll, tritt bei Verschiedenheit der UÜbergzanes- 
koeftizienten «a, und «a, auf den beiden Plattenseiten an Stelle der Gl. (15) «die Beziehung: 






ı) € - 
aA J 


B a,h + a,h | 
l 





Uber einiee andere Modifikationen, die dureh diesen (in den Anwenduneen häufieen) Umstand 
entstehen. sowie über die Übertraeune der Gleiehunzeen auf die Schale und ein Anwendunes- 







beispiel wird noch beriehtet werden. >49 


Über Ausbreitungs- und Anlaufvorgänge einer 
zähen Flüssigkeit mit zylindrischer Begrenzung. 
Von Wilhelm Müller in Aachen. 









I. Es soll sieh im folgenden um laminar-instationäre Bewegungen einer zähen Flüssıekeit 
handeln, die außer von der Zeit nur von einer Koordinate abhängen und die sıch als Lösungen 
der Stokes-Navierschen Differentialgleiehungen bestimmen lassen. Die beiden Hauptfälle 
sind die lineare rotationssvmmetrisch angeordnete und die ebene zirkulare Bewegung, deren 
zuzeordnete Geschwindiekeitsfunktionen den linearen (dimensionslos geschriebenen) Gleichungen 








(h). 







7 \ -- |» z \ N . ; A . s } (2) 
c) I” r or r" Er "A 4 ( / 







oenüzen. Man wird zweekmäßig diejenigen Bewegungen, bei denen das Druckgefälle ver- 
schwindet, als Ausbreitungs- (Diffusions-) Vorgänge, die unter der Einwirkung eines eingeprägten 
Druckeefälles zustande kommenden Bewerungen dagegen als Druck- oder Anlaufströmungen 
hezeiehnen. Neben den früher ausführlich behandelten freien Ausbreitungsvorgängen spielen 
besonders die erzwuneenen oder randgzeführten Ausbreitungsvorgänge eine Rolle, die wesentlich 
eekennzeiehnet: sind dureh eine vorgeschriebene Bewegung des Randes. 










>. Betrachten wir z. B. den der ersten Gruppe angehörenden Fall, daß die Flüssigkeit 


ebene kreisförmize Bahnen beschreibt, während die zylindrische Grenze mit der Ge- 
yır) 





schwindiekeit ("), ,. gl)» fit) rotiert. Unter der Voraussetzung, daß Jg ir) je 


also = == (r) einem stationären Zustande entspricht, setzen wir 
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Bestimmen wir also v, so, daß die Gleichung 


| v. au; oe, Ä 
F r yiır) Ji dr a a ee Ge Ve er  \ 
: - ! / Of 
' 2  #& , op 
erfüllt ist, so genügt die Summe r=",+r, der Gl. (2) mit ng. 0. Wenn z. B. die 
04 


Flüssigkeit das Innere eines Zylinders mit dem Radius a1 ausfüllt, hat man y(r)= or 
zu setzen. Um die entscheidende Lösung », zu finden, entwickelt man orf’(h) in eine Reihe 
mit Besselschen Funktionen .J,(4;r), wobei f als Konstante aufzufassen ist und die Para- 
meter /; als Wurzeln der Gleichung J,(/) 0 zu bestimmen sind. Führt man ferner für », 
die Reihenentwieklung 
un „ 
®. PB R;(t)e WET (Ar) 


N | 


ein, deren Koeffizienten von der Zeit abhängen, so erhält man dureh Einsetzen in die Gl. (#) 
für jeden Zeiger i eine Beziehung, aus der B; dureh eine einfache Quadratur sich berechnen 
läßt. Daraus entsteht dann unter der Voraussetzung der Konvergenz die Gesamtlösung 


% t 
\ v7 (4; Y) ):24 i 2: 
ı yir)- fit)+- 2 in _e ° "(det "dı (N) 
! A . Dr 1; ,(4;) \/ 
=) 1 


Setzen wir z.B. f)=1-e ‚ was einer anfangs beschleunigten und später angenähert 
konstanten Rotation des Zylinders entspricht, so ergibt sich 


SF. 
9 ? a .J (/;r) 24 2 
onril u 9m \ ei (ei " de 2..5: 
u) U ki Ki J,(/;) 
i | 
als Ausdruck für den Übergang in den stationären Bewegungszustand "= wr. Wächst «a’ 


unbegrenzt, so erhält man die Lösung für den Fall einer zur Zeit #=V0 plötzlich einsetzenden 
und konstant bleibenden Rotation des Zylinders. Nach derselben Methode lassen sieh auch 
die erzwungenen Sehwingungen berechnen, die von einer vorgegebenen Drehschwingung des 
/ylinders erzeugt werden. Ist der Zylinder von zwei konzentrischen Kreisen begrenzt, so Ist 
statt der Funktion ./, (7; r) eine aus Besselschen Funktionen erster und zweiter Art zusammen- 
sesetzte Zylinderfunktion zu verwenden. Wenn die Flüssigkeit endlich das Äußere des 
kreiszylinders ausfüllt, so erhalten wir statt der Summe ein Integral aus Hankelschen 
Funktionen, das für den Fall der plötzlichen Inbewegungsetzung des Zylinders bereits von 
S, Goldstein!) anzezeben Ist. 


.) 


3. Die andere wiehtige Gruppe stellen die Übergangs- und Anlaufvorgänge in einen 
geraden Rohr mit kreiszylindrischer Begrenzung, die unter der Wirkung eines konstanten 
oder zeitlich veränderlichen Druckgefälles sich entwickeln. Ich habe im besonderen in Er- 
weiterung der Ansätze von Szymanski’) und Vogelpohl’) den Fall einer axialen Strömung 
in einem Hohlzylinder mit den Radien r=b, r=a-=-1 als Lösung der Gl. (1) untersucht. 
Kür die Reihenentwieklung des Geschwindigekeitsausdrucks benutzt man hierbei zweekmäßiger 
Weise die lineare Verbindung der Besselschen Funktionen ./, und Y, erster und zweiter Art 


Y J(ir) Y,(Ar) S,(ır, Ab) 70 S,(Ar, Ja) R 
tar) or ——; bzw. Z (ir a (4 
Jh) Ylab) J,1zb)Y,(zb) | Jtza)Y,iza) 

wobeı die Parameter 4 als Wurzeln der Gleichung 
S,Aa,/b)=J, (ka) YAb)-Y,AaJ) Ad)=0 .. . 1.2.0202... (8) 


zu bestimmen sind. Daun verschwindet in Übereinstimmung mit der Haftbedineunge die 
Funktion Z,(ir) für die Grenzen r=b und r=a. Wenn b sich der Null nähert, so zeht 
Z,(4r) in die gewöhnliche Besselsche Funktion J,(Ar) über. Da in einer besonderen Ver- 
öffentlichung die Rechnungen ausführlich entwickelt und die Ergebnisse auch zahlenmäßig, 
insbesondere für den Fall, daß die Anfangsgeschwindigkeit und das Druckgefälle konstant 
sind, vorgeführt werden sollen, so kann auf eine Wiedergabe der Ergebnisse in diesem 
/usammenhange verzichtet werden. Für andere Begrenzungsformen, vor allem für den Fall 
des elliptischen @Querschnitts, sind die Rechnungen vorbereitet. DU) 


') Proe@e,. London Math. Soe., Vol. 34 (1932), S. 51. — ?) Journ. math. p. et appl., Paris 1932, 8.67. 3) ZAMM, 103%, 8. 46. 
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Die Randbedingungen der Hauptspannungs-Stromlinien. 


Von 






H. Neuber in München. 





ls gıbt bekanntlich eine Reihe von ebenen Spannungszuständen, bei denen die Haupt- 
spannungslinien zugleich als Stromlinien einer ebenen Potentialströmung aufgefaßt werden 
können. Die zugehörigen Beziehungen, welche den Zusammenhang zwischen Spannungen und 
Spannungslinien angeben, wurden zuerst von P. Nemenyi') abgeleitet. Neuerdings hat 
U. Wegner’) diese Frage sehr ausführlich behandelt. Beide Forscher haben sich jedoch 
mehr mit der allgemeinen Formulierung befaßt und insbesondere die Frage der Erfüllbarkeit 
vorzerebener Randbedinzungen noch so zut wie unbeantwortet gelassen. Diese Frage, welche 
für die Anwendbarkeit der Stromlinien -Spannungstheorie von grundlegender Bedeutung Ist, 
bildet den Gegenstand der vorliegenden Untersuchung. 









Sind @r, y kartesische Koordinaten und «,” die Koordinaten der Hauptspannungslinien, 
so können wir die Gleiehgewiehtsbedingungen in der Form 












00,,0, 0,0 h. 0 00, , 0, 00h, 0 
t— - on 2 - u - ° ö o . . . R (1) 
a7; h„ ou Or hı 0% 











schreiben. Hierin sind A, und Ah, die Verzerrungsfaktoren, welche sich aus 
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ergeben und den Abständen der Hauptspannungslinien proportional sind. ©, ist die Haupt- 
spannung normal zur Spannungslinie » = konst., entsprechend o,. 






Am lastfreien Rand » = konst verschwindet o,. Aus (1) folgt dann 
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Am lastfreien Rand ist die längs des Randes wirkende Spannung stets 
umgekehrt proportional dem Abstande der dem hand benachbarten 






Hauptspannungslinie. 





Zu diesem Satz, der ganz allgemein gilt, kommen bei den Hauptspannungslinien mit 
Stromliniencharakter noch weitere Bedingungen hinzu. Es bestehen dann zwischen «,r und 
“,9 die Cauchvy-Riemannschen Gleichungen, so daß h„=h,„=h wird (Isometrie). Aus (l) 
und der Verträglichkeitsbedingung 










ie: ee see 


folet dann dureh Integration die Nemenyische Lösung 
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Durch geeignete Wahl des Ursprungs können wir hierfür auch 












ar db 
h? 
setzen. 
Am lastfreien Rand » = konst verschwindet o,, so daß 6, der Bedingung 
ar +b en 
Ö, - (0 
h’ 





eenürzt. Andererseits muß (3) erfüllt sein. Eliminieren wir o, aus (7) und (3), so ergibt sich 
als Randbedingung für den Verzerrungsfaktor h: 


ar“ + l 


h 


1), ZAMM 13 (1933), S. 364. 
2, Ing.-Arch. 5 (1934), S. 449. 
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Bei einer lastfreien Hauptspannungs-Stromlinie gilt demnach eine doppelte Rand: 
bedingung: Einerseits muß die Stromfunktion längs des Randes konstant sein, andererseits 
muß ihr Gradient einem bestimmten Gesetz folgen. 

Setzen wir ferner (8) in (3) ein, so erhalten wir die Randbedingung für die im Rand 
wirkende Spannung: 


e® 


0, _- 


— ur : (+). 
ar’—+b 

Nurdann, wenn die im lastfreien Rand wirkende Spannung einen Verlauf 
hat, der diesem Gesetz entspricht, können die Hauptspannungslinien als 
Stromlinien aufgefaßt werden. 


Diese Bedingung ist, da sie vom Koordinatensystem frei ist, außerordentlich ansehauliech. 
Bei spannungsoptischen Messungen insbesondere läßt sich sehr einfach feststellen, ob diese 
Bedingung erfüllt ist, da Ja die Randspannung unmittelbar gemessen wird. 

Verschwindet die im Rand wirkende Spannung, so folgt aus (3): e=0. (8) liefert, 
mit e multipliziert, 


BEER : 3 ea a dh re 


En A u  - : | 


Beim ebenen Stromlinienspannungszustand kann am lastfreien hand 
dieim Rand wirkende Spannung nur dann verschwinden, wenn es sich um 
einen kreisförmigen Rand handelt. Für h besteht in diesem Falle keine Rand- 
bedinzung. | 

In einer ausführlicheren Arbeit des Verfassers’) werden diese Gesetzmäßigkeiten an 
Hand von Beispielen in bipolaren Koordinaten näher erläutert. Es handelt sich dabei z. T. 
um Probleme, welche bisher mit der Airyschen Spannungsfunktion gelöst wurden, für die 
nun aber mit Hilfe der Spannungs-Stromlinientheorie eine überraschend einfache Behandlung 
möglich wird, so z. B. der von P. Fillunger*) für die Bereehnung des Eisenbahn- 
lasthakens zugrunde gelegte Fall, ferner die gedrückte Walze; dann auch das 
allseitie gezogene Blech mit zwei kreisförmigen Löchern (Aırysche 
Spannungsfunktion von Th. Pöschl’) und ©. Weber®)). Schließlieh läßt sich auch leicht 
zeigen, daß es sich bei der symmetrischen Außenkerbe schwacher Krümmung 
bei Zugbeanspruchung (Airysche Spannungsfunktion vom Verfasser ’)) um einen Stromlinien- 
spannungszustand handelt, was für die Kerbforschung von Bedeutung ist. Sal 


Über eine direkte Methode zur Lösung von Aufgaben 
der Plastizitätstheorie. 


Von Th. Pöschl in Karlsruhe. 


Das mathematische Problem der Plastizität ist bekanntlich in den Spannungen statisch 
bestimmt, in den Verschiebungen aber statisch unbestimmt. Man gelangt zu Widersprüchen, 
wenn man für eine vorgegebene Aufgabe die Lösungen der Spannungsgleichungen allein hin- 
nimmt, ohne sich um die Verschiebungen zu kümmern, die zu jenen Spannungen gehören. 
Die vollständige Behandlung einer Aufgabe erfordert eine vollständigere Formulierung, wie 
sie im Anschluß an die bekannten grundlegenden Arbeiten — in sehr übersichtlicher Weise 
etwa in der Dissertation von R. Lohan (Berlin 1955) gegeben ist. 

Unter den üblichen Annahmen (ebene, beschleunigungsfreie Strömung ohne Massenkräfte) 
hat man die (sieben) Gleichungen 

are, eat, are eis 5, air 
für die (sieben) Unbekannten d (zwei Komponenten), $’” (drei Komponenten), p und 7 auf- 
zulösen:; ferner ist die Hinzunahme geeigneter Randbedingungen erforderlich. 

Als eine Methode zur Lösung dieser Gleichungen hat Frau H. Pollaczek-Geiringer 
in ihrem Stockholmer Vortrag (1930) die folgende verwendet: man geht von einem in statischer 
Hinsicht gelösten Problem aus, sucht die zu diesem gehörigen Verschiebungen und die diesen 


3) H. Neuber: Ing.-Arch. 6 (1955), S. 325. 
1) ZAMM 10 (1930), S. 221. 5) ZAMM 1 (1921), S. 174. 6) ZAMM 2 (1922), S. 185. ’) ZAMM 13 (1933), S. 430- 
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angepaßten Randbedingungen. Dabei wird die Gleichung für die Stromfunktion in der auf 
die Gleitlinien bezogenen Form verwendet, die sich erst aus der Lösung der Spannungsaufgabe 
ereibt. Es ist dabei nicht zu sehen, wie man vorzugehen hat, wenn es sich um die Lösune 
der plastischen Aufgabe für einen vorgegebenen Bereich handelt. 






lelı möchte deshalb für die direkte Lösung dieser Aufgabe einen Wee vorschlagen. der 
einem in der Potential- und Elastizitätstheorie seit langem bekannten Verfahren entspricht: 
er besteht darin, die Gl. (1) unmittelbar auf orthogonale krummliniee Koordinaten &. y„ zu 
translormieren, der Diverzenzzleiehung div d=0 in bekannter Weise durch eine Strom- 
funktion „> zu genügen und solehe Lösungen zu suchen, bei denen y' nur von einer Veränder- 
lichen (etwa ») abhängt; die Stromlinien oder die Linien des plastischen Fließens sollen 
»  konst. sein. In manchen Fällen wird gemäß der Beschaffenheit des Problems die Form 












der Stromlinien von vorneherein nahe liegen, so daß das Koordinatensysten, auf das die ganze 





Untersuchung bezogen wird, den geometrischen Bedingungen der Aufgabe angepaßt wird: auf 
diese Weise zelinet es dann in gewissen Fällen, zu genauen oder weniestens aneenäherten 
lLösunzen zu eelanzen. leh zebe hier die Ausführung dieses Gedankens für den Sonderfall 


isometriseher Koordinaten an, die Erweiterung auf beliebige, orthogonale krummliniee Ko- 
lolet 











ordinaten dann unmittelbar. 






Setzt man für das Bogzenelement 
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erfüllt dureh Einführung von y', sobald 















ist. Die Gleichung 7... Ak lautet, in den Verschiebungen ausgedrückt, 
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eine Gleiebune. die auch so zeschrieben werden kann 
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in der also die Hessesche Determinante von y nach ., y auftritt. 






Wenn man diese Gleichung auf die isometrischen Koordinaten £&, 7 transformiert, und 










cemäß dem oben dargelegten Verfahren die Funktion y' nur von „ abhängen läßt, also 
vs =0U setzt, so nimmt sie die Form an 
N 7,2 
um a” kin,” N INK NyNzylypnW x u), BL... (7) 
/ I, I I/a 7 Tr I 1 „ x x KU YTIKY / 7 / 7 Z Ö (»#. y) / 7 un 3,7% . . . . 











Die Gleichung grad p=1 P’ lautet in den Koordinaten 5, » 





Io» 0 0: O TE Iop Ö TE, OO 0: 


= rt re 
7 ss u (3) 
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Als erste Anwendune dieser Methode betrachten wir den Winkelraum zwischen zweı 
festen Grenzen und nehmen an. daß die Stromlinien mit den Radien zusammenfallen: die 
Stromfunktion hängt dann vom Polarwinkel „=g allein ab. Man setze also 


ge tn In: In» iq, 5 "COSQG @COSg, „ rsing eSsıng und u, zu f, 


dann nimmt die Gl. (7) die Form an 





“> " . e) ( 
e DT r I] > odeı 2 = Fr 7 ’ i ; r : } - ()). 


34 vavyvi’+4f’ 
Benutzt man die Gleichungen gradp= I W, W’=49' oder die Gl. (8), so folgt (dureh Elimi- 


nation von p) für f die Gleichung 
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Eine partikulare Lösung dieser Gleichung kann man unmittelbar angeben: f= Ae""", mit 
e+4=1/a? Sie führt zu folgenden Ergebnissen : 
. . .) D 

Z 4 l Ad 2, e(gq fo) r BR 7 . 3k 3 11 

/. ay3’ € ; &% 0} , Tr v5° Er }, 
vibt also konstante Spannungen, und 

| ck fo ) 
0, „€ ‚, =) (12). 


Die Konstante A stellt im wesentlichen die „Ergiebigkeit“ der im Koordinatenanfangspunkt 
sitzenden „Quelle* dar. Zur Festlegung von e ist noch eine Aussage über die Verteilung der 
Geschwindigkeit für ein bestimmtes r erforderlich. Die Gl. (10), die „ selbst nicht enthält, 
kann durch die Substitution f=y(f) auf eine Differentialgleichung 1. Ordnung zurückgeführt 
und vollständig integriert werden. 


Es ist zu bemerken, daß die sonst für den Keilraum angezrebene Lösung (siehe z. B. 
A. Nadai, Der bildsame Zustand der Werkstoffe, S. 149 u. f.) von der Form 


am —3ko+komsl, wı=h . » 2-3: 3% 050% KR 
für die Verschiebungen den Ansatz liefert 
Ur Fig ), ©, Fiq Bu. i a ; ı i b : i i ! - : ; (1 h), 


wobei die Funktion F vollständig willkürlich bleibt. Für den Keilraum zwischen festen 
Wänden muß man verlangen, daß v, im ganzen Bereiche einerlei Vorzeichen hat, während vr; 
für zwei verschiedene x verschwinden muß, was für die Ansätze (14) unmöglich ist. Diese 
l,ösung ist daher für den Keilraum nicht brauchbar. | 


Ein anderes Beispiel wäre das beim Walzen vorliegende Plastizitätsproblem; man müßte 
dann die Gl. (1) auf Kreiskoordinaten (bipolare Koordinaten) transformieren mittels der Funktion 


l, ic—z En 
a ee a 


a a 
strtn: i 
die Kreise &=-konst. bilden ein Büschel mit den „reellen“ Grundpunkten +ie, die Kreise 
; =konst. ein Büschel mit diesen als „imaginären“ Grundpunkten. Zwei Kreise 7=-+n, 
sollen die Walzen darstellen. Um zu einer nur von 7 abhängigen Stromfunktion = (1) 
zu gelangen, muß man £ als klein ansehen, oder cos > 1 setzen, was Ja den beim Walzen 
vorliegenden Verhältnissen in der Nähe des engsten Querschnittes entspricht. Die bezüglichen 
Rechnungen sind noch nicht abgeschlossen. 

Die Einführung anderer geeigneter Koordinatensysteme zur Lösung weiterer Aufgaben 
sowie die Erweiterung auf räumliche Probleme, insbesondere auf den rotationssymmetrischen 
Fall, ist naheliegend. 552 


Erzeugung von Oberflächenwellen durch schwingende Körper. 
Von M. Schuler in Göttingen. 
(Aus dem Institut für angewandte Mechanik der Universität Göttingen.) 


Zuerst wurde die Versuchsanordnung geschildert, mit der man die Tauchschwingungen 
eines Körpers an der Wasseroberfläche untersuchen konnte. Die ersten Versuche ergaben, 
daß bei kleinen Frequenzen und Amplituden laufende Wellen von der Frequenz des Erregers 
entstehen, während bei Steigerung der Frequenz und der Amplitude ein Umschlag zu stehenden 
Wellen eintritt. Dieser Wellenumschlag ist in dieser Zeitschrift 1933, S. 443 bis 446 behandelt. 
In dem Stuttgarter Vortrag wird die laufende Welle im stationären Falle und für das 
ebene Problem behandelt. Für eine schwingende Linie an der Oberfläche und für eine 
periodisch atmende (Quellinie in einer Tiefe A unter der Wasseroberfläche kann man die er- 
zeugten Wellen streng berechnen. Für eine an der Oberfläche des Wassers schwingende 
Platte von endlicher Breite wird eine Näherungsreehnung zur Ermittlung des erzeugten Wellen- 
systems angegeben. Für einen schwingenden Quader, dessen Boden eine Strecke h unter der 
Wasseroberfläche liegt und der endliche Breite besitzt, werden zwei Näherungsformeln für die 
Berechnung der erzeugten Wellen angegeben. Die durch diese Formeln gefundenen Näherungs- 
werte wurden mit den experimentellen Meßergebnissen verglichen. Es zeigte sich, daß die 
Näherungslösungen in bestimmten Frequenzbereichen brauchbar sind. a9) 
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Zur hydrodynamischen Theorie der Lagerreibung. 
Von @. Vogelpohl VDI in Berlin. 


Aus dem Institut für Technische Strömungesforschunge an der Technischen Hochschule Berlin.) 






Die bisherigen Ergebnisse der hydrodynamischen Theorie der Schmiermittelreibung 
sichern ihr noch nieht das Anwendungsgebiet, das ihr namentlich bei der Ölprüfung zukommen 
müßte. Die von Reynolds aufgestellte Differentialgleichung für den Druck p einer in der 
X-Z-Ebene liegenden Olschieht von der Stärke h lautet 








I Bu . „Muri a ee 
VE\NONx 07\n 07% ( 

wobei », die Zähigkeit des Schmiermittels und U die Geschwindigkeit der Begrenzung für 
y° h bezeichnet. Die Randbedingungen sind p=V0 für den Rand eines hier als rechteckig 
angenommenen Bereichs, die Seiten seien 2=0 und z=a bzw. z=0 und zb. 

Bei der Integration bereitet vor allem die Berücksichtigung der endlichen Lagerbreite 
Schwierigkeiten, so daß der Druckverlauf in der Achsriehtung meist von vornherein ange- 
nommen wird, entweder als Kosinuskurve (Gümbell, als gemeine Parabel (Schiebel, 
Karelitz) oder als Parabel mit gebrochenem Exponenten nach Versuchen zwischen 2,2 und 
27 (Nüeker) Eine genauere Berechnung, wie sie von Michell und Duffing für den 
ebenen Schmierkeil vorliegt, besteht für das endliche Zapfenlager nicht, ist aber zur Be- 
urteilung des Schmiervorganges ım wirklichen Lager unerläßlich. Nach dem gegenwärtigen 
Stand der Theorie ıst man nicht imstande, bei einem Vergleich zwischen Reehnune und 
Beobachtung zu beurteilen, ob auftretende Abweichungen den Voraussetzungen der Theorie oder 
Ihrer unvollkonmmenen mathematisch-reehnerischen Durchführung zugeschrieben werden müssen. 

Für eine weitergehende Behandlung war maßgebend, daß (1) eine selbstadjungierte 
Ditferentialgleiehung ıst, die sich aus einem Variationsproblem herleiten läßt: bei einer 
solehen Fassung ist dann das für die numerische Auswertung sehr geeienete Verfahren von 
kitz anwendbar. Die Frage nach dem physikalischen Sinn der Variationsaufgabe ergibt den 
interessanten Satz: Die Druckverteilung in der Schmierschicht eines Lagers 
stellt sich so ein, daß die dureh Reibung verzehrte Energie ein Minimum 
wird. Setzt man «die dissipierte Energie in bekannter Weise an und macht dieselben 
Voraussetzungen wie Reynolds über die Größenordnungen der auftretenden Differential- 
quotienten der Geschwindigkeitskomponenten, so erhält man unter Hinzunahme der Inkom- 
pressibilitätsbedingung 
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und als Ausdruck für das Eintreten eines auszezeichneten Wertes die Gl. (l) als Eulersche 
Gleiehune. Daß es sich um ein Minimum handelt, geht aus den Sätzen von Helmholtz 
und Korteweg über die Strömung einer zähen Flüssiekeit bei verschwindenden Trägheits- 
eliedern hervor, nach denen im stationären Zustand der Strömung die durch Reibung verzehrte 













nergie ein Minimum wird. 

Für die Auflösung des linearen Gleichungssystems zur Bestimmung der Beiwerte bei 
der Anwendung des Ritzschen Verfahrens wırd eine Art Orthogonalisierung angegeben, 
wobei aus dem zur Rechnung benutzten Funktionensyvstem Polynome gebildet werden. Mit 
diesen erhält man für jeden Beiwert eine lineare Gleiehung mit einer Unbekannten. Auf 
solehe Weise lassen sich die Unbekannten des Gleichungssvstems der Reihe nach berechnen, 
man hat die Möglichkeit der Probe nach der Bestimmung jeder neuen Unbekannten und einen 
uten Anhalt für den Grad der Konvergenz aus der Abnahme der Beiwerte. — Diese Ver- 
meidung der unmittelbaren Auflösung des Systems linearer Gleichungen zur Bestimmung der 
Beiwerte ist nicht auf den vorliegenden Fall beschränkt und ganz allgemein für das Ritzsche 
















Verfahren anwendbar. 

Die Reehnung wird durchgeführt unter Annahme einer nach Sinusprodukten fort- 
schreitenden Doppelreihe, die die Randbedingungen für eine rechteckige Gleitfläche erfüllt. 
Bei dem von Duffine behandelten Fall linear veränderlicher Zähizekeit ergibt sich eine sehr 
einfache Reehnung, schwieriger ist der von Michell gerechnete mit konstanter Zähigkeit 
des Schmiermittels. Ein Vergleich mit der exakten Lösung Michells zeigt die Brauchbar- 
keit des Ritzschen Verfahrens. Die Behandlung des technisch wichtigen Falles der Halb- 
schale wird in Aussicht gestellt. Auch eine Berücksichtigung der veränderlichen Zähigkeit 
oder der Durehbierzungen von Welle und Schale ist möglich, sie erfordert teilweise nur einen 
entsprechend größeren Aufwand an Rechenarbeit. Einzelheiten wird die endgültige Ver- 
4 





öffentlichung bringen. 
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Theoretischer Ansatz für die Schulerschen Umschlagwellen. 
Von €. Weber in Dresden. 


Die Schulerschen Umschlagwellen sind dargelegt in seinem Vortrage in Würzburg'). 

Gegeben sei ein Kanal mit rechteckigem Querschnitt. Die «-Achse gehe «er zum Kanal, 
(die „Achse von der Oberfläche an nach oben, die z-Achse in Längsrichtung des Kanales. An 
der Stelle 20 denken wir uns eine Querwand. Es soll gezeigt werden, daß es Flüssigkeits- 
bewegungen gibt, für die diese Querwand sieh periodisch hin und her bewegt, wobei keine 
von dieser Wand fortwandernde Kanalwellen, sondern die Scehulerschen Umschlagewellen 
auftreten. 


Die Bewegung sei eine Potentialbewegung einer idealen Flüssigkeit. Für die Potential- 
funktion P (a, y,2,t) gilt 


2 L3 Bi. 0 06 000 edel Lin > ER 
Aus ihr erhält man die Geschwindigkeit vd Grad». 
Die Randbedingungen sind 
a» ab 
vd für -=0 und zb. od) für y h. 
Or Ö „ 2 


Die obere Randbedingung ergibt sich wie folgt: Aus der Energiegleichung erhält man für 
ie Oberfläche, da dort der Druck p 0 ist, 


I /oP | 


Up | N .) 


2 SET TTE (2), 


aus der Geschwindigkeit eines Oberflächenteilchens 


t 
Ur | "mat dt- Yo: r : : (>). 


ı) 
Damit ein gewähltes PD, das der Gl]. (1) genügt. eine riehtiee Lösung ıst, müssen Gl. (2) 
und (3) für y dieselben Ausdrücke geben. Hierbei ist jedoch zu beachten, 


I. daß in der rechten Seite von (2) auch „ auftritt, welches zu beseitigen Ist, 


2. daß bei Inteeration von "mat nach Gl. (3) dieses » nicht dureh die veränderliehen 
Koordinaten .r, y,z auszudrücken ist, sondern dureh die Koordinaten «,,2,, die der 
Oberflächenpunkt zur Zeit ?==0 hatte. 


Der erste Ansatz für P (a, 9,2, t) ist 


pP, =AC0s uxe"yeosIe—"zsindcosrt. 


u ist so zu wählen, daß auf die Kanalbreite eine ganze Anzahl von Halbwellen kommt. 
Die Koeffizienten von y und z, u cos d und u sin d, sind so gewählt, daß die Differentialgleichung (1) 
erfüllt ist. Durch « ist im wesentlichen die Höhe der stehenden Wellen an der Stelle 2==0 
bestimmt, durch 9 der potentielle Abfall in z-Riehtung. Fürs erste sind beide unabhängig 
voneinander. 


Die Randbedingung für y h ist praktisch erfüllt, wenn 4 sehr groß Ist, was an- 
genommen werde (genauerer Ansatz durch Evo] u (y —h) cos"). 


Der Ansatz P =, kann die obere Randbedingung, Gl. (2) und (3) nur in erster Näherung 
befriedigen, d.h. y, wird gleich y,, wenn man nur die Glieder mit @ in erster Potenz berück- 
’ h e 1 | u cos u) 
sichtiet. Hierbei erhält man die Bedingungen y,- 0 und - = 
| ir 
Um die obere Randbedingung bis zu Gliedern mit a? einschließlich zu befriedigen, muß 
für PD außer y, noch eine zweite Potentialfunktion y, hinzugenommen werden, die von a’ ab- 
hängt. Man erhält 
a” u“ ke . ! 
(1,= L cos (2 u ysind)em?"?>sinssin2rt. 
2 ” u 
Bleibt man bei dieser Näherung, so kann man jetzt die Bewegung der Punkte berechnen, 
die für ?=0 in der (Juerwand 2-0 liegen. 


1), M. Schuler: Der Umschlag von Oberflächenwellen. Diese Zeitschrift 1035, Heft 6, S. 443. 
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Ihre Verschiebung in z-Riehtung wird 


atuısın d en Eh u. | 
\ | e" vsnscos ucsınv + | , sın’de2"yc0os9cos2 uxrcos2rt 
v R Y l 
a? u? 


— sın d cos? de?" veosd9 cos? yFt 
„s 


+- (Glieder ohne Ef). 

Wichtig ist das dritte Glied, das unabhängig von « ist, aber die doppelte Frequenz von 
y, hat. Ein Vergleich «dieses Gliedes mit dem zeitlichen Verlauf der Oberflächenverwölbung 
zeigt folgendes: 

Ist die Oberfläche eben, so sind die Punkte der Querwand zurückgewichen, ist die Ober- 
äche am stärksten gewellt, so sind sie vorgerückt. 

Das zweite Glied von £ kann dureh Überlagerung von Gliedern, wie sie im nächsten Vor- 
trage „Kanalwellen mit geringer Wellenhöhe* gebracht sind, beseitigt werden. Dieses gelingt 
aber nieht mit dem ersten Gliede. Um dieses zu beseitigen, ist ® bis zur dritten Näherung 
zu ermitteln, dann die Verschiebung [ der Querwand zu finden. Das Glied mit cos u x sin» t 
ist dann mit einem zweigliedrigen Ausdrucke behaftet, dessen erstes Glied wie vorher a in 
erster Potenz und dessen zweites Glied a’ enthält. Damit auch dieses Glied durch Über- 
lagerung weiterer Ausdrücke beseitigt werden kann, muß eine Beziehung zwischen « und 9 


| 


bestehen, die ermittelt werden kann. I.) 


Kanalwellen mit geringer Wellenhöhe. 
Von €. Weber in Dresden. 


(segeben sei ein Kanal mit rechteekigem (Juerscehnitt, Breite db, Höhe h. Ein recht- 
winkliges Koordinatensvstem ist so gelegt, daß der Koordinatenanfangspunkt auf der hinteren 
unteren Kante liegt. die „„-Achse quer zum Kanal, die „Achse nach oben und die z-Achse 
in Längsrichtung des Kanales geht. An der Stelle z2=0 befindet sich eine feste (Querwand, 
an der Stelle z=! eine bewegliche Querwand: die Geschwindigkeiten w in z-Riehtung der 
einzelnen Punkte dieser Wand sind gegeben. 

(Gesucht wird die Potentialbewegeunge einer idealen Flüssigkeit im Kanale, unter der 
Annahme, daß die Bewegungen klein sind, so daß die Bewegungsgleichungen ohne quadratische 
Glieder gelten. 

Die Geschwindigkeit dvd mit den Komponenten u, ©, w wird: 


— Br | Pol 0” \2 = 
'E oral Pla, z, PM), I» -+.— +:—)]P=0, 
„ 022 "dd ze) 
Für eine periodische Teilbewegung ist 
D— $bir,y.z)sin»t (oder cosrvf), | 
aD. OD. OD. (l). 
H sinrt, v= —sinvrt, w __sinvt, AP=0 
Ö&x Y 02 
Die Randbedingungen sind: 
i ap 
Für #0 und ‚5b (Seitenwände des Kanales) | " Ve 3 
O.X 
j ap 
für „=0 (Boden des Kanales) | N 
| y 
\D 2 
( ) 
für h (Oberfläche) = B=0 .,. Va CE e 
ey 4 
a» P 
für 20 (feste Querwand) E= Be Er a ee ie. Mi 
op 
Lan - - - .. mn r 6) 
für 21 (bewegte (uerwand) er Fi, (6). 
Fir, »p) ist gegeben. 
Untersucht werden drei Fälle: 
a) Die Querwand z = I steht still, Fir, 0; 
bi die Querwand 2-1 vollführt periodische Bewegungen, gesucht die periodische 


L,ösune: 
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ec) die Querwand z =]! steht für t<0 still, die Flüssigkeit ist hierbei in Ruhe; von 
t=0 bis t==#, sei eine beliebige Bewegung der Querwand gegeben. 
Fallal. Fie,=0. Für die Teilbewegung gilt der Ansatz 


m na une . 
Psinvt=cos x. cos BEER E : 5 “2 a 8 5 8 


b | 


BL E..,. 
Die Randbedinzungen sind bis auf (4) erfüllt. 


Differentialgleichung (1) gibt 


7 
. 


2 be = ” "| 
2? os (9). 
h b) 
Die Kandbedingung (4) gibt 
Pas äauh . » ..» . De an (9). 


Ist für *=0 an der Oberfläche D=0 (und wegen der Potentialbewegung überall v = 0), 
so wird das Potential der Gesamtbewerung 


Pi \ ’ m ı not 2 

D- Amm COS, 0 C08 7 2-60] Zu, y SI rn (40). 
N 2 
m, 


Die a,,,„ Sind aus der Oberflächenverwölbung für ?=0 zu bestimmen. Ist für 1-0 
jedoch D Z0, so treten noch Glieder mit cos» t hinzu. 

Fallb). Die Bewegung in z-Riehtung der Querwand zT ist w = Fler, y)sinrt. 

Das » ist jetzt gegeben. Wir machen für eine Teilbewegung den Ansatz 


m. 


EEE 
, UA RUM yE 0: 2 . tl) 


Pix, y,z)sinvt= cos 


Aus der Randbedingung (4) folgt für x 


‚2 ‚7p7 ür reelle 2 . . er 6 
4 | z | für imaginäres z., 
y’ ) x x \ | ur (l2b). 
tel h) also reelles |) | 


Aus der graphischen Lösung folgt ein reeller Wert z, und unendlich viele imaginäre 
Werte x,, 8: 


Aus der Differentialgleichung (1) erhält man für jeden #-Wert einen dazugehörigen /-Wert 
. = ME 


Jede Teillösung gibt für z=1 eine Bewegung mit der Geschwindigkeit 


ma i , 21: 
Im — C}, COS 7 2 prly) sn vt miıt 4 ol) Go] %Y; Fun (y) = «os ; Y, l; un .. 
Diese Funktionen y;, k=0,1,2... sind Lösungen einer Integralgleichung mit sym- 


metrischem Kern. Hieraus läßt sich zeigen, daß jede Funktion Fix, ,y) durch die Einzel- 
lösungen angenähert werden kann. 


Falle). Bis 20 sind Querwand und Flüssigkeit in Ruhe; von f--0 an bis t=14, 
ist die Bewegung der (Juerwand gegeben, wobei für 0 die Geschwindigkeit noch - 0 sein 
muß wegen des Anschlusses an die Zeit vorher. 


Die Bewegung von f=-0 bis t=f, wird beliebig so ergänzt, daß eine periodische 
Bewegung entsteht. Hierfür ist die periodische Lösung nach Fall b) zu finden. Für #--0 
erhält man jedoch eine Oberflächenverwölbung, die dureh Überlagerung einer Lösung nach 


Fall a) zu beseitigen ist. 296 
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Widerstandsuntersuchungen an Schiffen. 


Von Georg Weinblum in Berlin. 


Problemstellung. Die früher rein empiriseh-experimentell gerichteten Bemühungen, 
eünstige Schwimmkörperformen zu entwickeln, haben durch Michell, Havelock, Hogner 
und H. Warner einen starken Auftrieb von seiten der Theorie erhalten. Ausgehend von 
einer Arbeit in dieser Zeitschrift *) wird an einer Reihe kennzeichnender Beispiele ein weiterer 
experimenteller Nachweis der Tragfähigkeit der Theorie erbracht: alleemeinere Betrachtungen 
über einige für die Schiffstheorie grundlegende Minimalprobleme schließen sich an. — Bei 
der Behandlung der Schiffsformen werden die Hauptabmessungen L BT hier als konstant an- 
eesehen mit Ausnahme der Fälle, in denen es sich um eine Diskussion endlicher und unend- 
licher Tiefgänge handelt. Der Einfluß der Hauptabmessungen auf den Widerstand ist Gegen- 
stand besonderer Untersuechnngen, über die ich später berichten will. 


I. Modellversuche. 


ich habe gezeigt‘), wie man mit einfachen Polynomansätzen charakteristische Ände- 
rungen der Schiffsform darstellen und die entsprechenden Variationen des Wellenwiderstands 
erfassen kann. Zur Prüfung der Ergebnisse sind in der Preußischen Versuchsanstalt für 
Wasserbau und Sehiffbau 3 neue Versuchsserien durchgeführt worden. 


1. Scharfe Schiffe’. Wir eehen von zwei Auseanesmodellen der Grundform 
1 -a,(®—E E-9)-2]y(J)=-CWLR& 


mit «ddem Schärfeerad « 9/0556 und den Eiıntrittswerten —=?t—=( und 1 aus: wir 


! 


+») 


ändern Uen Schärfeerad der beiden Modelle, die sich dureh die Zusatzfunktion < 23,4: 0: 1> 
unterscheiden, dureh die Nlächenändernde Zusatzfunktion <2,3,4: 0,.1:0>. Die Versuche 
zeieen, daß Vorzeichen und Größenordnung der Änderungen des Wellenwiderstands auf Grund 
von Reehnung und Experiment gut übereinstimmen, sobald überhaupt eine nennenswerte 
Wellenbildung auftritt: ebenso erhält man im wesentlichen sich deekende Angaben über die 
Grenzen, bis zu denen eine Verringerung des Schärfegrades mit einer Widerstandsverringerung 
verbunden ist. Dureh Kombination eines Vorsehiffs mit verschiedenen Hinterschiffen kann 
auf den Einfluß der Zähiekeit geschlossen werden, insbesondere wieweit das Vorschiff für 
den Wellenwiderstand ausschlaggebend ist. Wie früher festgestellt. kommit bei hohen Froude- 
schen Zahlen (starker Wellenbildung) der Einfluß des Hinterschiffs immer mehr zur Geltung. 

Schlüsse über die praktische Verwendungsfähigkeit des Reduktionsfaktors von Havelock ') 
lassen sıch ziehen, 


2. Schiffe mit Endwulsten. Der Geltungsbereich der Theorie ist nicht so zut gewahrt 
wie im ersten Falle. Zwei Modellformen <2,4,6: 0,6> mit t=0 und f=1 werden ohne 
und mit einer Wulstform Iy, als" "ya, (£)*) untersucht; im Gegensatz zu bisherigen 
praktischen Ausführunren ordnen wir die Wulste zunächst sowohl am Bug wie am Heck An 
und danach erst nur am Bug oder Heck. Die „Tiefenverteilung* des Wulstes wird einmal 
proportional dem Hauptspant y, (£) = y(S) und sodann entsprechend der üblichen Ausführung 
zum Kiel hin verstärkt („versenkter Wulst*) , (£) == y (£) angesetzt. Ich nenne einige Ergeb- 
nisse der Rechnung und des Versuchs: a) y, ({) = yr(£). Bei höheren Froudeschen Zahlen ist 
die zum Hauptspant symmetrische Wulstanordnung günstig: der Bugwulst allein hat ın diesen 
Bereichen nur ungefähr die halbe Wirkung: Vorzeichen und Größenordnung des Wulsteinflusses 
auf den Wellenwiderstand wird von Froudeschen Zahlen $ > x 0,22 an von der Theorie richtig 
wiederrereben. b) Beim versenkten Wulst bestätigt sich die Rechnung ebenfalls: er ist 
unter den zemacehten Voraussetzungen schon vom 3. Buckel an (75 + 0.22) weniger wirksam : 
der Heekwulst scheint dureh Zähigkeitseinflüsse stärker paralysiert zu sein als unter a), so 
daß die Vorzüge der symmetrischen Anordnung erst von $ >05 zur Geltung kommen. 
Zusammenfassend: Die Theorie gibt nicht nur eine befriedigende Erklärung der Vorgänge, sie 
läßt sich auch als Hilfsmittel bei der praktischen Konstruktion anwenden. 


ı ZAMM 1935, S. 200. 
Erscheint in der Zeitschrift „Schiffbau“. 
Havelock, Proe. Royal Soe. 1935. 


ZAMM 1985, S. 2/5, Gl. (30). 














angew. 
I Mech. 


ıngen, 
ener 
l von 
ıterer 
uneen 

Bei 
nt an- 
ınend- 


reLeN- 


Ände- 
tands 
lt für 


wir 


: 1> 
suche 
rund 
werte 
r die 
erung 
kann 
f für 
ude- 
tung. 


ock *) 


wahrt 
ohne 
rieen 
»k an 
ıinmal 
rung 
‚reeh- 
n ıst 
lesen 
USSEs 
ichtig 
er Ist 
sam; 


1), SO 











Band 15, Heft 6 


Dezember 1935 








3. Völlige Schiffe. Nach Gesichtspunkten wie unter a) wurden Varianten für den Schärfe- 
orad und den Eintrittswert ? vorgenommen. Die Versuche zeigen, daß für >» 0,21 die 
Voraussagen der Theorie durchweg zutreffen und zu praktisch bedeutungsvollen Ergebnissen 
führen: der hohe Sehärfeerad bildet also keinen Hinderungserund für ihre Anwendung. Da- 
eeeen sind die Erzrebnisse unterhalb des dritten Buckels der Widerstandskurve (8 < » 0,20) 
zum Teil nieht eanz klar: eine besondere Untersuchung soll hierüber angestellt werden, da 
eine Lösune dieses Problems Wesentliches zur Konstruktion rationeller Frachtschifformen 
beitragen kann. 


II. Einige Minimalprobleme der Schiffstheorie. 

Für den Wellenwiderstand eines „Michell* schen Schiffes. eines getauchten Rotations- 
körpers (Havelock) und eines Drucksystems lieren Integraldarstellunzen vor, ebenso für den 
Reibungswiderstand eines Drehkörpers (Millikan, Kärmän-Moore). Es liegt nahe, hier- 
aus Formen eerinesten Widerstands abzuleiten. Isolierte Ansätze für den Wellenwiderstand 
haben nur Aussicht auf teehnische Bedeutung, solange die Voraussetzungen der Theorien 
einieermaßen eineehalten werden, die Zähiekeit keine entscheidenden Veränderungen hervor- 
ruft und die Formen praktisch zu verwirklichen sind. 

I. Unter den genannten Voraussetzungen habe ich aus dem Michellschen Integral A nach 
der Methode von Ritz günstige „Normalschiffe* abgeleitet’) (Wasserlinien mit nieht mehr 
als zwei Wendepunkten zwisehen Bug und Hauptspant, eine obere Grenze für F und eine 
untere für die Krümmunesradien). Die generelle Übereinstimmung mit Resultaten des 
Experiments war überraschend gut. Von Karman und Pavlenko haben gezeigt, daß die 
exakte Methode der Lösung des Variationsproblems für T-- und konstantes Volumen V in 
weiten Bereichen von % zu unsinnigen Formen führen kann: für einige 75 hat Pavlenko brauch- 
bare Lösungen angegeben. Im Verlauf einer Diskussion hat Herr Professor v. Karmäan einige 
Feststellungen gemacht und zwei Bedenken formaler Art gegen meine Formulierungen erhoben, 
die ich berühren will. Der entscheidende Passus in dem zusammenfassenden Berieht®) und 
in einer Zusehrift”) Herrn v. Kärmäns über die Konvergenzfrage läßt sich ım Hinblick auf 
die genannten Ergebnisse Pavlenkos ohne Einschränkungen nicht aufreehterhalten; eine Aus- 
sage über die Verhältnisse in dem von mir behandelten Fall Tendliech gestatten die bis jetzt 
miteeteilten Erzebnisse Herrn v. Karmäans nicht. Ich wıll aber zu meinen Ungunsten annehmen, 
daß auch hier ähnliche Schwierigkeiten auftreten. Eine bedinzungslose Identifizierung von 
Minimallösunzen des Michellschen Integrals mit Formen geringsten Widerstands muß zu Wider- 
sprüchen führen. Hinsichtlich der Terminologie „Ritzsches Verfahren“ würde ich Herrn 
Professor v. Kärmän beipflichten, wenn nicht in meiner Arbeit’) trotz aller einschneidenden 
Vorbehalte ein Variationsanspruch, Normalformen geringsten Widerstands zu approximieren, 
enthalten wäre: die geeienete Bezeichnung für die von Herrn v. Kärmän skizzierte Aufzabe’) 
wäre nach einem Vorschlag von Trefftz „Methode der optimalen Polynomlösungen*, die ich 
de faeto weiteehend verwende'!). leh bin sicher, daß Einverständnis über die Beibehaltung der 
Bezeichnung „Ritzsche Methode“ zu erzielen ist; letztere Methode behält ihren praktischen 
Wert. solange die ersten Glieder des Ansatzes vernünftige Lösungen ergeben, die der 
Versuch als günstige bestätigt, was bis jetzt der Fall gewesen ist, wobei es gleichgültig 
bleibt. ob die exakte Lösung brauchbar oder unbrauchbar ıst. Ich habe diesen Gedanken 
durch die Behauptung auszudrücken versucht, daß die Ritzsche Methode unter Zugrunde- 
legung von Polynomansätzen das Minimalproblem im Rahmen der Michellschen "Theorie 
für die Anwenduneen löst: der Einwand Herrn v. Kärmäns gegen diese zugegebener- 
weise mieht strenze Formulierung erscheint mir nicht ganz zutreffend’): geht man um ein 
Glied mit dem Parameter a,„,+, weiter, so wird in dem alleın interessierenden Falle einer 
fehlenden Konvergenz die neue Optimalform nach v. Kärmäan nur wenig von der früheren 
verschieden. ihr Widerstand aber etwas geringer sein: m. E. braucht aber die Änderung der 
Form nicht klein zu bleiben. der Geltungsbereich der Theorie wird mehr überschritten als 
[rüher, und man kann nur aussagen, daß der Integralwert R,+,<R,„ wird, der tatsächliche 
Widerstand (selbst in idealer Flüssigkeit) W,+, kann ungünstiger als W, werden.  Eın 
Kriterium, wann die Entwieklung mit Rücksicht auf die praktische Anwendbarkeit der 
Ergebnisse abgebrochen werden muß, gibt bis jetzt im wesentlichen nur der Versuch, wenn 
auch einiee Absehätzungen über die Fehlergerenzen des Michellschen Integrals von Havelock °) 
mitgeteilt sınd. 


Für eine einparametrige Schar läßt sich zeigen, daß der Minimalansatz bei üblichen T 


und höheren /$ noch zu vernünftigen Formen führt, wenn unter sonst gleichen Bedingungen 


5), ‚Schiffe geringsten Widerstands“, 3. Int. Kongr. Meeh., Stockholm 1939, 
6, v. Kärmän: IV. Int. Kongreß f. Techn. Mech, Cambridge 1954. 

”, v. Kärmän: Zusehrift Werft, Reederei, Hafen 1935, Heft 8. 

$, Havelock, Proc. Royal Soc. 1923, S. 571 und 1931, S. 40. 
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bei T=» sehon unbrauehbare Linien entstehen. Für ein Polynom mit endlicher Gliederzahl 
lassen sieh in erster Näherune in den Grenzfällen K>x und $ >00 einfache Zusammenhänge 
zwischen Widerstand, Völliekeit und Eintrittstanzente angeben. 


Sollten wider Erwarten im Bereich der praktisch wichtigen Froude-Zahlen % schon die 
ersten Glieder eines Ansatzes nach Ritz zu unsinnizeen Ergebnissen führen, so kann man 
versuchen, den Ansatz S(R-+-(ÜO0)=0 zu machen (C eine Konstante, O0 Oberfläche), der 
praktisch sinnvoll ıst (ein Minimalansatz für den Gesamtwiderstand scheitert an der Un- 
kenntnis der anderen „Widerstandsarten*). Richtiger ist es jedoch, in solchen Fällen die 

h 07; 
unzenaue Randbedineunge Michells # x e \ _ zu verbessern. Dahinzehende Bemühungen haben 
vorläufie gezeigt, daß die Abweichungen des so berechneten Inteerals von R bei üblichen 
einfachen Sehiffslinien nicht bedeutend sind; dagegen erhält man bei hohen % für ein- 
parametrige Scharen Optimalformen, die besser mit den auf experimentellem Wege abgeleiteten 
übereinzustimmen scheinen. 


Diese Ausführungen besaren nichts gegen den Wert der in Frage stehenden Uhnter- 
suchung Professor v. Karmäns, die eine willkommene Ergänzung der experimentellen Verfahren 
bedeutet; ebensowenig kann man die Arbeit v. Kärmäns gegen die von mir gewonnenen 
experimentell gesicherten praktischen Ergebnisse ins Treffen führen, wie das von dritter 
Seite zeschehen ist: letzteren Versuch betrachte ich hiermit als erlediet. Meine früher?) 
ausgesprochene Schlußfolgerung, daß die Theorie Michells geeignet erscheint, günstige 
Schiffslinien abzuleiten, bleibt in vollem Umfange bestehen, wobei, wie sehon dort hervor- 
vzehoben, die Erreichung des absoluten Optimums der Ausführung natürlich gänzlich außerhalb 
der Reichweite der bestehenden Theorie bleiben muß. 


>. Ähnliche Ansätze lassen sich zur Berechnung günstiger endlich tief getauchter 
Rotationskörper durchführen’), desgleichen zur Berechnung günstiger Drucksysteme an der 
Oberfläche. Bemühungen, die Lösungen von Millikan und v. Käarmän-Moore für die Be: 
stimmung von Drehkörpern geringsten Reibungswiderstands in entsprechender Weise aus- 


| 


zunutzen, sind zunächst mit großer Rechenarbeit verbunden. 557 


Angenäherte Bestimmung der Druckverteilung an einem Flügelprofil. 


Von F. Weinig ın Berlin. 


Bei nieht zu dieken und nicht zu sehr gewölbten Profilen genügt es, die Randbedingungen 
auf der Profilsehne statt auf der Profilberandung selbst zu erfüllen. Hierbei benutzt man als 
Profilsehne die Sehne der Profilmittellinie. Als Zuordnung der Profilberandungspunkte zu 
Sehnenpunkten empfiehlt sieh zunächst ihre Projektion auf die Mittellinie und deren nach- 
herize Projektion auf die Sehne. 


I 


Ist die Strömung um das Profil in der Ebene 765)=g-+ iy, so Ist 


Be t T 
iln 2 —_y—iln 
dz Wo 
wobeı y der Winkel der Strömunesriehtung eeeenüber (ler - ‚r-Achse, Ih ler Absolutwert der 
(weschwindiekeit ist. 


Bildet man die Sehne, deren Endpunkt 3= —-2 und 3= +2 sein mögen, auf den 


Kinheitskreis Ö = e’“ der ©-Ebene ab. so ıst 


und man kann schreiben 
yı-) 


9, „Rotationskörper reringesten Wellenw iderstands.* Vortrag, eehalten vor der (1es, ag. Math. ;erlin, 52: Il. 19:5 
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Dann ist 


v—ö+(», cosa- u, sin a) +(r, cos 2a 4 Hs sin a) 3 r 
und 
m | | 


In : (v, sina — u, cosa) + (r, sin 2a — u, cos 2a) 


NV oo N . p” 


= m 
Da v (a) gegeben ist, erhält man durch Fourieranalyse die v;, «; und kann damit In 2 berechnen. 
LH 


d ist hierbei der Anstellwinkel der Profilsehne gegen die Strömungsrichtung in unendlichen. 


Als Auftrieb- und Momentbeiwert bezogen auf 3-0 findet man mit Hilfe der Blasius- 
schen Formeln 


f . ’m | '» . 


«, muß identisch verschwinden, wenn die gemachten Vernachlässigungen nicht zu groß sind. 


Zwecekmäßig bestimmt man zuerst die Geschwindigkeitsverteilung bei glatter Umströmung 
(Index 0), d. h. für jene Anstellung, bei der der vordere Endpunkt der Mittellinie Ver- 
zweigungspunkt der Strömung ist. Man findet dann 


Im yn yn 
- | ’ ’ —— | ’ » \ ’ | ’ * *) 
DR 5 \r.da, =, \mcosada, 7 \meos2ada, 
u u u 


Weiter wiehtie ist der Einfluß eines Querruderausschlaes 5 (Index p). Das Profil möge 

v 5» | | 
also nunmehr an einer Stelle um den Winkel 5 abgeknickt sein. Die Rudertiefe sei f. Dann 
entspricht der Knickstelle ein Bildkreiswinkel 


4 
€ = 27% cos | | E 
| I 
Man findet 
[7 ) ) 
) ) n > £ < 
O3 E E. v3 pr ı SINT, v3 sın Dr. 

7 zı | Bu 


Bei der Untersuchung einer beliebigen Anstellung d,,. gegenüber der glatten Umströmung 
(Index zus) hat man zu beachten, daß der Bildkreiswinkel des vorderen Verzweigungspunktes 
um 20,,. vergrößert ist gegenüber dem beı glatter Umströmung. Man findet dann 


2) Ö Yu, + Bid, .;, RR „sin do 


zus ® 
Das Resultat für beliebige Anstellung und beliebigen Ruderausschlag findet man aus 


| | | ’ > ’ ’ ’ ’ > ı ’ 
n) 0, +03 t u; v,: Yo trıB trau: N, "oTtro? TrYazus: 


Die Geschwindigkeitsverteilung bei glatter Umströmung kann man statt durch die angegebene 
keihenentwicklung auch direkt erhalten aus 


In 
m, ‚a—äÖ 
In = (a) = \ ‚‚dinsin 
I v -. 
A=% 


Es empfiehlt sieh für die Auswertung dieses Inteerals ein Auszählverfahren : 
| u g 


Dr Be h MM 
Man legt ein Netz mit einer linearen »-Ordinaten- und einer In sin „ -Abszissenteilung 


mit dem Abszissenanfang in a auf die gezeichnete » (a)-Kurve und zählt die die » (a)-Kurve 
überdeckenden Rechtecke aus, bzw. schätzt deren Zahl ab. 


Für beliebige Anstellung und beliebigen Ruderausschlag ergibt sich 


our ac Alf, , BP, tsaörter/2f a 
= ’ S 4 98 \ l oO () oO 2 
we | ‚7 \ vd Insın > ! = In tg a)® — tg zjaf \60°° et tg 0,1.t18 TE Eng 
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CONSTANTIN CARATHEODORY, Prof. an (ler 
Universität München. Variationsrecehnunge 
Ind partielle Different ialeleichu neen 
Ordnune XI 107 8. m. 31 Fig. im 
leipzig und Berlin 1935. B. G. Teubner Ver 
s 22 M. 


erster 


ext. 
la®. Prei 
) insofern eine 
ersten Male die 
Beziehung (der 
Theorie der partiellen 
Vordergrund rückt 
wesentlich 


uch nimmt 
zum 
erwähnte enee 


Dieses | besondere 


Stellunz ein. als es schon self 
Jacobi wiederholt 
Variationsreehnunge zu (der 
Differentialeleichungen in 

und deshalb die Variationsreehnune in 
anderer Weise aufbaut. als es üblich ist. 
zerfällt in zwei 


(eeenstandes 


den 


trotz «(des 
zusammen 


Teile, die 


oreanisech 


Das Buch 


\ erschiedenen 


hängen: acht Kapitel über partielle Differential 
elejehuneen und ebenso viele über Variationsrech 
nune, Nur sind (diese etwas läneer. Zu Anfan«e 


stehen noch zwei Kapitel einleitender Art, in denen 
einige Teile der Analysis behandelt werden, die der 
Verf, braueht. die aber noch nieht alleemein zu dem 

\nfänzervorlesune zehören. >o 
etwas üben Konvereenz. normale Familien. 
eewöhnliche Differentinleleichungen, Felder von 
Kurven und die vollständigen Systeme nach Adolph 
\Maver, womit wir mitten in «der sind. 
lritte Kapitel behandelt dann die partiellen 
Differentialeleiehuneen erster Ordnung und die Cha 
wird weiter 
stark formaler Art 
Im vierten Kapitel 
Larrangeschen 


von 


einer 


stetioe 


Geerenstäanıd 


sehon Sache 


Das 
rakteristikentheorie. \uf «diese alles 
aufrebaut. die älteren 
im zehnten Kapitel. 
l’oissonsehen und die 
Klammern mit Anwendung auf (die 
Differentialeleiehungen. Es kommt ein Kapitel über 
Tensorreehnunge, «das auch organisch zum Ganzen 
| nichts anderes als die 
Transformationstheorie 


Fheorien 
] y , . % 
kommel ers! 


stehen die 
Systeme 


da es ja eigentlich 


Theorie S(0) 


eehört. 
für die wichtige 
bedeutet. es folgen dementsprechend auch Kapitel 
über kanonische Transformationen und Berührungs 
transformationen. Ein überaus klares Kapitel über 
las Pfaffsche Problem sehließt sieh an und vor (dem 
zehnten kommen dann noch die Funktionengruppen. 

Die Variationsreehnung wird durch ein Kapitel 


über gewöhnliche Maxima und Minima und quadra 


tische Formen eineeleitet. Es ist die wesentliche 
Stütze «(des foleenden. Denn bei dem Ausbau der 


Variationsreehnung erscheint vermöge einer sSpe- 
iellen Variation das Positivsein einer quadratischen 
Form als ausschlaeeebend. Von hier aus geht der 
Wer Hamiltonschen partiellen Differen 
tialeleiehun®e und zur Weierstraßschen E-Funktion. 
Das ist wohl das Bezeichnende Buches. 
Ferner. (dab die Eulersche Gleiehung zuerst als hin 
und zum als notwendige De 
lineune erscheint. 
Dieser Wer, der weder 
Standpunkt erste ist, hat 
Krstens können die 
und zweitens führt der 
Kapitel erstaunlich klar 
hei der Theorie L.a_eranzeschen 
hei Nebenbedinzunzen 


weiter zur 


(lieses 


reichende Schluß erst 
historisch noch vom naiven 
aus (den zwei Vorteile. 
Voraussetzungen etwas ermäßigt 
Were später, im 
und durehsichtig 
Problems. also 


werden 
jetzten 
des 
zum Ziele, 
Nachdem im dreizehnten Kapitel (die 
P’arameterdarstellung analog behandelt worden sind. 
£ im vierzehnten die positiv definiten Probleme. 


l’robleme in 


leen 
Das fünfzehnte bringt «die quadratischen Probleme 
ınd im Zusammenhang «damit die Theorie der 
weiten Variation. Dabei «die Sturm-Liouvilleschen 
Oszillationstheoreme und die Theorie der konjugierten 
Punkte. Kapitel XVI handelt von dem Ranıdwert 
problem und der Frage nach dem absoluten Mini 
num. Hier gibt es Zahl interessanter 
Beispiele zu merkwürdigen Vorkommnissen aller 
Art. Auch das Randwertproblem im großen ist nach 


. ) 
eine erobe 
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ausführlich behandelt. Die Sätze 8 und 9 
Kapitels mögen besonders hervorgehoben 
werden. Das siebzehnte Kapitel enthält noch die 
Theorie der geschlossenen Extremalen und der peri 
odischen Probleme mit Sätzen von Poincare. Hada- 
mard und Razmadze. 

Im ganzen ein starkes, konzentriert geschriebenes 
Buch, das noch als Lehrbuch bezeichnet werden 
kann, aber als kein zanz leichtes, Mit Absicht 
wurden die weometrischen Theorien Lies etwas 
zurückgestellt. da im selben Verlag ja schon das 
Buch von Engel und Faber da ist. Die formalen 
Methoden Jacobis und Cauchys überwiegen gegen 
das Geometrisehe, dessen Zurücktreten man im In 
Anfängers an einieen Stellen bedauern 
kann: z. B. bei der Charakteristikentheorie. deren 
Kinführune dem Anfänger als willkürlich erscheinen 
muß. Auch Seite 26 und ff. wäre etwas mehr an 
schauliche Behandlungsweise sicher manchem An 
fänzer anzenehm. Das Weerlassen der Summen- 
zeichen liebt Ref. gar nicht. es zibt den Formeln 
den Anschein einer falschen Leicehtiekeit. Verf, ist 
ersichtlich anderer Meinune, 

\ber das ist unwesentlich gegen den prachtvollen 
Inhalt. von dem neu Ist, anderes bis zur 


Hilbert 


lieses 


tPITESSP (les 


manches 


heutieen Forschune heranreicht. Die Theorie von 
lonelli wird z. B. weitgehend hineingearbeitet. Ein 


ausführliches Literaturverzeichnis von 190 Nummern 
ermöelieht die weitere Vertiefune. Sehr dankens 
wert sind die Ratschläge, die zur Benutzung der 
Literatur gegeben werden. Die Darstellung ist klar 


und selbstverständlieh einwandfrei. die Voraus 
setzunzen sind genau anzereben. Im ganzen wird 
auch eine erobe Vereinfaechune «dureh den neuen 


Aufbau erreicht: hat der Leser gewisse Schwierig- 


keiten überwunden. so wird er dureh die reichen 
Kreebnisse helohnt Dinze. (die SONST besonders 


schwer sind. werden hier verhältnismäßie leicht. 
Vor allem mare (den Lesern dieser Zeitschrift noch 
vesaet sein. daß die Beziehuneen (des behandelten 


Gerenstandes zur Mechanik überall durehleuehten. 
z. T. in Anwenduneen direkt zutawe treten. Als 
küstzeug zu einer vertieften Behandlung schwieriger 


Probleme aus der Mechanik wird das vorliegende 


Buch bald unentbehrlich sein. Hamel. 512 
HANS BOLZA, Ein neuer Wee zur Kr- 
[forschung und Darstellung volks- 


irtschaftlieher Vorränge. vier 8 


m. 68 Abb. Berlin 1935. Verlag Julius Springer. 
Preis kart. 4.80 M. 
Der Verfasser wirbt mit seinem Buche für die 


Kinführung exakter mathematisech-naturwissen 
schaftlicher Methoden zur Betrachtung wirtschaft 
licher Voreänge. Darüber hinaus will er das In- 
teresse der Mathematiker, Naturwissenschaftler uni 
Inzenieure für wirtschaftliche Fraren gewinnen, in- 
lem er die Anwendbarkeit der diesen zeläufigen 
Methoden auf Wirtschaftsfragen aufzeiet. Zu Be 
einn des Buches wird zezeirt. dab zwei Begriffe für 
die Beschreibung wirtschaftlicher Voreänge erunid 
lerend sind: „Menzeenfunktionen” (z. B. die Summe 
aller seit einem Anfangszeitpunkte erzeugten Mengen 
einer Ware) und deren  Differentialquotienten. 
„Leistungsfunktionen“ (z. B. die pro Zeiteinheit er- 
zeuete Menze). Die zeiehnerische Darstellung die- 
ser Funktionen erklärt und an Beispielen er- 
läutert. Des definiert (der Verfasser die 
umlaufenden ireendweleher Art 
Münzen. Zahlungesver- 


wird 
weiteren 
Zahlunesmittel 
Banknoten. Giraleeld. 


sprechen usw,) als Äquivalent der Differenzen von 
Warenmengen. Ferner werden die Begriffe der 


Kehrwert der 
und .Preis 


(les 


„Tausehzeit (selles” d. ı. der 
Umlaufseesehwindiekeit). des „Preises“ 
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niveaus“ definiert. Die Anwendung dieser Begriffe 
wird unter Zuziehune der zeiehnerischen Dar- 
stellung an Beispielen erläutert. Die letzten Ab 
schnitte des Buches gehen über das im Titel ge 
kennzeichnete Gebiet hinaus. Sie befassen sich mit 
dem Wesen des Buchgeldes. Besonders hervorge 
hoben wird die Krisenabhäneirkeit des Buchreldes. 
welche die Instabilität «der Wirtschaftslaere in 
Krisenzeiten zur Folge hat. Im letzten Kapitel 
macht der Verfasser Vorschläre. durch gesetzliche 
Maßnahmen diese Krisenabhängeiekeit zu beseitiren. 
eine Frare,. deren Lösung der Verfasser als die 
vordrinzlicehste Aufgabe der  Wirtschaftsführung 
ansieht. 

Der wirtschaftlich wenizer bewanderte Mathema- 
tiker wird beim Lesen des Buches zelegentlieh 
Schwieriekeiten haben und wird hier und da Be 
denken nieht unterdrücken können. So erscheint 
z. B. die Definition der Tauschzeit (S. 22) nieht ganz 
frei von Willkür zu sein: auf S. 30 bleibt der Vor 
ange der „Schaffune von Zahlunesmitteln"“ unklar. 
wenn man sieh an die vorher zegebene Definition 
ler Zahlungsmittel hält. Jedenfalls beizuptlichten 
ist der vom Verfasser auch an anderer Stelle ze 
äußerten Ansicht. daß die wirtschaftlichen Fragen 
mehr Interesse bei (den mathematisch-naturwissen 
sehaftlich Geschulten finden sollten. als es heute 
noch der Fall ist. Deshalb sei auf (las Buch hin 
eewiesen. Trefftz. 559 


Taschenbuch für den Maschinenbau, herausge- 
eben von H. DUBBEL. 6. Aufl. In 2 Bin. 
Berlin 1935. Julius Springer. Preis: 22,50 M. 


Die neue Auflage des Taschenbuchs hat die kenn- 
zeiehnenden Eigenarten der letzten Auflage dureh 
aus bewahrt. Das eilt vor allem von der Einfach- 
heit der Darstellung, die es einem Ingenieur erlaubt. 
auch dort etwas rasch nachzuschlagen und zu ver- 
stehen, wo es sieh nieht um sein engeres Fachgebiet 
handelt. Meistens sucht der Ingenieur in einem 
Taschenbuch Antwort auf eine ganz klar umrissene 
Fraee (z. B. Zahlenwerte),. Das Buch ist um so 
brauchbarer, je schneller es eine mögliehst eindeu 
tire Antwort auf die Frage gibt. Schwierigere 
erundsätzliche Fragen wird man im allgemeinen 
nieht an Hand von einem Taschenbuch. sondern 
eines Sonderwerkes behandeln. Ein Taschenbuch 
muß deshalb am besten eine möglichst umfassende 
Sammlung von einfachen Rezepten sein und braucht 
sich durehaus nieht zu schämen. wenn es nicht 
„wissenschaftlich“ in dem so oft falsch verstandenen 
Sinne ist. | 

Das Taschenbuch von Dubbel hat die genannten 
Erfordernisse in glücklicher Form erfüllt und wird 
sich deshalb sieher auch an den Technischen Hoch- 
schulen noch weiter «durchsetzen. In der neuen 
Auflage sind verschiedene Abschnitte entweder 
völlie oder zum Teil umgearbeitet worden. Auber- 
dem ist der gesamte Inhalt auf den jetzigen Stand 
der Technik gebracht. Dabei bringt das Buch die 
"ormeln und Zahlenwerte in einer Form, die es dem 
Benutzer fast immer erlaubt. die gesuchte Größe 
schnell zu bereehnen, ohne erst mehrere Seiten vor- 
her dureharbeiten zu müssen. 

Zum Schluß sei noch auf einige kleinere Mängel 
hingewiesen. Ableitungen zehören nicht in ein 
Taschenbuch, selbst wenn es sich um Hauptsätze 
handelt. Derartige Dinge werden ausführlich in 
Sonderwerken behandelt und nur dort auch dureh- 
vearbeitet. In dem Taschenbuch könnte dureh 
Fortlassen dieser Ableitungen mancher Platz für 
wichtieere Dinge freigemacht werden. Ferner er- 
scheint es ratsam, wenn in dem Abschnitt „Dynamik 
flüssirer Körper“ sämtliche Angaben über >trö 
mungesverluste in der Form der Gleichung auf 
Seite 341, unten. gebracht. d. h. von der Geschwin- 
diekeitshöhe abhängige gemacht werden. Das gilt 


vor allem von den an sieh sehr wertvollen Tafeln 
auf Seite 33D und 344. Dadurch entsteht unbedingt 
auch für den Ingenieur eine Vereinfachune, der es 
zewohnt war, z. B, mit äquivalenten Rohrlängen zu 
reehnen. 


Dresden. K. Sörensen. 52 


Dr. S.IGUCHI, Professor an der Hokkaido Kaiser- 
lichen Universität zu Sapporo, Japan. Eine 
lösung für die Berechnung der bieg 
samen rechteekigen Platten. 56 S. mit 
13 Textabb. u. 3 Taf. Berlin 1933, Verlage von 
Julius Springer. Preis 5 M. 


Der Verfasser berechnet orthogonal anisotrope 
rechteckige Platten (mit ungleicher Biegesteifig 
keit in den zwei orthogonalen Richtungen) und 
nimmt als Ausgang die Entwieklungen von M. T. 
Huber. Ausführlich untersucht werden die folgen- 
den, für die Anwendung wichtiren Fälle: Die auf 
vier Seiten freiaufliegende Platte bei hyılrostati- 
scher, gleiehmäßiger und partieller Belastung: die 
latte, bei der drei Seiten frei aufliegen und die 
vierte ganz frei ist: die auf vier Seiten voll ein- 
zespannte Platte: die Platte auf elastischem Unter 
erund. Ferner wird noch auf den Einfluß der 
Torsionssteifigkeit auf Durehbiegune und Biezungs 
moment hingewiesen und die Nachreehnune von 
Versuchen an kreuzweise bewehrten Eisenbeton 
platten kurz besprochen. Für den Benützer der 
Schrift sei vermerkt, daß die Kreis funktionen 
mit groben Anfangesbuchstaben bezeiehnet sind. 

Wien. J. Ratzersdorfer. 


Dr. teehn. JOHN-ERIK EKSTRÖM, Studien 
über dünne Schalen von rotationssymme- 
trischer Form und Belastune mit konstanter oder 
veränderlicher Wandstärke.  Ingeniörsvetenkaps- 
akedemiens,. Handlingar Nr. 121. Stockholm 1933. 
Svenska Bokhanleleentralen A.-B. 205 8. 


Das vorliegende Buch ist eine zusammenfassenide 
Darstellung der Bereehnunge von dünnen Schalen 
von rotationssymmetrischer Form und Belastung. 
Bei der exakten Berechnung weht der Verfasser 
von den Gleichungen Meißners aus, mit denen 
das Problem auf zwei symmetrische Differential 
eleichungen geführt wird. und untersucht die 
lösunzen von Meißner und seinen Schülern für 
sphärische Schalen, Ringflächenschalen und Kegel 
schalen mit konstanter Wandstärke und Keeel- 
schalen mit linear veränderlieher Wandstärke, 
Dann besprieht der Verfasser die Näherungsver- 
fahren von Bauersfeld-Geekeler. von 
Ritz und die Verwendung von Differenzengleiehun 


gen und zeigt die Bedingungen, unter denen jede Me 


thode zweckmäbig ist. — Hieran schließt sich die 
numerische Berechnung einer Kuppel. die aus einem 
am Scheitel offenen Boden eines Wasserbehälters 
besteht und die seitlich an vertikale zylindrische 
Wände angeschlossen ist. Die genauen Spannungen 
unterscheiden sieh erheblich von den Werten nach 
der gebräuchlichen Membrantheorie. bei der die 
Momente außer acht bleiben. Die Ergebnisse sin 
hierbei von der Größe der Poissonschen Zahl m 


wenie beeinflußt. so daß m x gesetzt werden 
darf. Zum Schluß untersucht der Verfasser sphä- 


rische Schalen mit veränderlicher Wandstärke und 
führt eine Änderung des Stärkeverhältnisses ein. 
bei der das Problem lösbar wird. In praktischen 
Fällen läßt sich dieser Verlöuf der Wandstärken 
infolee der Wahl von Koeffizienten dem wirklichen 
\erlauf nahebringen. 

Das interessante Buch muß allen sehr empfohlen 
werden. die sieh mit der Bereehnunz und Konstruk 
tion von Schalen befassen. 

Wien. J. Ratzersdorfer. 
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PIERRE DE FERMATS Abhandlungen 
über Maxima und Minima (1629). Aus dem 
l.ateinischen übers. u. m. Anm. versehen von Max 
MILLER. (Ostwalds Klassiker der exakten Wissen- 
schaften, Nr. 50 S. m. 19 Textfie. Leipzig 
1034. Akademische Verlarsees, m.b.H. Preis kart 
3.20 M. 

Lanze vor Erscheinen der ersten Arbeiten von 
Leibniz und Newton über die Infinitesimalrech 
nune hat Fermat die in diesem Heft zusammenge- 
stellten Arbeiten, die seit 1629 erschienen sind, ge 
schrieben, Er wendet darin ein ohne Beweis mitzeteil- 
tesVerfahren an, das auf das Nullsetzen der ersten 
\bleitung der entsprechenden Funktion einer passen 
oewählten Strecke der betreffenden Firur hinaus- 
kommt. Nach Methode werden sowohl die 
Bestimmune von Extremwerten behandelt in den 
beiden letzten Arbeiten z. B. der Were eines £e- 
brochenen Liehtstrahles nach dem Prinzip 
kleinsten Widerstandes bzw. der kürzesten Zeit 
wie aueh Konstruktionen von Kurventaneenten 2e- 
Iunden. Die nach «dem lateinischen Text der Veu- 
vres de Fermat (Paris 1891) hergestellte Übersetzung 
liest sich sehr eut, Wünschenswert wäre nur eine 
kurze historische Einleitune über die Stellung 
dieser Arbeiten in der Entwieklung der Mathematik. 
Fr. Willers. 500 
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Freibere Sı. 
SCHARDIN, Berlin. Das 
Schlierenverfahren. 
fürseine Anwendung und 

\uswertung (Forsehungsheft 


Dr.-Ine. HUBERT 
Toeplersche 
((‚rundlaren 
quantitative 


367. Beil. zu .„.Forsehune auf dem Gebiete des Ingze- 
nieurwesens“, Ausg. B, Bd. 5. Juli/’Aug. 1934). 
32 >. 79 Abb.. 5 Zahlentaf. Berlin 1934. VDI 
Verlae G.m. b. H. Preis 5M. für VDI-Mitel. 4.50 M. 


Wilhrend man beim Toeplerschen Schlierenverfah- 
ren bisher mehr oder weniger auf qualitative Be- 
trachtungen angewiesen war, versuchen die vor 
lieeenden,. auf Anreeune von ©. Cranz auseeführten 
Untersuchungen zu quantitativen Ergebnissen zu ge 


laneen. In der Arbeit werden die theoretischen 
(Grundlasen des Verfahrens soweit behandelt. daß 
man auf Grund der abseleiteten Formeln «#inmal 
einen zweckentsprechenden Aufbau der Versuchs 


vornehmen kann, andererseits aber 
Fällen aus der Liehtablenkung 
auf den Zustand des Objektes schließen kann. Wei 
ter wird eine Anzahl von Anwendungeseebieten so 
weit theoretisch besprochen. wie es für die \nwen 
dune des Schlierenverfahrens nötie ist. 
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anordnuneen 
auch in einfachen 


Freibere Sa. 

Dr. ROBERT HAUSSNER, 0. ö. Prof. (ler Mathe 
matik a. d. Un. Jena. Analytische Geome 
trie der Ebene. ?2. Aufl. 164 S. m. 60 Fie. 
Sammlung Göschen Bd. 65. Berlin und Leipzig 1934. 
Verlaxe Walter de Gruyter u. Co. Preis: Geb. 1,62 M 

Das Bändehen, eine eute Einführung in die ana 
der Ebene. behandelt den Stoff 
I'mfanze,. wie ihn etwa eine elementar 
einführende Vorlesung bringen wird. In 
(den \bsehnitten werden in der Hauptsache 
schiefwinklire eartesische Koortlinaten benutzt, die 
hier als Parallelkoordinaten bezeiehnet werden, eine 
von Plücker eineeführte Bezeichnung, die heute 
nieht mehr zu empfehlen ist. da man unter Parallel 


Ivtische Geometrie 


in dem oe 
haltene, 


ersten 


13:35 


koordinaten neuerdines meist etwas anderes ver 
steht. In diesen Absehnitten werden Gerade und 
Kurven zweiter Ordnung einzehend behandelt. Die 


weiteren \hsehnitte brineen unter  Benutzun«e 


von Linienkoortdinaten die Behandlune der Kurven 
‚weiter Klasse. Im letzten Absehnitt werden dann 
schließlich  Ähnliehkeitspunkte und -achsen von 





Kreisen und das appollonische Berührungsproblem 
kurz besprochen. 
Freiberg (Sa. Fr. A. Willers. 523 


dr. ROBERT HAUSSNER, 0. ö. Prof. der Mathe- 
matik a. d. Un. Jena, Analytische Geome- 
trie des Raumes. 132 8. m. 36 Fig. Samm- 
lung Göschen Bd. 89. Berlin u, Leipzig 1935, Ver- 
lag Walter de Gruyter u. Co. Preis: Geb. 1,62 M. 
Nach Behandlung von Punkt. Ebene und Geraden 
In inhomogenen Punkt- und Ebenenkoordinaten 
unter stetiger Hervorhebung des Dualitätsprinzipes 
werden homogene Koorilinaten eingeführt und unter 
Benutzung dieser Koordinaten die Flächen zweiter 


Ordnung und Klasse eingehend behandelt. Das 
Bändehen, das eine recht brauchbare Einführung 
in die analytische Geometrie des Raumes gibt. 


schließt mit der Behandlung der Fokaleigenschaften 
der Flächen zweiter Ordnung und ihrer Erzeugung 
nach Mae Cullach. Für eine Neuauflage wäre die 
Beseitieung einiger störender Druckfehler und die 
Neuzeiehnunge einiger nicht einwandfreier Figuren 
wiinschenswert. 


Freiberg (Sa.) Fr. A. Willers. 523 
Ferner sind bei der Schriftleitung folgende 


bücher (ausführliche 


bleibt 


Prof. Dr. A. THUM, VDI, und Dr.-Ing. F. WUN- 
DERLICH, VDI, Dauerbiegefestigkeit 
von Konstruktionsteilen an Einspan- 
nunzeen, Nabensitzen und ähnlichen 
Kraftangriffstellen. (Mitteilungen der Ma- 


einrerangen Besprechung 


vorbehalten): 


terialprüfungsanstalt an der Techn. Hochschule 
Darmstadt. hrse, von Prof. Dr. A. Thum, VDI. 
Heft 5.) VI S2 S. m. 94 Abh. und 10 Zahlentafeln. 


Berlin 1934. VDI-Verlae G.m.b.H. Preis brosch. 


750 M 

Dr.-Ing. HANS H. BLEICH, Die Berechnung 
verankerter Hängebrücken. 101 S. m. 
17 Abb. Wien 1935, Verlag Julius Springer. Preis 
kart. 12 M. 


Dr.-Ing. HORST MÜLLER, Privatdozent a. d. 
Techn. Hochseh. Hannover. Führer durch die 
technische Mechanik, eine neuartige Über- 
sieht über ihre Grundlagen. Methoden und Ergeb- 


118 S. m. 166 
Preis 


nisse für Studium und Praxis. VIII + 
Textabb. Berlin 1935. Verlag Julius Springer. 
broseh. 850 M. 

0, ZARISKI, Alzebraie Surfaces. (Ergeb- 
nisse der Mathematik und ihrer Grenzgebiete. hrsg. 
von der Schriftleitung des ..Zentralblatt für Mathe- 
matik“. 3. Bd. 5.) 198 S. Berlin 1935. Verlag 
Julius Springer. Preis brosch. 22.75 M. 

Enceyklopädie der mathematischen Wissenschaften 
mit Einschluß ihrer Anwendungen. Bd. IV 1 I. 
Heft 4. XIIT + 804 8. Leipzig 1935, Verlag B. G. 
Teubner. Preis geh. 12.20 M, 

Dr.-Ing. CHR. BREUER, Führungsgzetriebe. 


Praktische Getriebeteehnik. hrsg. von Privatdozent 


Dr.-Inz. Kurt Rauh, Teehn. Hochseh., Heft 1.) 30 8. 
m. 128 Abb. Berlin 195. VDI-Verlag. Preis 
350 M. 

Dr.-Ine. A, REINHARDT, VDI. Dresden, Zum 
ÄAhnliehkeitsresetz für Hohlraum- 


bildluneen.und Dr.-Ing. H. KRANZ, VDI, Essen. 

StrömungeinSpiralgehäusen. (Forschungs- 

heft 370. Beilage zu „Forschung auf dem Gebiete 

Inzenieurwesens“, Ause. B. Bd. 6, Jan./Febr. 

31 S. m. 64 Abb. Berlin 1935. VDI-Verlag 
Preis 5 M. 
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1435.) 


m.b. H. 
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10. Deutscher Physiker- und Mathematikerta 
zu Pyrmont 


10. bis 15. September 1934 


J- 


Eröffnungsansprache des Vorsitzenden Dr. Karl Mey 


Is ıst mir eine große Ehre, die Tagung der Deutschen Physiker und Mathematiker eröffnen 
‚u dürfen. Ich heiße Sıe alle herzlich wıllkommen, insbesondere die Mitglieder der veranstaltenden 
(sesellschaften, der Deutschen Physikalischen Gesellschaft, der Deutschen Gesellschaft für technische 
Physik, der Gesellschaft für angewandte Mathematik und Mechanik und der Deutschen Mathematiker- 
Vereinigung. Ihnen und unsern Gästen wünsche ıch einen guten Verlauf der Tagung. Die große Zahl 
er Besucher, mehr als 300, zu denen ın einigen Tagen noch die Geophysiker hinzutreten werden, 
Ibertrifft unsere Vorschätzung erheblich und zeigt, daß diese seit Jahren hhebgewordene Einriehtung der 
Phvsiker- und Mathematikertagung noch nichts von ıhrer Zugkraft verloren hat und auch beı äußereı 
Ungunst der Zeiten zahlreiche Genossen ımmer wıeder zusammenführt zum Austausch ıhrer Gedanken 
ınd zur Klärung der Fortschritte, die ım abgelaufenen Jahre ın unseren Wissenschaften erreicht worden 
ınd. Ich freue mich, daß wir auch wıeder Vertreter unserer ausländischen Ortsgruppen hier begrüßen 
können, insbesondere aus Österreich und der Tschechoslowakei und schließe in diese Begrüßung auch ein 
unsere Freunde aus Holland und den andern Ländern, die mit großer Regelmäßigkeit hier erscheinen und 
über deren Ausbleiben wır uns wundern würden, wenn es einmal stattfinden sollte. Ich danke insbesondere 
en Herren aus Holland für die persönliche Mitwirkung durch ıhre Vorträge. Unsere Gesellschaften haben 
a viele Mitglieder im Auslande, bei der Physikalischen Gesellschaft beträgt diese Zahl sogar ein Drittel 
es Bestandes, und dıe Pflege der wissenschaftlichen und persönlichen Beziehungen zu diesen ausländischen 
lachgenossen ist immer eine Hauptsorge, aber auch Hauptfreude unserer Gesellschaften gewesen. Sie Ist 


n der Gegenwart besonders wichtig und wird auch von den höchsten Stellen unserer Reichsregierung 


anerkannt und gefördert. Wenn heute wegen mancherlei Umständen unsere auswärtigen Mitglieder nicht 


‚o zahlreich vertreten sind, wie ın andern Jahren, so mögen die erschienenen Mitglieder doch sicher sein, 
laß wır sie um so herzlicher willkommen heißen, und ich bitte sie, unsere Grüße und Gefühle denen zu 
(Ibermitteln, die nicht haben teilnehmen können, die wir aber vielleicht im nächsten Jahre um so zahl- 
eicher in unserer Mitte sehen werden. Ich bin überzeugt, daß die Eindrücke, die Sie von dieser Tagung 
ıitnehmen werden, durchaus freundlich sind und daß Sie davon gern in der Heimat berichten werden. 

ls ıst mır ebenfalls eine Ehre, die Herren des Bades Pyrmont, den Herrn Bürgermeister und den Herrn 
Nurdirektor hier begrüßen zu dürfen. Wir haben mit großer Freude die Einladung Pyrmont’s angenommen 


ınd der erste Eindruck hier ist für alle. die dieses herrliche Bad noch nicht kannten. überraschend sewesen. 





Ortes, die herrlichen Gärten und Parkanlagen und die prachtvollen Gebäude und Räume, 
unsere Veranstaltung zur Verlürune stellen. fassen unsere Tagune ın einen schönen äußeren 
ich danke Ihnen herzlich ım voraus für alles Entzegenkommen und für Ihre Bemühungen 

ls ist Übung geworden, bei wissenschaftlichen Tagungen ein gut Teil der Zeit auf Erörterungen übe: 
ie Stellung der betreffenden Wissenschaft ım neuen Staat. auf den Nachweis ıhrer Bedeutung und Volks 
verbundenheit zu verwenden: wir möchten aber in unserem Kreis bei der alten Übung bleiben und di: 
liche und fachliche Arbeıt für uns sprechen lassen. Daß wır Physiker und Mathematiker alle unserm 
Führer und Reıchskanzler folgen und ıhn ın seiner harten Arbeit und Mühe. ın seinen hohen Zielen für 
Deutschlands Wohl mit aller Kraft, jeder an seinem Platz, unterstützen, ıst selbstverständlich. Wir haben 
dieses Gelöbnis ım vorıgen Jahre ın Würzburg abgelegt und erneuern es ın dieser Stunde. Die Physik ıst 
akte Wissenschaft, sıe legt es ihren Anhängern nahe und ermöglicht es ıhnen, besser als durch Wort: 
oeleistete Arbeit und durch die sıchtbaren Fortschritte zu beweisen, daß sıe dıe richtige Eın 
zum Staatseanzen und zum Volkswohl haben 
Ich hofte, daß Ihnen schon die Wahl der beiden Hauptgruppen unserer Verhandlungsgegenstände das 
zeiven wird. indem neben der grundsätzlichen Forschung auch die Praxis und ıhre Gegenwarts 
aufgaben zu ıhrem Rechte kommen werden. Die „Physik der tiefen Temperaturen‘ wird ın einer Reıh: 
von Vorträgen behandelt werden. die neue Beobachtungen und Aufschlüsse für die Kenntnis der Natuı 
lıetern. \Wır freuen uns ganz besonders, daß hierbei mehrere Vorführungen höchster Experimentierkunst 
stattfinden werden und dabeı auch der schwere Wasserstoff. vor wenigen Jahren noch gänzlich unbekannt. 
Ihnen sogar ın flüssıger und fester Form vorgeführt werden wird. 

Das andere Hauptthema „Physik und Werkstoft' greift ıns volle Leben und wırd über viele techniscl 
wichtige Fortschritte berichten. Mögen Fernerstehende vielleicht etwas enttäuscht sein. wenn diese Bericht: 
‘ade unmittelbar zur Abfassung von aufsehenerregenden Zeitungsartikeln eignen werden: 
Phvsiker weiß, was es bedeutet, wenn Materialkonstanten um so und so viel haben gesteigert werden 
können, seı | Härte oder Zähiekeit eine S Werkstoftes. die elektrische stiokeit 1. oder wenn ın deı 
Auffassung des metallischen oder Glaszustandes oder ın der Theorie und Praxıs der Lichtstrahlung von 
Gasentladungen oder ın der Feststellung feinster Beimengungen bei Werkstoffen um nur einiges zu 
nennen ein sicherer Schritt voran getan worden ıst. Man möge hierbei bedenken. daß vieles. geradı 
weil es auf einem Gebiet heet. das gegenwärtige nıcht nur wissenschaftlich sondern auch volkswirtschaftlich 
wichtige ist. sich in diesem Augenblick noch nicht zur Veröffentlichung eignet. Es fallen darunter besonders 
viele Forschungen, die in den Laboratorien der Großindustrie nach jahrelanger Arbeit vor der Vollendung 
Ste hen und uns aul dem Ge bıet der \Werkstoffe. Fasern und vieler anderer (sevenstände des allgemeinen 
sroßen Verbrauchs mehr auf eigene Füße stellen sollen. und die dadurch für unsere gesamte Volkswirt 
schaft in den nächsten Jahren von großer Bedeutung sein werden. Hier sind außerordentliche Fortschritt: 

ım Entstehen. über die wir vielleicht auf den künftigen Tagungen ausführlich berichten können. 
Ob der Phvsiker zur Rohstofffrage viel beitragen könne, ob das nıcht mehr Sache des Chemikers sei, 


wird häufig gefragt. Gestatten Sie mir dazu aus meinem Arbeitsgebiet ein kleines Beispiel, um auch den 


Fernerstehenden einen Begriff zu geben, wie in jedem Fach und beı jedem Gegenstand diese Überlegung 


angestellt und zum Erfolge geführt werden kann. Ich meıne die Ihnen allen gut bekannten Rundfunk 
röhren. In ihnen ist Nickel in Form von Blechen und Drähten enthalten, wofür jährlich ungefähr einund 
dreiviertel Millionen Mark ins Ausland wandern. ein überraschend hoher Betrag. Für das Nickel unter den 
vielen bekannten Metallen und Metallesierungen einen Ersatzstoff zu finden, scheint leicht. ist es aber 
nicht. denn es wird gefordert, daß das Metall geschmeidig und gut formbar und dech auch ın gewissem 
Grade steif ist. daß es an der Luft nicht anläuft, daß es hohe Temperaturen verträgt ohne brüchig zu werden 
oder zu zerstäuben und vor allem, daß es leicht zu entgasen und auch noch billig ıst. Gelingt es nicht. all 
diese Eigenschaften und noch einige andere bei dem Ersatzmetall voll zu erreichen, so bleibt noch Änderung 
des gesamten Aufbaues der Radioröhren möglich, um zum Ziel zu gelangen, und man sieht leicht, daß 


hier dem Phvsiker eine auch vom Gesichtspunkt der Volkswirtschaft aus lohnende Aufgabe gestellt ist. 





mühungen 





1 


uneen uber 
und Volks 


en und di 


ılle unserm 


Zielen für 
Wir haben 
Physik ıst 
urch Wort« 


chtige Eın 


ıstände das 
‚evenwarts 
ıner Reıhi 
der Natuı 
entierkunst 


unbekannt. 


e technise| 
se Bericht: 
:n werden: 
ert werden 
venn ın deı 
ıhlung von 
einiges zu 
les. gerad: 
"tschaftlich 
" besonders 
Vollendung 
IIgemeınen 
Volkswirt 
'ortschritt: 
en. 

mikers sei. 
n auch den 
"berlegung 
Rundfunk 
‚hr einund 
unter den 
ist es aber 
1 GEWISSEM 
zu werden 
nicht, all 
ı Änderung 
leicht, daß 
sestellt ıst. 


nd Räume. 


en aubßberen 


wie ın diesem Beispiel wird die Aufgabe sıch für den Physiker ın sehr vielen andern Fällen darstellen 

ist im allgemeinen nicht wahrscheinlich, daß für jeden der Stoffe, deren Ersetzung durch einheimisch: 
‚ünschenswert ist. ein ın allen Punkten entsprechender Ersatz gefunden werden kann. Dann ıst unteı 
mständen schon gleichzeitige mit der Entwicklung des neuen Stoffes dıe weıtere Aufgabe zu lösen, dıe 
Konstruktionen, Verwendungen oder Herstellungsverfahren der behandelten Gegenstände so umzuformen, 
aß auch mit dem neuen Werkstoff mindestens die gleichen Leistungen der Gegenstände erreicht 
erden. Ich will darüber nichts weiteres sagen; die Industrie und ıhre Physiker kennen ıhre Aufgaben 
nd sınd am Werk. Manches wird ın den nachfolgenden Vorträgen noch näher auseinandergesetzt 


( rden. 


Ich könnte viellei ht mit dem. was ıch vesagt habe. den Kındrus k erwecken, als ob oevenüber dei 
\nwendung der Physik ın der Technik und ıhrer gegenwärtigen großen nationalen Bedeutung die Leistun: 
es theoretisch oder experimentell rein grundsätzlich forschenden Phvsikers geringer als früher einzuschätzen 
äre. Das ıst nıcht so. Wie gerade unsere Tagungen zeigen, arbeiten beide Teile Hand ın Hand und müssen 
‚ arbeiten. denn gerade die Lösung der bedeutendsten technischen Aufgaben führt immer wıeder zurück 

vielen Grundfragen unserer Wissenschaft über den Aufbau der Stoffe, dıe Natur des metallischen Zu 
tandes,. der Atome und Elektronen, über die Theorie günstiger Kreisprozesse und anderes. Ich glaub: 
‚war sagen zu dürfen, daß die theoretische Forschung jetzt auch ın den großen Laboratorien der Industrie, 
ie sich immer gewaltiger entfalten und bedeutungsvolle allgemeine Ergebnisse erzielen, eine Heimstätt: 
efunden hat. Die verantwortlichen Leiter haben dort erkannt. daß der äußerste Fortschritt sıch nicht 
ur durch Experimentieren erreichen läßt, sondern daß die Gedanken, die dabeı führen, die höchste Kennt 
nıs und Vorstellung über das, was hinter den Erscheinungen liegt, zur Grundlage haben müssen. Zugleich 
ıber zeigen auch dıe Herren an den Hochschulen immer mehr Neigung, weitgehend und gern den Änrezungen 
ınd Sorgen der andern Seite zu folgen, und es scheint sıch eine ım Laufe der Jahre ımmer besser und 
evensreicher werdende Zusammenarbeit zu entwickeln. Dies wırd erleichtert durch die gegenüber älteren 
/.«eiten viel nähere Berührung mit technischen Dingen und Fragen. der sıch auch die theoretisch veranlagten 
üngeren Physiker nıcht haben entziehen können. Rang- und ähnliche Streitigkeiten sind im Schwinden 
ınd haben keinen Platz mehr, die Zeit ıst dazu zu ernst. Die Physik ıst ja eine von den Wissenschaften 
lie ın ihren Forschungen besonders günstig gestellt sınd. Sıe stellt ihre Fragen an dıe Natur und bekomm! 
ie durch das Mittel des Experimentes beantwortet. Die Art der Fragestellung hängt freilich weıtgehen(! 
von den theoretischen Vorstellungen ab. Sınd diese abwegig. so pflegt sıch das nach nıcht sehr langeı 
Zeit durch ıhre unmittelbare Widerlegung aus dem Experiment oder durch Unfruchtbarkeit zu neueı 
"ragestellung an dıe Natur zu oftenbaren. Diese Reinigungsprozesse sind durchaus natürlich und würden 
sich auch ruhig abspielen, wenn nicht aus Beweggründen, die gänzlich außerhalb physikalischer Gedanken 
ange liegen, sie öfter zum Tummelplatz eigensüchtiger Bestrebungen gemacht würden. 

Hs ıst deshalb hocherfreulich. daß der Führer ın seiner großen Rede kürzlieh aueh ın Kulturdiıngen 
entschiedene Klarheit gefordert und mit nicht mißzuverstehender Eindeutigkeit erklärt hat. daß dazu 
vor allem gehöre die Unabhängigkeit und Reinheit der wissenschaftlichen Forschung und ıhre Befreiung 
von Beeinflussungen und Beweggründen, die nicht auf dem Drang nach Wahrheit und ıhrer Anwendung 
‚um allgemeinen Wohle beruhen. Das schließt ein. daß künftig persönliche und nieht dem allgemeinen 
Wohle dienende Bestrebungen und Gefühle zurückgedrängt werden müssen, um dıe großen Aufgaben ın 
voller Klarheit durchführen zu können. Der Führer hat ın dıeser Rede auch von der Zahl und 
Kühnheit der Erfindungen und Entdeckungen gesprochen, denen dıe Jahrtausende 
vorher nichts Vergleichbares zur Seite stellen könnten. Er hat dıe Physik als ecıne Grob 
macht bezeichnet, „die die nicht wissenschaftlichen Ansıchten,. vermeintlichen Erkennt 
nisse und Vorurteile zerbricht. die Hand in Hand mit der Technik und Chemie dıe Welt 
kenntnis dauernd weitet und die Schätze des Erdballs mobilisiert für einen Aufstieg der 
\lenschheit, der in seinem Tempo fast beängstigend wirkt‘. Kein von seinem Fach begeisterter 


Physiker könnte eine höhere und schönere Auffassung von seiner Wissenschaft aussprechen, als es unser 





Führer und Reichskanzler mit diesen Worten vetan hat. Und doch hört man nicht selten. daß das Zeitalter 
dieser exakten Wissenschaft vorüber seı. daß Biologie. Eugenik und andere Zweige dıe Führung übernehmen 
würden. Neiıdlos werden wır dem Aufblühen neuer Zweige der Forschung zusehen und freudız mithelfen. 
wo wir können; aber wer ım Leben der Physik steht und sıeht, welche Fülle neuer Erkenntnisse und Auf- 
saben sich von Jahr zu Jahr ergibt, und zwar von grundlegender Art ich erinnere nur an das Auffinden 
sanz neuer Stoffe oder der Höhenstrahlung, dıe uns ın ungeahnter Weise mit den fernsten Welten ver- 
bindet kann kaum der Ansıcht seın, wenn er nıcht selbst schon müde ıst, daß das Wesentliche ın der 
Physik getan sei. Diese herrliche Wissenschaft steht erst am Anfang und ıst noch jung und kräftıg genug, 
um uns alle noch Gewaltiges miterarbeiten und miterleben zu lassen. In den Worten des Führers. der die 


'hysik als eine Großmacht bezeichnet hat, legt eine hohe Ehrung für uns. Ich möchte nur wünschen, 


daß alle. die diese Großmacht vertreten. sıch auch der ıhnen damit auferleeten Verpflichtung ımmer bewußt 


sein mögen. Unsere physikalischen und mathematischen Gesellschaften werden es tun und ın sachlicher 
Arbeit, unter Einsatz aller Kräfte, ıhrer Wissenschaft dienen und darauf auch die Hoffnung und das Recht 
auf weiteres Blühen gründen. Sıe werden auch weıter den Zusammenhang mit den Wissenschaften und 
(Gelehrten in andern Ländern pflegen, wie es Ja von unserem Führer gefordert worden ist. Wir werden uns 
treuen und stolz sein, wenn aus unserer Arbeit sichtbare Ergebnisse entspringen und diese zum Wohle 
unseres Volkes dienen. Wir sınd ın diesem Willen einıg und wollen das, darum bitte ich Sıe, bekräftigen, 
indem wir auf unser hebes deutsches Vaterland und auf den Reıchskanzler und Führer, der ın harter Arbeit 


und Sorsoe uns ın eine schöne Zukunft zu führen sucht. eın dreifaches Sıee-Heil ausbrineen. 



























ıs Zeitalter 
‚ernehmen 
mithelfen. 
und Auf 

Auftinden 
‚elten ver 

che ın der 
Itıe genug, 
0, der dıe 
wünschen, 
ier bewußt 
sachlieher 
das Recht 
ratten und 
verden uns 
um Wohle 
ekräftigen, 


rter Arbeit 





ZEITSCHRIFT FÜR 
ANGEWANDTE MATHEMATIK 


92| 39 
UNIVERSITY MICROFILMS 





